Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,580,194 T→C 8.1% intergenic (‑189/+21) pcaP ← / ← galP‑II porin/porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,580,1940TC 8.1% 832.2 / 39.0 270intergenic (‑189/+21)pcaP/galP‑IIporin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (123/125);  new base C (11/11);  total (134/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGCACCAC  >  NC_002947/1580045‑1580343
                                                                                                                                                     |                                                                                                                                                     
gATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTAt                                                                                                                                                       >  1:1077996/1‑150 (MQ=255)
gATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTAt                                                                                                                                                       >  2:4690084/1‑150 (MQ=255)
  tAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATcc                                                                                                                                                     <  2:3957081/150‑1 (MQ=255)
      gATCAGGAGAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGa                                                                                                                                                 <  2:3369061/150‑1 (MQ=255)
       aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGaa                                                                                                                                                <  1:1437912/150‑1 (MQ=255)
       aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGaa                                                                                                                                                <  1:4551258/150‑1 (MQ=255)
         cAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGaata                                                                                                                                              <  1:4831420/150‑3 (MQ=255)
         cAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAaga                                                                                                                                              >  2:752410/1‑150 (MQ=255)
          aGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAatag                                                                                                                                             <  2:5375111/150‑1 (MQ=255)
           ggATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAgg                                                                                                                                            <  1:808079/150‑1 (MQ=255)
            gATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                           <  2:1280579/150‑1 (MQ=255)
            gATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                           >  2:1423832/1‑150 (MQ=255)
             aTAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAGGcc                                                                                                                                          <  1:1593132/150‑1 (MQ=255)
                gCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAt                                                                                                                                       >  2:3301052/1‑150 (MQ=255)
                 ccTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAtt                                                                                                                                      <  1:1168688/150‑1 (MQ=255)
                   tGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                                    <  1:4410993/150‑1 (MQ=255)
                   tGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                                    <  1:1009813/150‑1 (MQ=255)
                     cAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGt                                                                                                                                  <  2:3615181/150‑1 (MQ=255)
                     cAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGt                                                                                                                                  <  2:1051804/150‑1 (MQ=255)
                         caccacCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcaca                                                                                                                              <  2:1572212/150‑1 (MQ=255)
                            caccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaa                                                                                                                           <  1:4482640/150‑1 (MQ=255)
                             accaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaag                                                                                                                          >  2:2184321/1‑150 (MQ=255)
                             accaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaag                                                                                                                          >  2:5078412/1‑150 (MQ=255)
                               catcatCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAt                                                                                                                                       <  1:684512/135‑1 (MQ=255)
                                   atcaGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAg                                                                                                                    <  2:4561117/150‑1 (MQ=255)
                                    tcagctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                                   >  2:2354164/1‑150 (MQ=255)
                                    tcagctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGc                                                                                                                   <  2:1133591/150‑1 (MQ=255)
                                         tcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAatag                                                                                                                                             <  2:5183550/119‑1 (MQ=38)
                                           agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAc                                                                                                            >  2:3480104/1‑150 (MQ=255)
                                           agcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAc                                                                                                            <  1:1509528/150‑1 (MQ=255)
                                             cGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTc                                                                                                          <  2:2862491/150‑1 (MQ=255)
                                              gAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcac                                                                                                                               <  2:4984641/128‑1 (MQ=255)
                                              gAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcac                                                                                                                               >  1:4984641/1‑128 (MQ=255)
                                                aTGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAGGCCagac                                                                                                                                      >  1:289565/1‑116 (MQ=37)
                                                aTGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTcaca                                                                                                       >  2:5109440/1‑150 (MQ=255)
                                                        aGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCg                                                                                               >  1:5419207/1‑150 (MQ=255)
                                                           tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGtt                                                                                            >  2:912159/1‑150 (MQ=255)
                                                           tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAACCGGtt                                                                                            >  1:520317/1‑150 (MQ=255)
                                                             gCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATaa                                                                                                             <  2:2287335/131‑1 (MQ=255)
                                                             gCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATaa                                                                                                             >  1:2287335/1‑131 (MQ=255)
                                                                 tGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATc                                                                                      <  1:3439748/150‑1 (MQ=255)
                                                                   ggCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGa                                                                                    >  1:5565803/1‑150 (MQ=255)
                                                                     cAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCGAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAAc                                                                                  >  1:3136043/1‑150 (MQ=255)
                                                                         tCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGt                                                                              <  1:2374386/150‑1 (MQ=255)
                                                                         tCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGt                                                                              >  1:3788488/1‑150 (MQ=255)
                                                                         tCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGt                                                                              >  1:2611643/1‑150 (MQ=255)
                                                                          cATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGtt                                                                             <  1:5595178/150‑1 (MQ=255)
                                                                           aTGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTg                                                                            <  1:4088505/150‑1 (MQ=255)
                                                                             ggCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAatag                                                                                                                                             <  1:2835817/83‑1 (MQ=25)
                                                                             ggCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGt                                                                          <  2:2815486/150‑1 (MQ=255)
                                                                                aGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCt                                                                       >  2:2664733/1‑150 (MQ=255)
                                                                                  ggTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGa                                                                     >  1:2572942/1‑150 (MQ=255)
                                                                                  ggTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGGGTTCTCATCGAACGCGTTGGTGCTGaa                                                                    >  1:373964/1‑150 (MQ=255)
                                                                                   gTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcac                                                                                                                               >  2:2928738/1‑91 (MQ=255)
                                                                                   gTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcac                                                                                                                               <  1:2928738/91‑1 (MQ=255)
                                                                                   gTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                                    <  2:2120542/150‑1 (MQ=255)
                                                                                   gTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                                    >  1:27079/1‑150 (MQ=255)
                                                                                    ttttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAg                                                                   >  1:2325939/1‑150 (MQ=255)
                                                                                     tttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGt                                                                  <  1:5057773/150‑1 (MQ=255)
                                                                                         atcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTcgcg                                                              >  2:2536629/1‑150 (MQ=255)
                                                                                          tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGa                                                             <  1:1725936/150‑1 (MQ=255)
                                                                                          tcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGa                                                             <  2:3000071/150‑1 (MQ=255)
                                                                                           caCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGAc                                                            <  2:3898273/150‑1 (MQ=255)
                                                                                                tGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATg                                                       <  2:4966196/150‑1 (MQ=255)
                                                                                                     ccAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGa                                                  <  1:4140334/150‑1 (MQ=255)
                                                                                                     ccAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGa                                                  >  2:4168992/1‑150 (MQ=255)
                                                                                                          gatgCCTTCAACAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                             <  2:5361010/150‑1 (MQ=255)
                                                                                                          gatgCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                             <  1:3275139/150‑1 (MQ=255)
                                                                                                          gatgCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTCgc                                             <  1:773599/150‑1 (MQ=255)
                                                                                                               cTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                        >  2:1966306/1‑150 (MQ=255)
                                                                                                               cTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                                        >  2:4850545/1‑150 (MQ=255)
                                                                                                                ttCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAtt                                       >  1:2097879/1‑150 (MQ=255)
                                                                                                                      agggagggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCg                                 >  1:4358133/1‑150 (MQ=255)
                                                                                                                         gagggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACa                              >  1:5132823/1‑150 (MQ=255)
                                                                                                                          agggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAg                             <  2:2615442/150‑1 (MQ=255)
                                                                                                                          agggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAg                             <  1:2052964/150‑1 (MQ=255)
                                                                                                                           gggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAgg                            <  2:1774488/150‑1 (MQ=255)
                                                                                                                           actcaattGGCCTCTTCGCGGATAAACCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAgg                            >  2:4338413/8‑150 (MQ=38)
                                                                                                                               gTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCt                        <  2:1877777/150‑1 (MQ=255)
                                                                                                                                attGGCCTCTTCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCtt                       >  1:2984111/3‑150 (MQ=39)
                                                                                                                                  tGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa                     >  1:1249371/1‑150 (MQ=255)
                                                                                                                                  tGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGa                     >  2:644341/1‑150 (MQ=255)
                                                                                                                                   ggCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGaa                    <  1:260570/150‑1 (MQ=255)
                                                                                                                                    gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAAcgc                                                                                <  1:3878916/89‑1 (MQ=255)
                                                                                                                                    gCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAg                   >  1:3039282/1‑150 (MQ=255)
                                                                                                                                      cTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCt                 >  1:3430145/1‑150 (MQ=255)
                                                                                                                                       tcttCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTg                >  2:2010752/3‑150 (MQ=255)
                                                                                                                                       tATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTg                >  1:3901322/1‑150 (MQ=255)
                                                                                                                                         ttCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGcc              >  1:4369195/1‑150 (MQ=255)
                                                                                                                                          tCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc             >  1:2472347/1‑150 (MQ=255)
                                                                                                                                          tCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc             <  2:1686863/150‑1 (MQ=255)
                                                                                                                                            gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGa           >  2:3855468/1‑150 (MQ=255)
                                                                                                                                            gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGa           >  2:1975134/1‑150 (MQ=255)
                                                                                                                                            gcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGa           >  1:4933584/1‑150 (MQ=255)
                                                                                                                                             cgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGAc          >  1:2546080/1‑150 (MQ=255)
                                                                                                                                               ggTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTg        <  1:4690084/150‑1 (MQ=255)
                                                                                                                                                  ttATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcac     <  1:1423840/150‑1 (MQ=255)
                                                                                                                                                    aTCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcacca   <  2:4415899/150‑1 (MQ=255)
                                                                                                                                                     tCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcaccac  >  1:3764763/1‑150 (MQ=255)
                                                                                                                                                     tCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGcaccac  >  1:3165418/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCCGACTGCACCAC  >  NC_002947/1580045‑1580343

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: