Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,144 G→A 5.2% intergenic (+22/+79) PP_1661 → / ← PP_1662 dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,1440GA 5.2% 358.6 / 5.9 116intergenic (+22/+79)PP_1661/PP_1662dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (59/51);  new base A (2/4);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 119 total aligned reads displayed.

GCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAT  >  NC_002947/1859002‑1859293
                                                                                                                                              |                                                                                                                                                     
gcgTGGCCAGCAATGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTcc                                                                                                                                                >  1:301365/1‑150 (MQ=255)
      ccAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAAt                                                                                                                                          >  1:1382777/1‑150 (MQ=255)
       cAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATg                                                                                                                                         >  2:1523649/1‑150 (MQ=255)
                     gCAGCTGTGCATCCATGACGAACTGGACGGCCGTTACGGGCGTCGCGTGATCGATAACCTGGTCAGCCTGTCTGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCCGCAATGCCGCAGCATTTTgg                                                                                                                           >  2:1834944/1‑150 (MQ=255)
                     gCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCTGCGGCAttgccg                                                                                                                           <  1:568322/150‑5 (MQ=255)
                     gCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgg                                                                                                                           >  1:2246198/1‑150 (MQ=255)
                     gCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCGGAGCAtcgtgg                                                                                                                           <  2:2000566/150‑6 (MQ=38)
                      cAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCGGAGCAtcgtggg                                                                                                                          <  2:1773081/150‑7 (MQ=38)
                       aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgggg                                                                                                                         >  1:332921/1‑150 (MQ=255)
                         cTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGtt                                                                                                                       >  2:1282639/1‑150 (MQ=255)
                             gCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATAGCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGcca                                                                                                                   <  2:849330/150‑3 (MQ=37)
                              cATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGTGGGGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                  <  1:1021105/150‑1 (MQ=255)
                              cATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                  <  2:1936750/150‑1 (MQ=255)
                               aTCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCTGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTggg                                                                                                                          >  2:1752812/1‑141 (MQ=255)
                               aTCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCTGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTggg                                                                                                                          <  1:1752812/141‑1 (MQ=255)
                               aTCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCAc                                                                                                                 >  2:244065/1‑150 (MQ=255)
                                          ccTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGc                                                                                                      <  1:1518038/150‑1 (MQ=255)
                                          ccTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGc                                                                                                      <  2:2369612/150‑1 (MQ=255)
                                           cTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCg                                                                                                     >  2:2033709/1‑150 (MQ=255)
                                           cTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTAGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCg                                                                                                     >  1:441402/1‑150 (MQ=255)
                                            tGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGa                                                                                                    <  2:2246198/150‑1 (MQ=255)
                                            tGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGa                                                                                                    <  1:20445/150‑1 (MQ=255)
                                                cGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCa                                                                                                >  2:2283938/1‑150 (MQ=255)
                                                   gCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                             <  1:1523649/150‑1 (MQ=255)
                                                     gTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGt                                                                                           <  1:1282639/150‑1 (MQ=255)
                                                       tACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTc                                                                                         <  1:2283938/150‑1 (MQ=255)
                                                          ggggCTGGGCTGAGCGATAGCCTGGGCAGCCTGTGGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAGCCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTGACTGAGCCCAGGCGAGTGAGGGGTTcggc                                                                                      <  2:1740371/150‑3 (MQ=255)
                                                          ggggCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgcc                                                                                      >  2:301557/1‑150 (MQ=255)
                                                               tGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgggg                                                                                            >  1:900727/1‑139 (MQ=255)
                                                               tGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgggg                                                                                            <  2:900727/139‑1 (MQ=255)
                                                               tGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAAt                                                                                 <  2:332921/150‑1 (MQ=255)
                                                               tGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAAt                                                                                 >  2:577698/1‑150 (MQ=255)
                                                                  ccTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATggg                                                                              >  2:411496/1‑150 (MQ=255)
                                                                   cTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGc                                                                             >  1:1063239/1‑150 (MQ=255)
                                                                    tGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCa                                                                            >  1:301870/1‑150 (MQ=255)
                                                                    tGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCa                                                                            >  2:1564905/1‑150 (MQ=255)
                                                                     gAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAg                                                                           <  1:1937603/150‑1 (MQ=255)
                                                                      aGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                          >  1:1051342/1‑150 (MQ=255)
                                                                      aGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCACCGATGCTCCGGCAATGCCGCAGCATTGTTTGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                          <  1:296487/150‑1 (MQ=255)
                                                                       gCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGGTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAggg                                                                         <  2:1127861/150‑1 (MQ=255)
                                                                         gATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTc                                                                       >  2:951188/1‑150 (MQ=255)
                                                                           tAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAt                                                                     >  2:2082257/1‑150 (MQ=255)
                                                                             cccggggCTGCCTGTCGGTAGAGGTGCGTGTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGGATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtg                                                                   <  2:333300/146‑1 (MQ=255)
                                                                              ccTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgt                                                                  >  1:2227439/1‑150 (MQ=255)
                                                                              ccTGGGCAGCCTGTCGGTAGAGATGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgt                                                                  >  2:2363529/1‑150 (MQ=255)
                                                                               cTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCt                                                                                                                    >  2:1181285/1‑99 (MQ=255)
                                                                               cTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCt                                                                                                                    <  1:1181285/99‑1 (MQ=255)
                                                                               cTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtg                                                                 <  1:1206370/150‑1 (MQ=255)
                                                                                tGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgt                                                                <  2:2167308/150‑1 (MQ=255)
                                                                                 gggCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                               <  1:951188/150‑1 (MQ=255)
                                                                                     aGCCTGTCGTTATAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGgcgc                                                           >  2:602713/1‑150 (MQ=255)
                                                                                        cTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCt                                                        >  2:2084405/1‑150 (MQ=255)
                                                                                             ggTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                   <  2:301870/150‑1 (MQ=255)
                                                                                              gTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGtttt                                                  <  2:1051342/150‑1 (MQ=255)
                                                                                               tAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGTTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTc                                                 >  2:185185/1‑150 (MQ=255)
                                                                                                agagGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                                >  1:1146993/1‑150 (MQ=255)
                                                                                                   ggTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAgg                                             <  2:1382777/150‑1 (MQ=255)
                                                                                                    gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGt                                            >  2:21058/1‑150 (MQ=255)
                                                                                                    gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGt                                            >  1:1075513/1‑150 (MQ=255)
                                                                                                         tttCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCa                                       >  2:1759704/1‑150 (MQ=255)
                                                                                                          ttCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCaa                                      <  2:1585780/150‑1 (MQ=255)
                                                                                                          ttCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCaa                                      <  1:1759704/150‑1 (MQ=255)
                                                                                                              gggCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTc                                  >  1:1780265/1‑150 (MQ=255)
                                                                                                                gCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCaa                                                                                  >  1:2009737/1‑100 (MQ=255)
                                                                                                                gCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCaa                                                                                  <  2:2009737/100‑1 (MQ=255)
                                                                                                                gCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCaa                                >  2:2012175/1‑150 (MQ=255)
                                                                                                                   gtgcAGGCCCTGTTATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGCGTTGCCCACTGCGCCCAGGCGAGGCAGGGTTTCGCCGCAATGGGCAGGGGCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAAcgc                             >  2:211592/1‑150 (MQ=255)
                                                                                                                    tgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGtt                                                    >  2:1732190/1‑126 (MQ=255)
                                                                                                                    tgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGtt                                                    <  1:1732190/126‑1 (MQ=255)
                                                                                                                     gtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                           <  2:1063239/150‑1 (MQ=255)
                                                                                                                        ggCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTAt                        >  1:489154/1‑150 (MQ=255)
                                                                                                                         gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTAtt                       >  1:368346/1‑150 (MQ=255)
                                                                                                                          cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTc                                                                                                 <  1:1677655/75‑1 (MQ=255)
                                                                                                                          cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                      <  2:2337389/150‑1 (MQ=255)
                                                                                                                          cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                      >  2:555091/1‑150 (MQ=255)
                                                                                                                           ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                     >  1:1182810/1‑150 (MQ=255)
                                                                                                                           ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                     >  2:77930/1‑150 (MQ=255)
                                                                                                                           ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGGCATTGTGTGACGCGCTTGTTGACCAGGTCTTCAAGCGCAACGCGGTATTgg                     >  2:1109378/1‑150 (MQ=255)
                                                                                                                            cTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgt                    <  1:2363529/150‑1 (MQ=255)
                                                                                                                            cTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATAGTGTGACGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgt                    <  1:1161395/150‑1 (MQ=255)
                                                                                                                             tGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGt                                                                                                                        <  1:1009585/49‑1 (MQ=255)
                                                                                                                             tGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGt                                                                                                                        >  2:1009585/1‑49 (MQ=255)
                                                                                                                             tGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtg                   >  2:1998115/1‑150 (MQ=255)
                                                                                                                             tGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGAGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtg                   >  2:1647697/1‑150 (MQ=255)
                                                                                                                              gTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtgt                  <  2:1933551/150‑1 (MQ=255)
                                                                                                                              gTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCGGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtgt                  >  1:1933551/1‑150 (MQ=255)
                                                                                                                                   ggCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc             <  1:555091/150‑1 (MQ=255)
                                                                                                                                   ggCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc             >  2:586104/1‑150 (MQ=255)
                                                                                                                                    gCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCg            >  2:8567/1‑150 (MQ=255)
                                                                                                                                    gCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCg            >  2:591188/1‑150 (MQ=38)
                                                                                                                                      aaCCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                           >  2:1584014/1‑133 (MQ=255)
                                                                                                                                      aaCCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                           <  1:1584014/133‑1 (MQ=255)
                                                                                                                                         cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc       <  1:244065/150‑1 (MQ=255)
                                                                                                                                          ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                               <  1:1908791/93‑1 (MQ=255)
                                                                                                                                          ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                               >  2:1908791/1‑93 (MQ=255)
                                                                                                                                            aCGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAtt                                                                    >  1:1065495/1‑86 (MQ=255)
                                                                                                                                            aCGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGtt    <  2:1780265/150‑1 (MQ=255)
                                                                                                                                             cGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTa   <  1:774636/150‑1 (MQ=255)
                                                                                                                                              gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGTTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAt  <  1:185185/150‑1 (MQ=255)
                                                                                                                                              gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAt  <  1:577698/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAT  >  NC_002947/1859002‑1859293

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