Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 746,976 C→A 6.3% intergenic (‑100/‑694) PP_0638 ← / → PP_5448 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947746,9760CA 6.3% 201.5 / 2.9 63intergenic (‑100/‑694)PP_0638/PP_5448transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base A (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

GGCTTGGGGTAGAGGTAGACTTTTTCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCAATGACGTCACC  >  NC_002947/746830‑747125
                                                                                                                                                  |                                                                                                                                                     
ggCTTGGTGTAGAGGTAGACTTTTTCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGACCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCa                                                                                                                                                    <  2:1182819/150‑1 (MQ=255)
       ggtagaggtagaCTTTTTCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGccacca                                                                                                                                             <  2:1764323/150‑1 (MQ=255)
         tagaggtagaCTTTTTCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATa                                                                                                                                           <  1:1501224/150‑1 (MQ=255)
                      tttCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCCGTTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGc                                                                                                                              >  1:898410/1‑150 (MQ=255)
                      tttCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACACCAGCCACCATAGCTCAATGGCTGc                                                                                                                              >  1:2335183/1‑150 (MQ=255)
                          gACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTc                                                                                                                          >  1:214633/1‑150 (MQ=255)
                             tttGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCa                                                                                                                                                    >  1:1864661/1‑121 (MQ=255)
                             tttGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCa                                                                                                                                                    <  2:1864661/121‑1 (MQ=255)
                             tttGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGc                                                                                                                       >  1:1302032/1‑150 (MQ=255)
                             gtggCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGc                                                                                                                       >  2:2212213/4‑150 (MQ=255)
                               tGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCt                                                                                                                     >  1:1100265/1‑150 (MQ=255)
                                 cTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTcc                                                                                                                   <  2:1100265/150‑1 (MQ=255)
                                  tGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCa                                                                                                                  >  1:717056/1‑150 (MQ=255)
                                    tCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAAc                                                                                                                >  2:1671923/1‑150 (MQ=255)
                                       gggCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATATGGGGAGTATGGATCACTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCt                                                                                                             >  2:573435/1‑150 (MQ=255)
                                          cgcATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATTTCTCCGTTCAGCCTCCAACCCTCAg                                                                                                          <  1:573435/150‑1 (MQ=255)
                                          cgcATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGACGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAg                                                                                                          <  2:1693448/150‑1 (MQ=255)
                                           gcATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGc                                                                                                         <  1:1671923/150‑1 (MQ=255)
                                              tcatGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCggc                                                                                                      >  2:2276267/1‑150 (MQ=255)
                                                  ggCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACACCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCggcag                                                                                                    <  2:772204/148‑1 (MQ=255)
                                                  ggCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACACCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCggcag                                                                                                    >  1:772204/1‑148 (MQ=255)
                                                          ggctCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGCGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGt                                                                                          <  1:2276267/150‑1 (MQ=255)
                                                           gctCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTg                                                                                         >  1:1675851/1‑150 (MQ=255)
                                                             tCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGat                                                                                       >  1:2191375/1‑150 (MQ=255)
                                                               cAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGatat                                                                                     <  2:717056/150‑1 (MQ=255)
                                                               cAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGatat                                                                                     >  1:1596374/1‑150 (MQ=255)
                                                                agaaagaaATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATc                                                                                    >  2:198762/1‑150 (MQ=255)
                                                                agaaagaaATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATc                                                                                    >  1:681440/1‑150 (MQ=255)
                                                                 gaaagaaaTCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCg                                                                                   <  2:1675851/150‑1 (MQ=255)
                                                                          cGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGAtt                                                                          <  1:2212213/150‑1 (MQ=255)
                                                                            ggAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAATCCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCaa                                                                                                                 <  1:2301087/109‑1 (MQ=255)
                                                                            ggAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCaa                                                                                                                 >  2:2301087/1‑109 (MQ=255)
                                                                             gAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCa                                                                                                           <  1:475199/114‑1 (MQ=255)
                                                                             gAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCa                                                                                                           >  2:475199/1‑114 (MQ=255)
                                                                               gCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTg                                                                     >  2:1494365/1‑150 (MQ=255)
                                                                                  caGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTgaga                                                                  <  2:1555515/150‑1 (MQ=255)
                                                                                   aGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGAt                                                                 >  1:522527/1‑150 (MQ=255)
                                                                                      aTCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCa                                                              >  2:1683037/1‑150 (MQ=255)
                                                                                         gggAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAg                                                           <  2:681440/150‑1 (MQ=255)
                                                                                                aGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAAttt                                                    >  1:266951/1‑150 (MQ=255)
                                                                                                   tGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCgg                                                 <  2:1596374/150‑1 (MQ=255)
                                                                                                   tGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCgg                                                 >  2:425853/1‑150 (MQ=255)
                                                                                                                tgtgGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGc                                                                                              >  2:717372/1‑92 (MQ=255)
                                                                                                                tgtgGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGc                                                                                              <  1:717372/92‑1 (MQ=255)
                                                                                                                   ggggATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCac                                 >  2:1939057/1‑150 (MQ=255)
                                                                                                                      gatgatGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGTTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCaccga                              <  2:214633/150‑1 (MQ=255)
                                                                                                                      gatgatGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCaccga                              <  1:1939057/150‑1 (MQ=255)
                                                                                                                                cccTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGt                    <  2:1302032/150‑1 (MQ=255)
                                                                                                                                 ccTTACGCTTAGAAAACACCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTg                   <  2:2335183/150‑1 (MQ=255)
                                                                                                                                  cTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCt                                                                      <  1:453541/98‑1 (MQ=255)
                                                                                                                                  cTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCt                                                                      >  2:453541/1‑98 (MQ=255)
                                                                                                                                  cTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTaa                                             >  1:1813851/1‑123 (MQ=255)
                                                                                                                                  cTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTaa                                             <  2:1813851/123‑1 (MQ=255)
                                                                                                                                     aCGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTgg               <  1:860800/150‑1 (MQ=255)
                                                                                                                                      cGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGc              <  2:266951/150‑1 (MQ=255)
                                                                                                                                       gCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCa             >  1:210093/1‑150 (MQ=255)
                                                                                                                                        cTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCaa            <  1:198762/150‑1 (MQ=255)
                                                                                                                                           aGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTga                                                                    >  2:2018629/1‑91 (MQ=255)
                                                                                                                                           aGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTga                                                                    <  1:2018629/91‑1 (MQ=255)
                                                                                                                                            gAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCAATGAc        <  1:1683037/150‑1 (MQ=255)
                                                                                                                                            aaaacCCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGGGATGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCAATGAc        <  1:25209/149‑1 (MQ=255)
                                                                                                                                              aaaCCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTg                      >  1:1572021/1‑134 (MQ=255)
                                                                                                                                              aaaCCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTg                      <  2:1572021/134‑1 (MQ=255)
                                                                                                                                                 ccccAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCAATGACGTCAc   >  1:379621/1‑150 (MQ=255)
                                                                                                                                                  cccAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCAATGACGTCAcc  >  1:1430748/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCTTGGGGTAGAGGTAGACTTTTTCGACTTTGCTGTCGGGGCGCATCATGGCGGGCTGGCTCCAGAAAGAAATCGGGAGCACAGCATCGGGAATCAGGTGGACGGTTGGAATGTGGGGATGATGGATCCCTTACGCTTAGAAAACCCCAGCCACCATAGCTCAATGGCTGCGTTCAGCCTCCAACCCTCAGCGGCAGGGCAGCTCGTGATATCGCTCACGATTGGCTGAGATTCATAGGAAATTTCGGGTAAGAGGGGCGGCACCGACCGTTCTGGTGCTGGCAATGACGTCACC  >  NC_002947/746830‑747125

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