Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,673,133 C→T 5.1% W6* (TGG→TAG)  mdtI ← multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,673,1330CT 5.1% 537.2 / 12.0 155W6* (TGG→TAG) mdtImultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base C (64/83);  new base T (4/4);  total (68/87)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 159 total aligned reads displayed.

CTTGAGAAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACA  >  NC_000913/1672987‑1673281
                                                                                                                                                  |                                                                                                                                                    
cTTGAGAAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCAtt                                                                                                                                                   <  2:3053642/150‑1 (MQ=255)
cTTGAGAAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCAtt                                                                                                                                                   <  2:1711561/150‑1 (MQ=255)
    agaAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCaaa                                                                                                                                               <  2:824894/150‑1 (MQ=255)
    agaAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCaaa                                                                                                                                               <  2:2098128/150‑1 (MQ=255)
     gaAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAAc                                                                                                                                              <  1:1880538/150‑1 (MQ=255)
       aaGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATCTTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTg                                                                                                                                            <  2:3237473/150‑1 (MQ=255)
       aaGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTg                                                                                                                                            <  2:3949145/150‑1 (MQ=255)
         gcgcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTgcg                                                                                                                                          <  2:3602439/150‑1 (MQ=255)
         gcgcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTgcg                                                                                                                                          >  1:209392/1‑150 (MQ=255)
          cgcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTgcgc                                                                                                                                         >  1:1037119/1‑150 (MQ=255)
          cgcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTgcgc                                                                                                                                         >  1:2943532/1‑150 (MQ=255)
          cgcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTACAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTgcgc                                                                                                                                         <  2:2679255/150‑1 (MQ=255)
           gcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGcc                                                                                                                                        >  1:3869519/1‑150 (MQ=255)
           gcACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGcc                                                                                                                                        >  2:2266222/1‑150 (MQ=255)
             aCTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCAt                                                                                                                                      >  2:3028303/1‑150 (MQ=255)
             aCTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCAt                                                                                                                                      >  2:2477372/1‑150 (MQ=255)
               tAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATgg                                                                                                                                    >  1:2498893/1‑150 (MQ=255)
               tAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATgg                                                                                                                                    >  1:2800036/1‑150 (MQ=255)
               gacaGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTGTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATgg                                                                                                                                    <  1:2400450/147‑1 (MQ=255)
                  aGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTc                                                                                                                                 <  2:298547/150‑1 (MQ=255)
                  aGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTc                                                                                                                                 >  2:1076328/1‑150 (MQ=255)
                  aGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTc                                                                                                                                 <  2:1106951/150‑1 (MQ=255)
                  aGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTc                                                                                                                                 <  2:3685703/150‑1 (MQ=255)
                   ggCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCa                                                                                                                                >  1:193069/1‑150 (MQ=255)
                     cagccagcACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCAcc                                                                                                                              >  1:3119261/1‑150 (MQ=255)
                       gccagcACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTc                                                                                                                            <  2:3907220/150‑1 (MQ=255)
                         cagcACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCa                                                                                                                          <  2:2176879/150‑1 (MQ=255)
                          agcACCGCCGCCAGTGAGAGCAAGCCAAATATATTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAg                                                                                                                         <  1:3226942/150‑1 (MQ=255)
                              ccgccgCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCa                                                                                                                     <  1:2495269/150‑1 (MQ=255)
                              ccgccgCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCa                                                                                                                     <  1:3561359/150‑1 (MQ=255)
                              ccgccgCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCa                                                                                                                     <  1:2772150/150‑1 (MQ=255)
                               cgccgcCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAg                                                                                                                    <  1:2976108/150‑1 (MQ=255)
                                gccgccAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAgg                                                                                                                   <  1:3223918/150‑1 (MQ=255)
                                     cAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTAc                                                                                                              <  1:147525/150‑1 (MQ=255)
                                     cAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTAc                                                                                                              >  1:148814/1‑150 (MQ=255)
                                              cAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTAc                                                                                                     >  1:1149038/1‑150 (MQ=255)
                                              cAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTAc                                                                                                     >  1:367390/1‑150 (MQ=255)
                                                aGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACgg                                                                                                   <  2:3751027/150‑1 (MQ=255)
                                                   cAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCGGGTTTACGCGCTTTACGGGTa                                                                                                >  1:423124/1‑150 (MQ=255)
                                                          tttGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGAtt                                                                                         >  2:1159330/1‑150 (MQ=255)
                                                               gacgaAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATc                                                                                    >  1:3834310/1‑150 (MQ=255)
                                                                 cgaAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCaa                                                                                  >  2:3598869/1‑150 (MQ=255)
                                                                 cgaAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCaa                                                                                  <  2:2101260/150‑1 (MQ=255)
                                                                 cgaAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCaa                                                                                  >  1:193279/1‑150 (MQ=255)
                                                                   aaaGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAaca                                                                                <  1:1609873/150‑1 (MQ=255)
                                                                    aaGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCCTTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAacac                                                                               <  1:3097945/150‑1 (MQ=255)
                                                                     aGCCGTCAGAAAATTTCAAAAAGACGTTATCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAAcaca                                                                              <  1:871238/150‑1 (MQ=255)
                                                                       ccGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCTATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAAt                                                                            >  1:3273923/1‑150 (MQ=255)
                                                                       ccGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCTATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAAt                                                                            >  1:1302361/1‑150 (MQ=255)
                                                                       ccGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCTATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAAt                                                                            >  1:3273927/1‑150 (MQ=255)
                                                                        cGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATc                                                                           >  1:2940019/1‑150 (MQ=255)
                                                                           cAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCg                                                                        <  1:957522/150‑1 (MQ=255)
                                                                           cAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCg                                                                        <  1:3951437/150‑1 (MQ=255)
                                                                            aGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCgg                                                                       >  1:3857161/1‑150 (MQ=255)
                                                                             gAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGc                                                                      <  1:1912802/150‑1 (MQ=255)
                                                                             gAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGc                                                                      <  2:2084190/150‑1 (MQ=255)
                                                                              aaaaTTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCa                                                                                                                          >  1:1398685/1‑97 (MQ=255)
                                                                              aaaaTTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCa                                                                                                                          <  2:1398685/97‑1 (MQ=255)
                                                                               aaaTTTCAAAAAGACGTGAGCAACGATTTCCAGCCCGATTGCCACTGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGa                                                                    <  2:378348/150‑1 (MQ=255)
                                                                                aaTTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGAc                                                                   >  2:2734346/1‑150 (MQ=255)
                                                                                   ttCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAg                                                                >  1:2113602/1‑150 (MQ=255)
                                                                                       aaaaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTCCCTGATTTTATCAACACAATCCCGGCGACCAGggtg                                                            >  2:2071309/1‑150 (MQ=255)
                                                                                       aaaaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGggtg                                                            >  2:2255427/1‑150 (MQ=255)
                                                                                        acaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGTTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGggtgg                                                           <  2:298793/148‑1 (MQ=255)
                                                                                            aCGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTaa                                                       >  1:2101147/1‑150 (MQ=255)
                                                                                                 aGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCgg                                                  <  2:2988030/150‑1 (MQ=255)
                                                                                                  gCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGc                                                 <  2:3119261/150‑1 (MQ=255)
                                                                                                    aaCGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCaa                                               >  1:1805787/1‑150 (MQ=255)
                                                                                                     aCGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAAt                                              <  2:3415000/150‑1 (MQ=255)
                                                                                                       gATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAAttt                                            <  1:3678666/150‑1 (MQ=255)
                                                                                                         tttCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtc                                          <  2:687376/150‑1 (MQ=255)
                                                                                                             cAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtcatca                                      <  2:236122/150‑1 (MQ=255)
                                                                                                              aGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAg                                     >  1:3917763/1‑150 (MQ=255)
                                                                                                               gCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGc                                    >  1:3318416/1‑150 (MQ=255)
                                                                                                               gCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGc                                    >  2:3329341/1‑150 (MQ=255)
                                                                                                                    agtGCCAATGCCAGCCAGGCGGCGTGAACCCATTCACACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCTACACAATCCCGGCGACCAGGGTGGTTATCCCGGCTATTTTCATCAGCGATaa                               <  2:2358322/148‑1 (MQ=255)
                                                                                                                     ttGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATaaa                              <  1:2235930/150‑1 (MQ=255)
                                                                                                                     ttGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATaaa                              >  2:3479963/1‑150 (MQ=255)
                                                                                                                        ccAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACtt                           >  1:104214/1‑150 (MQ=255)
                                                                                                                        ccAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACtt                           >  2:1435103/1‑150 (MQ=255)
                                                                                                                         cAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACttt                          >  1:442879/1‑150 (MQ=255)
                                                                                                                         cAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACttt                          >  2:656491/1‑150 (MQ=255)
                                                                                                                          aaTGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTtc                         <  2:3883382/150‑1 (MQ=255)
                                                                                                                           aTGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTtcg                        >  1:857609/1‑150 (MQ=255)
                                                                                                                             gccagccaGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTtcgtc                      <  1:780227/150‑1 (MQ=255)
                                                                                                                              ccagccagGCGGCGTGACCCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTtcgtcg                     >  1:992261/1‑150 (MQ=255)
                                                                                                                                agccagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGaa                   >  1:716046/1‑150 (MQ=255)
                                                                                                                                 gccagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAAt                  <  2:3111775/150‑1 (MQ=255)
                                                                                                                                 gccagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAAt                  <  2:1718112/150‑1 (MQ=255)
                                                                                                                                    agGCGGCGTGAACCTATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAAc               <  2:1302356/150‑1 (MQ=255)
                                                                                                                                    agGCGGCGAGAACCTATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAAc               <  2:3273923/150‑1 (MQ=255)
                                                                                                                                        ggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaaa           >  2:1511264/1‑150 (MQ=255)
                                                                                                                                        ggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaaa           >  2:1511258/1‑150 (MQ=255)
                                                                                                                                        ggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaaa           <  1:3466895/150‑1 (MQ=255)
                                                                                                                                        ggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaaa           <  1:2889680/150‑1 (MQ=255)
                                                                                                                                         gcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAAc          <  2:482367/150‑1 (MQ=255)
                                                                                                                                         gcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAAc          <  2:482753/150‑1 (MQ=255)
                                                                                                                                               aCCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTaaa    >  1:803619/1‑150 (MQ=255)
                                                                                                                                                cccATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAAc   >  2:1771184/1‑150 (MQ=255)
                                                                                                                                                cccATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAAc   >  1:1113313/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTGAGAAAGCGCACTAAAGGCAGCCAGCACCGCCGCCAGGGAGAGCAAGCCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACA  >  NC_000913/1672987‑1673281

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑