Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,897,245 G→T 6.4% intergenic (‑56/+88) PP_mr24 ← / ← rdgC mini‑ykkC regulatory region/recombination‑associated protein RdgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,897,2450GT 6.4% 298.1 / 6.3 94intergenic (‑56/+88)PP_mr24/rdgCmini‑ykkC regulatory region/recombination‑associated protein RdgC
Reads supporting (aligned to +/- strand):  ref base G (46/42);  new base T (4/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

AACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAATTC  >  NC_002947/1897096‑1897392
                                                                                                                                                     |                                                                                                                                                   
aaCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAgg                                                                                                                                                     <  2:639145/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                    <  1:1780461/150‑1 (MQ=255)
  cAGGATGATCCAGGACATGAGTGAGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAgggg                                                                                                                                                   <  1:951366/150‑1 (MQ=255)
   aGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                    >  1:89379/1‑148 (MQ=255)
   aGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                    <  2:89379/148‑1 (MQ=255)
     gatgatCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTa                                                                                                                                                >  2:723778/1‑150 (MQ=255)
     gatgatCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTa                                                                                                                                                <  2:1345425/150‑1 (MQ=255)
        gatCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTAAAAAAGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCa                                                                                                                                             >  1:1206078/1‑150 (MQ=255)
         atCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTcaccac                                                                                                                                            >  1:257044/1‑150 (MQ=255)
          tCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGTAAGGGGTTACCACt                                                                                                                                           >  1:941201/1‑150 (MQ=255)
              ggACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCa                                                                                                                                       <  2:1634429/150‑1 (MQ=255)
                  aTGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGgcg                                                                                                                                   >  2:1261000/1‑150 (MQ=255)
                   tGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGTGGTAACCCCTTGCAGgcgc                                                                                                                                  <  1:1758711/150‑1 (MQ=255)
                     aGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCgg                                                                                                                                <  2:1541214/150‑1 (MQ=255)
                     aGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGACTGCAAGGGGTTACCACTTGCAGGCGCgg                                                                                                                                <  1:1180939/150‑1 (MQ=255)
                         tgtCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGttt                                                                                                                            >  2:1124415/1‑150 (MQ=255)
                           tctcCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTAc                                                                                                                          >  2:1788952/1‑150 (MQ=255)
                            ctcCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTAcc                                                                                                                         >  1:85560/1‑150 (MQ=255)
                                 tCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCa                                                                                                                                       >  2:2369784/1‑131 (MQ=255)
                                 tCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCa                                                                                                                                       <  1:2369784/131‑1 (MQ=255)
                                  cGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGaa                                                                                                                   <  1:332147/150‑1 (MQ=255)
                                    tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAaga                                                                                                                 >  1:2006164/1‑150 (MQ=255)
                                       gAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACt                                                                                                                                           >  2:1810100/1‑121 (MQ=255)
                                       gAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACt                                                                                                                                           <  1:1810100/121‑1 (MQ=255)
                                        aaTGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGcc                                                                                                             >  1:1644199/1‑150 (MQ=255)
                                               gggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTg                                                                                                      <  2:870832/150‑1 (MQ=255)
                                               gggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTg                                                                                                      >  1:1284611/1‑150 (MQ=255)
                                                ggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTgg                                                                                                     <  2:1644199/150‑1 (MQ=255)
                                                ggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTgg                                                                                                     >  1:2296498/1‑150 (MQ=255)
                                                ggCCGTCCCCGGAGTGGGCACGCGGTGCGTCGTGAGGTCGTCCGCACCGTGGGCACCATGGTGCCTGTGTGCGGCCTTGCGGGCAAGACCGCGCCGGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTgg                                                                                                     <  2:340908/150‑1 (MQ=255)
                                                  ccGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                            <  1:85679/141‑1 (MQ=255)
                                                  ccGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                            >  2:85679/1‑141 (MQ=255)
                                                        ccGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAAc                                                                                             <  1:1884055/150‑1 (MQ=255)
                                                        ccGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAAc                                                                                             >  1:2264203/1‑150 (MQ=255)
                                                          ggAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGa                                                                                           <  2:85560/150‑1 (MQ=255)
                                                          ggAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGa                                                                                           <  1:2033608/150‑1 (MQ=255)
                                                          ggAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTAGCAGAAAACGa                                                                                           <  1:1835621/150‑1 (MQ=255)
                                                               ggTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTCACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAg                                                                                      <  2:257044/150‑1 (MQ=255)
                                                                atCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTACCTGTCTGCGGCCTTGCGGGCAGGACCGCGCCTGCAAGGGGTTACCACTTGCCGGCGCGGGTGTACCCGCGACGACGACGGAGGTGGCGCAAAACGATGCacg                                                                                     >  1:881617/2‑148 (MQ=255)
                                                                 tCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGc                                                                                    <  1:1261000/150‑1 (MQ=255)
                                                                   atgcgGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGAGTCAGGCCa                                                                                  >  2:861468/3‑150 (MQ=255)
                                                                            acggtGAGGACGACCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTa                                                                                                                                                >  2:8136/4‑79 (MQ=255)
                                                                            acggtGAGGACGACCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTa                                                                                                                                                <  1:8136/76‑1 (MQ=255)
                                                                               gtGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCAttt                                                                      >  1:2166218/1‑150 (MQ=255)
                                                                                tGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGCTTTTCCCGCTAAGATGCCGGTGCTCGCAGAAAACGATTCAGGCCAGCACCTGGAtttt                                                                     >  2:1144811/1‑149 (MQ=255)
                                                                                   ggTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATa                                                                  <  1:1018045/150‑1 (MQ=255)
                                                                                   ggTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATa                                                                  >  1:1891188/1‑150 (MQ=255)
                                                                                    gtcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAAc                                                                                             >  1:1045356/1‑122 (MQ=255)
                                                                                    gtcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAAc                                                                                             <  2:1045356/122‑1 (MQ=255)
                                                                                    gtcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAc                                                                 <  2:1284611/150‑1 (MQ=255)
                                                                                    gtcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAc                                                                 <  2:1163420/150‑1 (MQ=255)
                                                                                     tcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAcc                                                                >  1:1491110/1‑150 (MQ=255)
                                                                                       gtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCaa                                                              >  1:1526400/1‑150 (MQ=255)
                                                                                       gtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGGGCTGGCAGAAAACGATTCAGGCCAGCGCCTGCATGTCATACCaa                                                              >  1:153401/1‑150 (MQ=255)
                                                                                        tcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAAt                                                             >  1:999038/1‑150 (MQ=255)
                                                                                                gTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAACAGGAGTCAGGCCAGCACCTGCATTTCATACCAATCAGAcccc                                                     >  1:1840936/1‑150 (MQ=255)
                                                                                                acaggCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGAcccc                                                     >  2:2345446/4‑150 (MQ=255)
                                                                                                 tGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTg                                                                                                                                         <  1:1280046/65‑1 (MQ=255)
                                                                                                 tGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTg                                                                                                                                         >  2:1280046/1‑65 (MQ=255)
                                                                                                  aggCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                           >  1:1480219/2‑142 (MQ=255)
                                                                                                  aggCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                           <  2:1480219/141‑1 (MQ=255)
                                                                                                    gCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAc                                                                 >  2:1177053/1‑134 (MQ=255)
                                                                                                    gCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAc                                                                 <  1:1177053/134‑1 (MQ=255)
                                                                                                    gCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGc                                                 <  1:2192750/150‑1 (MQ=255)
                                                                                                        cATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGaa                                             >  2:1268530/1‑150 (MQ=255)
                                                                                                        cATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGaa                                             >  1:1375008/1‑150 (MQ=255)
                                                                                                            gTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAAtctc                                         <  1:2345446/150‑1 (MQ=255)
                                                                                                               ccTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGTAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTc                                      <  2:941201/150‑1 (MQ=255)
                                                                                                                 tGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCACGGAACTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTcac                                    <  1:723778/150‑1 (MQ=255)
                                                                                                                   tCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTcacca                                  <  2:1896868/150‑1 (MQ=255)
                                                                                                                         gCCTTGCGGGCAAGACCGCGCTTGCACGGGAATGTTTCTGGGAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAgc                            <  2:1840936/150‑1 (MQ=255)
                                                                                                                                 ggCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                           <  2:1772135/111‑1 (MQ=255)
                                                                                                                                 ggCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                           >  1:1772135/1‑111 (MQ=255)
                                                                                                                                 ggCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                           <  2:566486/111‑1 (MQ=255)
                                                                                                                                 ggCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                           >  1:566486/1‑111 (MQ=255)
                                                                                                                                 ggCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGaa                    <  1:1788952/150‑1 (MQ=255)
                                                                                                                                  gCAAGTCCGCTCCTGCAAGGAGTTACCACTTGCAGGCGCGGGCTTACCAGCGAATAGGCCGGTGCTGGCAGAAAACGATTAAGGCCAGCACCTGCATTTCATACCAATCCGACCCCTCGCGGAACCCCTTCACCACCCCCCTCCCCGAAc                   >  2:862207/1‑150 (MQ=255)
                                                                                                                                  gCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAAc                   >  2:1929225/1‑150 (MQ=255)
                                                                                                                                   cAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCCAGCGGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTGCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACa                  >  1:1636513/1‑150 (MQ=255)
                                                                                                                                    aaGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCtt                                       <  1:2281559/128‑1 (MQ=255)
                                                                                                                                    aaGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCtt                                       >  2:2281559/1‑128 (MQ=255)
                                                                                                                                      gACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGcc               >  2:441272/1‑149 (MQ=255)
                                                                                                                                       aCCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCAc              <  2:2006164/150‑1 (MQ=255)
                                                                                                                                        ccGCGCCTGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCAcc             >  2:1364389/1‑150 (MQ=255)
                                                                                                                                         cgcgCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCg            <  2:905475/150‑1 (MQ=255)
                                                                                                                                             gcggCAAGGGGTTACCAGTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGCTTCAGGCCAGCACCTGCATTTCGTACGACTCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAg        <  2:579820/147‑1 (MQ=255)
                                                                                                                                              cTGCAAGTGGTAACCCCTTGCAGGCGCGGTCTTGCCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAgg       >  1:2026010/1‑150 (MQ=38)
                                                                                                                                               tGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATa                                                                  >  2:1896419/1‑90 (MQ=37)
                                                                                                                                               tGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATa                                                                  <  1:1896419/90‑1 (MQ=37)
                                                                                                                                               tGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGa      >  2:1526416/1‑150 (MQ=255)
                                                                                                                                               tGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGa      >  1:946755/1‑150 (MQ=255)
                                                                                                                                               tGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGa      >  1:5035/1‑150 (MQ=255)
                                                                                                                                                gCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGaa     <  2:2264203/150‑1 (MQ=255)
                                                                                                                                                 cAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAAt    >  1:1055138/1‑150 (MQ=255)
                                                                                                                                                  aaGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCtt                                       <  2:1012063/114‑1 (MQ=255)
                                                                                                                                                  aaGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCtt                                       >  1:1012063/1‑114 (MQ=255)
                                                                                                                                                  aaGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGACtt   >  2:1450023/1‑150 (MQ=255)
                                                                                                                                                   aGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGTAATCGCTTCACCACCCAGCGTCTCGAACAGCACCGGCAGGAATTc  <  2:97079/150‑1 (MQ=255)
                                                                                                                                                   aGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGACTGTTGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAACCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAATTc  <  1:1124415/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAATTC  >  NC_002947/1897096‑1897392

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: