Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,079 T→G 14.8% 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,0790TG14.8% 360.4 / 10.3 168intergenic (+40/+112)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (57/80);  new base G (27/2);  total (85/83)
Fisher's exact test for biased strand distribution p-value = 1.91e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

TGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCA  >  NC_002947/667930‑668228
                                                                                                                                                     |                                                                                                                                                     
tGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGT                                                                                                                                                       <  1:1142128‑M1/150‑1 (MQ=255)
   ccccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGC                                                                                                                                                    <  2:2122333‑M1/150‑1 (MQ=255)
   ccccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGC                                                                                                                                                    <  2:1541533‑M1/150‑1 (MQ=255)
    cccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCG                                                                                                                                                   >  2:3723885‑M1/1‑150 (MQ=255)
           aCCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAggg                                                                                                                                             <  2:4147644/149‑1 (MQ=255)
              tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGATGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGc                                                                                                                                         >  2:1664846/1‑150 (MQ=255)
              tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGGGCGCGAGGGCCGc                                                                                                                                         >  1:735996/1‑150 (MQ=255)
                      gccgccgccGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                                 <  2:3671058/150‑1 (MQ=255)
                       ccgccgccgCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                                <  2:3628679/150‑1 (MQ=255)
                       ccgccgccgCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                                <  2:3913400/150‑1 (MQ=255)
                            gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                      >  1:414315/1‑139 (MQ=255)
                            gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                           <  2:3437984/150‑5 (MQ=255)
                             ccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGTGGCGCTAGGGCCGCCAATCGGccccccca                                                                                                                          >  1:1766181/1‑146 (MQ=255)
                               gcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTTCGGCGCGATAGGGCCCCAAAGCGGCCCCCACaaa                                                                                                                        >  2:15912/1‑149 (MQ=255)
                                cAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGCGCGGTAGGGCCGCCAAGCGGGCCCGCCaaat                                                                                                                       >  1:286812/1‑148 (MQ=255)
                                cAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGCGGGGCGGAAGGGGCGCAAAGCGGCCCCGACaaaa                                                                                                                       >  1:286816/1‑148 (MQ=255)
                                   gTCGAATCTTGGTACCGGCCGAACTTGCTGGGGGTGGGCTTTGGGAATGGGTAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTAGGGGTAAAGTGTGGGAGCGGCCTTGCGTGGCTTTTGGGCCGCAAAGGGGCCCCGACAATTCaa                                                                                                                    <  1:3484442/150‑1 (MQ=255)
                                   gTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCTTGCGTGGCGGTGGGGGCGCCGGGCGGCCCCGACaagtgat                                                                                                                    >  1:629234/1‑145 (MQ=255)
                                     cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATAAAgc                                                                                                                  >  2:2223705/1‑150 (MQ=255)
                                     cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGGGGGGCCGCCACGCGGCCCCGCCCAGTCAAgc                                                                                                                  >  1:3377350/1‑150 (MQ=255)
                                        aTCGTTGTACCGGCCGAACCTGCTGGCGGTGGGCGTTGGGCTTGGGGAGATGAGGGTTAAAGCCCGATAGAGACCCGATCTACGGGTAATGTGTGGGAGCGGCCTTGCGTCGCTATAGGGCCGAAAAGCGGCCCCGACAATTCAAgcggc                                                                                                               <  1:1282075/150‑1 (MQ=25)
                                        aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCCAAGCGGCCCCCCCACTTCAAgcggc                                                                                                               >  1:3338359/1‑150 (MQ=255)
                                          cGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGa                                                                                                             <  2:332792/150‑1 (MQ=255)
                                           gATGTACGGGCCTAACCTGCTGGCGGTGGGCGTTGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGGGGCCTTGCGTCGCTATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                            <  1:3853847/150‑1 (MQ=255)
                                           gATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                            <  2:2790439/150‑1 (MQ=255)
                                            aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                           <  2:1133288/150‑1 (MQ=255)
                                            aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                           >  2:2639768/1‑150 (MQ=255)
                                            aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                           <  2:2669877/150‑1 (MQ=255)
                                            aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGGGGCCCGGTAGGGGCCCCACGGGGGGCCGACAAGTCGAgcggcccac                                                                                                           >  1:4224297/1‑146 (MQ=255)
                                             tGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgct                                                                                                          <  1:3729696‑M1/150‑41 (MQ=255)
                                              gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                         <  2:664450/150‑1 (MQ=255)
                                                 ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCACCAGTCAAGCGGCGAAGCTGaa                                                                                                      >  1:3747798/1‑150 (MQ=39)
                                                 ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGTCGCGGTTGGGCCGCCACGCGGCCCCGACAATTGAAGCGGCGACGCGGaa                                                                                                      >  1:2546690/1‑150 (MQ=255)
                                                   ggCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAGAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCTAGGGGCGCCAAGCGGCCCCGACAAATCAAGCGGCGCAGCTGagag                                                                                                    >  1:977525/1‑150 (MQ=255)
                                                    gCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGGGCGGTAGGGCCGCAAAGCGGCCCCCCCAAATCAACCGGCGAAGCAGAAAGc                                                                                                   >  1:2329152/1‑150 (MQ=255)
                                                       gacccTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                <  2:3759867/147‑1 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGGGGGCTGGGGCCGCCCAGCGGGCCCCCCAAATGCAGCGGCGAAGGGGGAAGccgg                                                                                                >  1:3067923/1‑148 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGGGCGCTAGGGCCGCCAAGCGGCGCCGCCCAATCAAGCGGCCCAGCGGGAAGccgg                                                                                                >  1:1456918/1‑148 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGCCAATTCAAGCGGCGCAGCTGCAAGccgg                                                                                                >  2:2464607/1‑148 (MQ=39)
                                                        aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  2:2500319/150‑1 (MQ=255)
                                                               aGGCCGTGGGTGATGGGCAAGGGGAGATTACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGAGTGGGAGCGGCCTTGAGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTCGGgcagag                                                                                        >  2:3727420/1‑146 (MQ=38)
                                                               aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgct                                                                                        <  1:1690073‑M1/150‑59 (MQ=255)
                                                               aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                        <  2:3200969/150‑1 (MQ=255)
                                                               aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                        <  2:4269284/150‑1 (MQ=255)
                                                                   cGTGGGCGCTGGGCAAGGGGAGATGACGGTGAACGCCCGCGAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGc                                                                                    <  2:4014102/150‑1 (MQ=255)
                                                                    gTGGGCTCTGGTCATGGGGAGAGGTGGGTGACCGCCCGCTGGAGAACGGCTGGCCGCGGAAATTGTGGTTGCGGCCTTGGTACGCGATAGGGCCGCAAGGCGGGCCCGACATTTCAAGCTGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  2:2099809/150‑1 (MQ=255)
                                                                    gTGGGCGATGGTCAAGGGGATATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCTGAAAGTTTGGGAGCGTCATTGATTCGCGATAGGGCCTCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGACTGGGGCTGCTGTGCa                                                                                   <  2:1912218/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:2861293/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGGGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:839588/150‑1 (MQ=255)
                                                                      gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                 <  2:4043823/150‑1 (MQ=255)
                                                                      gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                 <  2:3953850/150‑1 (MQ=255)
                                                                       ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                <  2:754883/150‑1 (MQ=255)
                                                                           aTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctat                                                                            <  2:412953‑M1/150‑71 (MQ=255)
                                                                            tGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGGGTAAAGTGTGGGTGCGGCCTTGCGCCGCTTTCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                           >  1:510398‑M2/1‑148 (MQ=255)
                                                                            tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatc                                                                           <  1:1846408‑M1/150‑72 (MQ=255)
                                                                            tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                           >  2:1412839‑M2/1‑148 (MQ=255)
                                                                            tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTGGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                           >  2:3435005‑M2/1‑148 (MQ=255)
                                                                             gggCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcg                                                                          <  1:3332763‑M1/150‑73 (MQ=255)
                                                                                 acggggagATGACGGTGACCGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                      >  2:1924432‑M2/3‑143 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCACAATTCAATCGGGGCACCGGAAAGCCAGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                      <  2:2858110‑M2/150‑8 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgac                                                                      <  1:1631946‑M1/150‑77 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                      <  2:2999241‑M2/150‑8 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCAGGGGCTGCGGTGCAGCCCAGtcgcgtc                                                                      <  2:2200852‑M2/150‑8 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCCCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                      <  2:4468727‑M2/150‑8 (MQ=255)
                                                                                    gggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgca                                                                   <  2:2567316‑M2/150‑11 (MQ=255)
                                                                                    gggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCCAAAGGGCCCCAAACCCGCCCCCACCAATCAAGCCGCCAACCGGAAAGCCAGGGGCGGCGGTGCCCCCCAGGCGCCACGcc                                                                   >  2:61050/1‑149 (MQ=255)
                                                                                     ggAGATGACGGTTAAAGCCCGCTCGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                                  >  1:1939108‑M2/1‑139 (MQ=255)
                                                                                     ggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                                  >  2:3539137‑M2/1‑139 (MQ=255)
                                                                                       agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                >  2:366592‑M2/1‑137 (MQ=255)
                                                                                                tGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCGGCCCATAGGGCCGCAAACCCGCCCCCCAAATTAAAGCGGCGGAGCGGAAAGCCAGGGGCTGGGGGGCCGCCCAGTCGCGACCCCCGGCCGCACCca                                                       >  2:122044/1‑150 (MQ=255)
                                                                                                  aaaGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGCGGTCGGGCCGCCACGCGGCCCCCCCCAAGCCAGCGGCGAACCTGAAAGCCGGGGGCTGCGGGGCCGCCCTGTCGCGAGGCCAGGCCGCGCCcacc                                                     >  1:2134723/1‑149 (MQ=255)
                                                                                                    aGCCCGCTAGTTCCCCGCTGTAGGGGTAATGTTTGGGTGGGGCGTTGGGTCGCTTTTGGGCCCAAACGGGGCCCCGAAATTTCACGCGGCTAAGCTGAAAGCCTGGGGCTGATTTGCTGCCCTGTCGCGAGGCAAGGCCGCTCCCACaaa                                                   <  2:642933/150‑1 (MQ=255)
                                                                                                    aGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGCGATAGGGGCGCCAAGCGGCCCCGACAAATCAAGCGGCGGAGCGGAAATCCTGGGGCTGCGGTGGAGCCCGGTCGCGCCGCCAGGCCGCACCCACaaa                                                   >  2:4286274/1‑150 (MQ=255)
                                                                                                                   cGCTCTACGCGGACTGTGTGGGAGCGGCTTTGCTTCGCTATCGGGCCGAACGGCGGCCCCGACACTACACGGGGCTACGCTTAAAGCGGGGGGCTGCTGTGCAGCCCTGGCGCGACGAACGGCCGCTCCCACAAAGACCGCGCAGATCAg                                    <  2:2861997/150‑1 (MQ=255)
                                                                                                                    gCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGGTAGGGCCGCAAAGCGGCCCCGGCAAATCAATCGGCGAAGCGGTAAGCCGGGGGCTGCGGTTCCGCCCGGTCGCCACGCAAGGCCGCACCCACAAAAAACGCCCacaagagt                                   >  2:897886/1‑143 (MQ=19)
                                                                                                                      tctACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCCCCCAATCAAGCGGCGCACCGGAAAGCCGGGGGCTGCGGTGCAGCCCAGTCGCGACCCAAGGCCGCGCCCACCAAGACCCCCCACATCAgggg                                 >  2:2100047/1‑150 (MQ=255)
                                                                                                                        tACGCGGAAAGTGTGGGAGCGGGCTTGCGGCGCGGTAGGGGCCCCAAGCCGCCCCGCCAAATCAAGCGGCGCATCGGCAAGCCTGGGGCGGGTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCCAAGACCGCGCCGAACCGTGGAt                               >  2:4189728/1‑150 (MQ=255)
                                                                                                                        tACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGGGGAGCCCGGTCGCGACGCAAGGCCGCTCCCACAAAGACCCCGCAGATGAGTggct                               >  2:127219/1‑148 (MQ=25)
                                                                                                                         aCGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCCCAAAGCCGCCCCGCCAAATAAAGCGGCGAAGCGGAAAGCCGGGGGCTGCTGTGCAGCCCAGTCGCCACCCAAAGCCGCGCCCACAAACAACGCGCAGATCAGTGGCTc                              >  2:2008962/1‑150 (MQ=255)
                                                                                                                           gcgGAAAGTGTGGGAGCGGCCTTTCGTCGCGGTATGGCCCCCAAGCGGCCCCCACAAATCAAGCGGCGGAGCGGAAAGCCGGGGGGGGCGGTGCCGCCCGGGCGCCGCGCAAGGCCCCTCCCCCAAAGACCCCCCCGAACAGTgggtaac                            >  2:461388/1‑145 (MQ=255)
                                                                                                                             ggAAAGTGTGGGAGCGGCCGTTCGTCGCGAAATGGCCGCAAAGCGGCCCCCACAAATAAAGCCGCGAACCAGAAAGCCAGGGGGGGCTGTGCAGCCCAGTCGCGACGCCAAGCCCCTCCCACAAATACCCCGCAGATCAGTGGATCAacc                          >  2:3392529/1‑150 (MQ=255)
                                                                                                                               tctttgtgggagCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                       <  2:1271088‑M2/99‑7 (MQ=255)
                                                                                                                                     tgGGAGCGGCCATGCGTCGCCATAGGGCCGCAAAGCGGCCCCCACAAATCAAACGGCGAAGCTGAAAGCCAGGGGCTGCTGTGCAGCCCAGTCGCGACCCAAGGCCGCTCCCACAAAGACCCCGCAGATAACTGGATAACCCAGCacatc                  >  2:1803289/1‑150 (MQ=25)
                                                                                                                                     ggggaGGGGCCTGGCGTCGCGCctgggcggcaccgcagcccccggctttccgctccgcggcttgaattttgggggccgctttgcggccctagcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttca                  <  2:121742‑M1/149‑129 (MQ=255)
                                                                                                                                       ggAGCGGCCTTGCGTCGCTATCGGGCCGCAACGCGGCCCCTACAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAc                <  2:1766181/150‑1 (MQ=255)
                                                                                                                                            ggCCTTGCGTCGGGATAGGGCCGCAACGCGGCCCCGACAATTCACGCGGCGAAGCTGAATGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcaca           <  2:3377332/150‑1 (MQ=255)
                                                                                                                                             gCCTTGCGGCGCGATGGGGCCTAAACGCGCCCCCGCAAATACACGCGGCGAAGCTCAACGCCTGGGGCTGCTGTGCCGCCCGGTCGCGAGGAATGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg          <  2:2315795/150‑1 (MQ=255)
                                                                                                                                              ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc         >  2:2637212/1‑150 (MQ=255)
                                                                                                                                              ccTAGCGTCGCGAcagggctgcacagcagccccaggctttcagcatcgccgcttgaattgtcggggccgctatgccgccctaacgcgacgcaaggccgctcccacactttccgcgtagagcgcgtcactaccgggctttcaccgtcatct         >  2:366768‑M1/1‑13 (MQ=255)
                                                                                                                                                ttGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgcgtttgggcccgaaagcgacgcaaggccgcgcccccactttccccgtagagcggggctctagcgggcgttaaccggcaactcc       >  1:1495633‑M1/1‑11 (MQ=255)
                                                                                                                                                ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc       >  2:4245825/1‑150 (MQ=255)
                                                                                                                                                ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc       >  2:4099979/1‑150 (MQ=255)
                                                                                                                                                    gTCGCGAcagggctgcacagcagcccctggctttccggtgccccgattgaattgtgggggccgctttgcggccctatcgcgacgcaagcccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccctt   <  2:2858110‑M1/150‑144 (MQ=255)
                                                                                                                                                    gTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgatttgtcggggccgctttggggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccctt   <  2:4468727‑M1/150‑144 (MQ=255)
                                                                                                                                                    gTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgatttgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccctt   <  2:3343301‑M1/150‑144 (MQ=255)
                                                                                                                                                    gTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccctt   <  2:2999241‑M1/150‑144 (MQ=255)
                                                                                                                                                    gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCAACACCCACAGCCCCAGc   >  1:4320172/1‑150 (MQ=255)
                                                                                                                                                     tCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttg  >  2:2184759‑M1/1‑6 (MQ=255)
                                                                                                                                                     tCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacaaaga                                                 >  1:1271088‑M1/1‑6 (MQ=255)
                                                                                                                                                     tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCa  >  2:2903589/1‑150 (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: