Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,079 T→G 13.7% intergenic (+40/+112) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,0790TG13.7% 518.9 / 21.6 237intergenic (+40/+112)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (56/138);  new base G (34/4);  total (90/147)
Fisher's exact test for biased strand distribution p-value = 2.01e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

TGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCA  >  NC_002947/667930‑668228
                                                                                                                                                     |                                                                                                                                                     
   ccccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGC                                                                                                                                                    >  2:5680352‑M1/1‑150 (MQ=255)
     ccGGGCACCTACCGCGTGCCGCCGCCGCCGGTCGACTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGA                                                                                                                                                  <  1:5675961‑M1/150‑1 (MQ=255)
     ccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGA                                                                                                                                                  <  2:4075375‑M1/150‑1 (MQ=255)
        ggCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                               <  1:2086792/150‑1 (MQ=255)
        ggCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                               <  1:6866641/150‑1 (MQ=255)
            ccTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                           <  2:5360277/150‑1 (MQ=255)
             cTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCg                                                                                                                                          >  1:3363275/1‑150 (MQ=255)
              tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTGGCGATGGGGCCGc                                                                                                                                         >  1:7251433/1‑150 (MQ=255)
              tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAtagggctgc                                                                                                                                         <  2:4612694‑M1/150‑10 (MQ=255)
               aCCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGGTAGGGCCGCa                                                                                                                                        >  2:3251487/1‑150 (MQ=255)
                   cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                    <  2:8787190/150‑1 (MQ=255)
                   cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                    <  2:5750114/150‑1 (MQ=255)
                    gTGCCGCCGCCGCAGGTCGAATCGAGGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGACAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCg                                                                                                                                   <  2:7659876/150‑1 (MQ=255)
                      gccgccgccGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                                 <  1:3719032/150‑1 (MQ=255)
                       ccgccgccgCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagcc                                                                                                                                <  2:4696948‑M1/150‑19 (MQ=255)
                       ccgccgccgCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                                <  2:2514227/150‑1 (MQ=255)
                        cgccgccgcAGGTCGACTCGAGGTACCGGCCGAACCTGCAGGCCGTGGGCGCTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                               <  2:5708047/150‑1 (MQ=255)
                         gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGCGCGGTCGGGCCGCAAAGCGGcccc                                                                                                                              >  1:8809914/1‑150 (MQ=255)
                          ccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGcccta                                                                                                                             <  1:386299/150‑3 (MQ=255)
                           cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctat                                                                                                                            <  2:2481614/150‑4 (MQ=255)
                               gcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACCAt                                                                                                                        >  2:6372640/1‑150 (MQ=255)
                                 aGGTCGAACCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCATGGGGAGTTGACGGTGAAAGCCCGATTGAGACCCGCTGTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTc                                                                                                                      <  1:7083518/150‑1 (MQ=38)
                                      gactCGATGTACCGGCCGACCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcg                                                                                                                 <  2:7104983/147‑1 (MQ=255)
                                      gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCCAAGCGGCCCCGACAATTCAAgcg                                                                                                                 >  2:8521170/1‑150 (MQ=255)
                                       aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTTGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                <  1:8702533/150‑1 (MQ=255)
                                       aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGGGCGGTTGGGCCCCCAACCGGCCCCCCCCAATCAAgccg                                                                                                                >  1:7469873/1‑148 (MQ=255)
                                       aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCGGGGCGGCGCGGTGGGGCGGCCGAGCGGCGCCGACCAttacaccgg                                                                                                                >  1:2503843/1‑143 (MQ=255)
                                       aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                >  2:8205307/1‑150 (MQ=255)
                                        aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCGCAAAGCGGCCCCGACCATTCAAgcggc                                                                                                               >  2:848578/1‑150 (MQ=255)
                                           gATGTACGGGCCGACCCTGCCGGCGGTGGGCGATGGGCATGGGGTTATGACGGTGAAAGCCCGCTAGAGACCCGATCTACGCGTAACGTGTGGGAGGGGCCTTGCTTCGCTTTTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                            <  1:6932078/150‑1 (MQ=255)
                                            ttGTACCGGCCGAACCTGCAGGCGGTGGGCGTTGGGCATGGGTATATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                           <  1:1351867/149‑1 (MQ=38)
                                              gTACCGGCCGAACCTGCAGGCGGGGGGCGAGGGGCCTGGGGATATGCGGGTGACCGCCCGCTAGCGACCCGCTCTCCGCGGAAAGTGGGGGCGGGGCCTTGCGTCGCGATCGGGCCGAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                         <  1:6611502/150‑1 (MQ=255)
                                              gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                         <  2:2985344/150‑1 (MQ=255)
                                              gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                         <  2:7318515/150‑1 (MQ=255)
                                              gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                         <  2:654187/150‑1 (MQ=255)
                                              gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                         <  2:7892158/150‑1 (MQ=255)
                                               gaCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTg                                                                                                        <  2:5262133/149‑1 (MQ=255)
                                                aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGCGATAGGGCCGCCCAGCGGCCCCGACCAATCAAGCGGCGAAGCGGa                                                                                                       >  2:5875168/1‑150 (MQ=255)
                                                aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCGGa                                                                                                       >  2:452238/1‑150 (MQ=255)
                                                aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCCCAAAGCCGCCCCGACAATACAATCGgcggagcgga                                                                                                       >  2:83580/1‑150 (MQ=255)
                                                 ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCCAAGCGGCCCCCACAATTCAAGCGGCGAAGCGGaa                                                                                                      >  2:2207322/1‑150 (MQ=255)
                                                   ggCCGAACCTGCTGGCCGTGGGCGTTGGGCATGGGGAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGTGCGGCCTTGCGTCGCTTTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAg                                                                                                    <  1:5183488/150‑1 (MQ=255)
                                                       gaccgtGCTGGCGGTGGGCGTTGGGCAAGGGGAGATGAGGGTTAAAGCCCGCTCGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                <  1:5426650/145‑1 (MQ=37)
                                                       gacccTGCAGGCCTTGGGCTTTGGGAAAGGGTTGATGACGGTTAAAGCCCTCTAGAGACCCGCTCTCCGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCTATTGGGCCCAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                <  1:7714657/147‑1 (MQ=37)
                                                        aaCCTGCAGGCGGTGGGCGTTGGGAATGGGTATTTGACGTTTAAAGCCCGCTTGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCGTTGCGTCGCGATTGGGCCGAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  1:6182187/150‑1 (MQ=255)
                                                        aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  2:1271207/150‑1 (MQ=255)
                                                        aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  2:1913650/150‑1 (MQ=255)
                                                          ccTGCTGGCGGTGGGCGAGGGGCATGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGGGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                             <  1:4480142/150‑1 (MQ=37)
                                                          ccTGCTGGCCGTGGGCGTGGGGCATGGGGTGCTGACGGTGAAAGCCCGCTAGAGCCCCGCTCTAGGGGGAATGTGTGGGCGCGGCCTTGCGTCGCGATTGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                             <  1:965880/150‑1 (MQ=255)
                                                               aGGCCGTGGGCGATGGGCAAGGGGAGATGACTGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCGCGCGATAGGGCCGCAAAGCCGCACCGACAAATCAAGCGCCCAACCAGACAGACTGGGgcagca                                                                                        >  2:7452607/1‑149 (MQ=255)
                                                               aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGCGGCGCGCTTGGGCCGCCAAGCGGCCCCGCCCCGTGACGCGGCGAAGCGGAAAGCCTGGGgcggct                                                                                        >  1:3617728/1‑150 (MQ=255)
                                                                  ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGGGCGGTAGGGCCGCCCCGCGGCCCCGCCCAAAAAAGCCGCGCAGCGGCGAGCCGGGggcggcgggg                                                                                     >  1:2715387/1‑146 (MQ=255)
                                                                   cGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGc                                                                                    <  2:4968250/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:8618335/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGAAGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGAGTGGGAGCGGCCATGCCGCGCCATAGCGCCCCCCAGCCGCCCCCACAAATCAACCGGCGCACCGGCAACCCTGGGGCGGCtgtgaa                                                                                   >  2:3805193/1‑148 (MQ=255)
                                                                      gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggc                                                                                 <  2:4264194‑M1/150‑66 (MQ=255)
                                                                       ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcc                                                                                <  1:4078085‑M1/150‑67 (MQ=255)
                                                                       ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                <  2:7058323/150‑1 (MQ=255)
                                                                         ggtggggCACGGGGTGATGACGTTGAAAGCCCGCTTGAGACCCGCTCTACGCGGAACGTGTGGGCGCGGCCTTGCGTCGCGATCGGGCCGAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                              <  1:2990447/146‑1 (MQ=255)
                                                                          gATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCACAGCAGcccca                                                                             <  2:226726/150‑3 (MQ=255)
                                                                            tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatc                                                                           <  2:8751868‑M1/150‑72 (MQ=255)
                                                                            tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggcagcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatc                                                                           <  1:2758201‑M1/150‑72 (MQ=255)
                                                                             gggCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGGGTTAGGGCCGCAAAGCGGCCCCGTCAAAGGAAGCGGCGAAGGGGAAAGCCGGGGGCTGCTGTGGAGcccaggcg                                                                          >  1:6464384/1‑145 (MQ=255)
                                                                                cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCAGGCGTGGCGATAGGGCCGCCAAGCGGCCCCGACCATACAAGCGGCCAAGCGGTACGCCAGGGGCGGCGGTGCCGCCCAGGCGCGa                                                                       >  1:1027646/1‑150 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCTCGCGATAGGGCCGCCAAGCGGCCCCCACAAATCAAGCGGCGAACCGGAAAGCCTGGGGCTGCAGTGCAGCCCGGtcgcgac                                                                      <  2:4440916‑M2/150‑8 (MQ=255)
                                                                                    gggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAACCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgca                                                                   <  2:8424198‑M2/150‑11 (MQ=255)
                                                                                     gttgaTGCCGTTCAAAGCCCGATATATCCCCGATATACGCGGAATGTGTGGTTGGGTCGTTGCTTCGCGATAGGGCCCCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCaa                                                                  <  1:1608485/147‑1 (MQ=255)
                                                                                     ggtgttgacggTTAAAGCCCGCTAGAGACCCGCTATACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                                  >  1:8363029‑M2/6‑139 (MQ=255)
                                                                                      gagaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaag                                                                 <  1:8346003‑M1/150‑82 (MQ=255)
                                                                                      gagaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaag                                                                 >  2:2995081‑M2/1‑138 (MQ=255)
                                                                                        gttgaGGGTCAACGCCCGCTATACACCGGCTCTACGGGTAACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggc                                                               >  1:1751967‑M2/3‑136 (MQ=255)
                                                                                           gACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATAAAGCGGCGAAGCTGAAAGCCAGGGGCTGCGGTGCAGCCCAGtcgcgacgccaggccgc                                                            <  2:2588190‑M2/150‑18 (MQ=255)
                                                                                           gACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacccaaggccgc                                                            <  2:2588624‑M2/150‑18 (MQ=255)
                                                                                              ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGGCCCGACAAATATAGCGGCGCAGCTGTAAGCCTGGGGCTGCTGTGCAGCCCGGtggcggcgcaaggccgctcc                                                         <  2:2567195‑M2/150‑21 (MQ=255)
                                                                                                    aGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcc                                                                                <  2:5246947‑M1/121‑67 (MQ=255)
                                                                                                           ttgagcCCCGCTCTCCGGGGACTGTGTGGGCGCGGCCTGGCGCCGCTATTGGGCCCCAACGCGGCCCCGACAATTCACGGGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGATCCCACAAAGACcgcg                                            <  2:1704375/148‑1 (MQ=255)
                                                                                                             gagaCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCgcgcc                                          >  2:8062925/1‑149 (MQ=255)
                                                                                                               caCCCGCTCTCCGCGGACTGTGGGGGCGCGGCCTTGCGTCGCGATGGGGCCGAAACGCGGCCCCGACACTTCACGCGGCGACGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGTAGa                                        <  2:8242546/149‑1 (MQ=25)
                                                                                                                aCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATCAAGCGGCGGAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCAAAGACCgcgccgct                                       >  2:3398330/1‑146 (MQ=37)
                                                                                                                        tACGCGGATTGTGGGGGCGCGGCGTGGGGTCGCTATCGGGCCCCAAGGCGCCCCCGACAATTCACGCGGCGACGCTGAACGCGTGGGGCTGGGGTGCCGCCGTGCCGCGCGGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGAt                               <  2:2096383/150‑1 (MQ=255)
                                                                                                                          cgcgGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCCGCCCCGACAAATCAAGCGGCCCAGCTGAAAGCCGGGGGCTGCGGGGCAGCCCAGTCGCGACGCAAGGCCGCGCCCACCAAGACCGCGCAGATCAGTGGATCa                             >  2:3982047/1‑150 (MQ=37)
                                                                                                                                 ttggtgggaGCGGCCTTGCGTCGCGActgggctgcccagcagccccaggctttcagcttcgccgctttaattgtcggggccgctgtgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggct                      <  2:2039686‑M1/147‑125 (MQ=255)
                                                                                                                                  gtgtGGGAGCGGGCTTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCGACAATTAAAGCGGCGAAGCTGAAAGCCAGGGGCTGCTGTGGCGCCCTGTCGCGACGCAAGGCCGCTCCCCCAAAGACCGCGCACATCAGTGGATCACCCAGCtc                     >  2:4888678/1‑150 (MQ=37)
                                                                                                                                   tgtgGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggcttt                    <  2:7855435‑M1/150‑127 (MQ=255)
                                                                                                                                     tgGGCGGGGCCTTGGGTCGCTTTCGGGCCGCAACGGGGCCCCGACAATTCACGCGGCCAAGCTGAAAGCCTGGGGCTGCTTTGCAGCCCTTTGGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtcatc                  <  2:1359164/150‑1 (MQ=25)
                                                                                                                                       ggcgcggCCTGGCGTCGCGAcagggctgcaaagccgcccccggcttacagcttcgccgctggaattgtcggggccgctttgcggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcacc                <  2:5209047‑M1/147‑131 (MQ=255)
                                                                                                                                        gcgCGGCCTTGCGCCGCTATCGGGCCGAAACGCGGCCCCGACAATTAACGCGGCGACGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAcc               <  2:4346061/148‑1 (MQ=37)
                                                                                                                                         cgcggcgTTGCGTCGCGTACGGCCCGAACCCCGGCCCCGACAATCCCACCGCGGCCGCTGACACCCGGGGCCGGCGGTCCCGCCCGGTCCCCACGCACGCCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc              <  2:8284486/149‑1 (MQ=255)
                                                                                                                                         aGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc              >  2:8412343/1‑150 (MQ=255)
                                                                                                                                           cGGCCTTGCGTCGCGActgggctgcacagcagcccccggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtca            <  2:4509690‑M1/150‑135 (MQ=255)
                                                                                                                                             gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatc          <  1:6001027‑M2/150‑68 (MQ=255)
                                                                                                                                              ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgcacccacactttccgcgtagagcgggtgtatagcgggctttcaacgtcatct         >  1:2691811‑M1/1‑13 (MQ=255)
                                                                                                                                              ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcccgatcgcgacgcaaggcccctcccacactttccgcggagagcgggtctctagcgggctttcaccggcctca         >  1:946702‑M1/1‑13 (MQ=255)
                                                                                                                                              ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctgtgcggccctatcgcgacgcaaggccgctcccacactgtccgcgtagagcgggtctctagcgggctttcaccgtcatct         >  2:3507025‑M1/1‑13 (MQ=255)
                                                                                                                                               cTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctc        >  1:1399461‑M1/1‑12 (MQ=255)
                                                                                                                                               cTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttacgcgtagagcgggtctctagcgggctttcacagtcatctc        >  1:3594363‑M1/1‑12 (MQ=255)
                                                                                                                                                    gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccctt   <  1:2219150‑M2/150‑75 (MQ=255)
                                                                                                                                                    gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc   >  1:5825597/1‑150 (MQ=255)
                                                                                                                                                    gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc   >  1:5825311/1‑150 (MQ=255)
                                                                                                                                                     tCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttg  >  2:8052571‑M1/1‑6 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCA  >  NC_002947/667930‑668228

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