Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,167,912 +A coding (725/1344 nt) oprB‑I → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,167,9101.A91.8% 1420.9 / 114.6 439coding (723/1344 nt)oprB‑Icarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base . (19/17);  new base A (198/205);  total (217/222)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 443 total aligned reads displayed.

ATCTGGTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCAGGCGCTGACA  >  NC_002947/1167764‑1168054
                                                                                                                                                   |                                                                                                                                                 
 tCTGGTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCaaag                                                                                                                                              >  1:4772396/1‑148 (MQ=255)
  cTGGTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAagg                                                                                                                                             <  1:9617671/150‑1 (MQ=255)
   tGGTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAaggg                                                                                                                                            >  1:6166949/1‑150 (MQ=255)
    ggTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAagggg                                                                                                                                           <  1:5489524/150‑1 (MQ=255)
     gTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAaggggg                                                                                                                                          <  2:8979950/150‑2 (MQ=255)
     gTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAaggggg                                                                                                                                          >  1:10035466/1‑149 (MQ=255)
      tACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAagggggc                                                                                                                                         >  1:4973544/1‑148 (MQ=255)
         aaCTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGAt                                                                                                                                       <  2:2424467/150‑1 (MQ=255)
         aaCTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGAt                                                                                                                                       >  1:4076743/1‑150 (MQ=255)
          aCTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATc                                                                                                                                      >  1:1188566/1‑150 (MQ=255)
              gCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGc                                                                                                                                  <  2:9148428/150‑1 (MQ=255)
              gCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGc                                                                                                                                  >  1:1070783/1‑150 (MQ=255)
                tgaccggCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCg                                                                                                                                <  2:7605950/144‑1 (MQ=255)
                cGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCg                                                                                                                                >  1:7244661/1‑150 (MQ=255)
                  gTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGt                                                                                                                              <  1:3602345/150‑1 (MQ=255)
                  accggCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGt                                                                                                                              <  2:3709018/146‑1 (MQ=255)
                     aGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACTGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGa                                                                                                                           >  1:2762514/1‑150 (MQ=255)
                          gTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCgg                                                                                                                      >  2:7118458/1‑150 (MQ=255)
                           tGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGgt                                                                                                                     <  2:8484766/150‑1 (MQ=255)
                           ggggCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGgt                                                                                                                     <  2:8118878/149‑1 (MQ=255)
                             ggCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGgtgt                                                                                                                   >  1:8549474/1‑150 (MQ=255)
                              gcgcTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGtgtg                                                                                                                  >  1:2991732/1‑150 (MQ=255)
                                 cTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCGTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTc                                                                                                               <  2:306142/150‑1 (MQ=255)
                                     gggTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGccc                                                                                                           >  1:7430178/1‑150 (MQ=255)
                                      ggTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCa                                                                                                          >  2:5681297/1‑150 (MQ=255)
                                       gTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCaa                                                                                                         <  2:11330007/150‑1 (MQ=255)
                                        tGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAg                                                                                                        >  1:8250502/1‑150 (MQ=255)
                                        tGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAg                                                                                                        >  1:11300966/1‑150 (MQ=255)
                                         gAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAgg                                                                                                       >  2:11070360/1‑150 (MQ=255)
                                            gTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCa                                                                                                    >  2:6711688/1‑150 (MQ=255)
                                                 acaTCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGc                                                                                               <  1:6260483/150‑1 (MQ=255)
                                                   ctCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCt                                                                                             <  2:11303028/149‑1 (MQ=255)
                                                   aTCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGATCAGAACCCTTCCAACCTGGAAACCGGCAACGGATTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCAGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCt                                                                                             <  1:5268562/150‑1 (MQ=255)
                                                    tCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTg                                                                                            >  1:6034022/1‑150 (MQ=255)
                                                     cACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGc                                                                                           <  1:1713084/150‑1 (MQ=255)
                                                            gAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGa                                                                                    <  1:1485234/150‑1 (MQ=255)
                                                             aGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAg                                                                                   <  2:4973181/150‑1 (MQ=255)
                                                                 cttcGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGAACTTGCCGGTGGACGCGGTGTGTTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTAcc                                                                               >  2:2229313/1‑150 (MQ=38)
                                                                  ttcGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCg                                                                              >  1:11102518/1‑150 (MQ=255)
                                                                     gTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCt                                                                           >  2:1701780/1‑150 (MQ=255)
                                                                       aCAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTgg                                                                         >  2:10874195/1‑150 (MQ=255)
                                                                       aCAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTgg                                                                         >  1:4508713/1‑150 (MQ=255)
                                                                       aCAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTgg                                                                         <  2:11496511/150‑1 (MQ=255)
                                                                         aGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATAGCCTGCCGGGCGAGTACCGCCTGGGt                                                                       >  1:2133264/1‑150 (MQ=255)
                                                                           gTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTta                                                                     <  1:637234/150‑1 (MQ=255)
                                                                           gTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTta                                                                     <  2:2788673/150‑1 (MQ=255)
                                                                              ggggCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtacta                                                                  <  2:1294074/150‑1 (MQ=255)
                                                                              ggggCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtacta                                                                  >  2:2741351/1‑150 (MQ=255)
                                                                                 gCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtactacta                                                               <  1:9844846/150‑1 (MQ=255)
                                                                                 gCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtactacta                                                               <  1:3994640/150‑1 (MQ=255)
                                                                                  ccTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCACAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtactactac                                                              <  1:11537794/150‑1 (MQ=255)
                                                                                   cTTCGAGCAGAACCCTTCCAACCTTGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACa                                                             <  1:735491/150‑1 (MQ=255)
                                                                                   cTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACa                                                             <  2:4508713/150‑1 (MQ=255)
                                                                                   cTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACa                                                             <  1:10454478/150‑1 (MQ=255)
                                                                                    ttCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAg                                                            >  2:11610320/1‑150 (MQ=255)
                                                                                     tCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGc                                                           >  1:1846574/1‑150 (MQ=255)
                                                                                       gAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGt                                                                                                                              >  2:6273984/1‑81 (MQ=255)
                                                                                       gAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTAAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGctc                                                         >  2:10732187/1‑148 (MQ=255)
                                                                                         gCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGc                                                                             <  1:10180740/128‑1 (MQ=255)
                                                                                         gCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGc                                                                             >  2:10180740/1‑128 (MQ=255)
                                                                                         gCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACgg                                                       >  1:6590383/1‑150 (MQ=255)
                                                                                         gCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAACACgg                                                       <  1:5681297/150‑1 (MQ=255)
                                                                                             aaCCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCaa                                                   >  1:2252801/1‑150 (MQ=255)
                                                                                             aaCCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCaa                                                   <  2:3241315/150‑1 (MQ=255)
                                                                                             aaCCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCaa                                                   <  2:4170049/150‑1 (MQ=255)
                                                                                              aCCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtacta                                                                  <  1:9443904/134‑1 (MQ=255)
                                                                                              aCCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTtacta                                                                  >  2:9443904/1‑134 (MQ=255)
                                                                                               cccTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCg                                                                                                                       <  2:10588941/80‑1 (MQ=255)
                                                                                               cccTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCg                                                                                                                       >  1:10588941/1‑80 (MQ=255)
                                                                                                ccTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGc                                                <  2:6541084/150‑1 (MQ=255)
                                                                                                ccTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGc                                                <  1:2683813/150‑1 (MQ=255)
                                                                                                ccTTCCAACCTGGAAACCGGCAACGGCTTCAAGCCCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGc                                                <  1:1211758/150‑1 (MQ=255)
                                                                                                 cTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCt                                               <  2:9792691/150‑1 (MQ=255)
                                                                                                  ttCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGGAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTg                                              <  2:4772396/150‑1 (MQ=255)
                                                                                                  ttCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGc                                                                                      >  1:4700045/1‑110 (MQ=255)
                                                                                                  ttCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTggg                                                                        >  1:7542860/1‑124 (MQ=255)
                                                                                                  ttCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTg                                              >  1:1553775/1‑150 (MQ=255)
                                                                                                   tCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGa                                             <  2:8755258/150‑1 (MQ=255)
                                                                                                    ccAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGAc                                            <  2:7871123/150‑1 (MQ=255)
                                                                                                      aaCCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGa                                          >  2:5981684/1‑150 (MQ=255)
                                                                                                      aaCCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGa                                          >  1:7743678/1‑150 (MQ=255)
                                                                                                       aCCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGAt                                         <  1:3325267/150‑1 (MQ=255)
                                                                                                       aCCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGAt                                         >  1:11485200/1‑150 (MQ=255)
                                                                                                       aCCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGAt                                         <  1:10870180/150‑1 (MQ=255)
                                                                                                           ggAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATgtgt                                     <  2:3948657/150‑1 (MQ=255)
                                                                                                           ggAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATgtgt                                     <  1:2823903/150‑1 (MQ=255)
                                                                                                              aaCCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTacg                                  <  1:6747781/150‑1 (MQ=255)
                                                                                                                 cGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTacgacga                              >  2:6322720/1‑150 (MQ=255)
                                                                                                                  ggCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTacgacga                              <  2:7430178/150‑1 (MQ=255)
                                                                                                                     aaCGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTTGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGt                           >  1:9255416/1‑150 (MQ=255)
                                                                                                                     aaCGGCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGt                           >  2:1421136/1‑150 (MQ=255)
                                                                                                                        ggCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCaa                        >  1:6920021/1‑150 (MQ=255)
                                                                                                                        ggCTTCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCaa                        >  2:2072457/1‑150 (MQ=255)
                                                                                                                           ttCAAGCTCAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACgg                     >  2:2558022/1‑150 (MQ=255)
                                                                                                                                   cAGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGc             <  1:8845327/150‑1 (MQ=255)
                                                                                                                                    aGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCa            <  1:2423134/150‑1 (MQ=255)
                                                                                                                                    aGCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCa            >  2:3548054/1‑150 (MQ=255)
                                                                                                                                     gCGGCAGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCAg           >  2:2825574/1‑150 (MQ=255)
                                                                                                                                          aGTGGTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCAGGCGCt      <  2:6773891/150‑1 (MQ=255)
                                                                                                                                             ggTACCAAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCAGGCGCTGAc   <  1:1356878/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCTGGTACAACTGGCCGGTCAGCCAGTGGGCGCTGCGGGTGAAGTACAACATCACGCCGGAGTTCTTCGTACAGGTCGGGGCCTTCGAGCAGAACCCTTCCAACCTGGAAACCGGCAACGGCTTCAAGCTCAGCGGCAGTGGTACCAAGGGGGCGATCTTGCCGGTGGAAGCGGTGTGGTCGCCCAAGGTCAATGGCCTGCCGGGCGAGTACCGCCTGGGTTACTACTACAGCACGGCCAAGGCTGACGATGTGTACGACGACGTCAACGGCAACCCGCAGGCGCTGACA  >  NC_002947/1167764‑1168054

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: