Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,707 C→T 7.0% A315V (GCC→GTC)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947542,7070CT 7.0% 568.5 / 21.9 187A315V (GCC→GTC) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (94/79);  new base T (8/5);  total (103/84)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 192 total aligned reads displayed.

GAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGCGTCTGCAC  >  NC_002947/542558‑542855
                                                                                                                                                     |                                                                                                                                                    
gAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCgcgc                                                                                                                                                      >  2:763650/1‑150 (MQ=255)
  aCGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCAGGAAGCTGTTGACGCGCTGTTCGATAACGGCCGTCGTGGCCGCGTCa                                                                                                                                                    >  2:3998874/1‑150 (MQ=255)
    gACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATc                                                                                                                                                  >  1:456227/1‑150 (MQ=255)
     aCCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCa                                                                                                                                                 >  2:2433147/1‑150 (MQ=255)
     aCCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCa                                                                                                                                                 >  2:1240338/1‑150 (MQ=255)
       cTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACt                                                                                                                                               <  1:5055636/150‑1 (MQ=255)
         gTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTgg                                                                                                                                             <  1:917128/150‑1 (MQ=255)
            tCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTc                                                                                                                                          >  1:1269628/1‑150 (MQ=255)
            tCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTc                                                                                                                                          >  2:2800942/1‑150 (MQ=255)
            tCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTc                                                                                                                                          >  2:1590203/1‑150 (MQ=255)
              gCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGa                                                                                                                                        <  2:1792537/150‑1 (MQ=255)
               ccGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGaa                                                                                                                                       <  1:371923/150‑1 (MQ=255)
                cGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAAc                                                                                                                                      <  2:3234331/150‑1 (MQ=255)
                 gggTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACa                                                                                                                                     >  1:2208395/1‑150 (MQ=255)
                  ggTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACaa                                                                                                                                    >  1:1577628/1‑150 (MQ=255)
                  ggTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACaa                                                                                                                                    <  2:1162870/150‑1 (MQ=255)
                  ggTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACaa                                                                                                                                    <  2:3378352/150‑1 (MQ=255)
                   gTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAg                                                                                                                                   >  1:4809203/1‑150 (MQ=255)
                    tGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCAGGAAGCTGTTGACGCGCTGTTCGATAACGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGc                                                                                                                                  >  1:979111/1‑150 (MQ=255)
                    tGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGc                                                                                                                                  <  1:6084/150‑1 (MQ=255)
                         aaCCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCg                                                                                                                             >  2:1430753/1‑150 (MQ=255)
                            cGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTg                                                                                                                          >  2:4468860/1‑150 (MQ=255)
                              tAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGaa                                                                                                                        <  1:1973698/150‑1 (MQ=255)
                                acaacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGt                                                                                                                      >  2:1564614/1‑150 (MQ=255)
                                  aacCGTCTGAGGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTcc                                                                                                                    >  2:1100439/1‑150 (MQ=255)
                                   acCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTccc                                                                                                                   <  1:2549256/150‑1 (MQ=255)
                                    ccGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCt                                                                                                                  <  2:271232/150‑1 (MQ=255)
                                      gTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTgg                                                                                                                >  2:3169342/1‑150 (MQ=255)
                                           aaGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGTTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGAc                                                                                                           <  2:3899643/150‑1 (MQ=255)
                                           aaGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGAc                                                                                                           >  1:3966857/1‑150 (MQ=255)
                                                ccTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGAt                                                                                                      <  2:4984088/150‑1 (MQ=255)
                                                ccTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGAt                                                                                                      <  2:2413291/150‑1 (MQ=255)
                                                 ctgctgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATc                                                                                                     <  2:1102704/150‑1 (MQ=255)
                                                 ctgctgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATc                                                                                                     <  1:880708/150‑1 (MQ=255)
                                                     tgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACCTGATCAAAg                                                                                                 >  1:4486592/1‑150 (MQ=255)
                                                      ggATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAgg                                                                                                >  1:2993191/1‑150 (MQ=255)
                                                         tCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTaa                                                                                             <  2:1691124/150‑1 (MQ=255)
                                                         tCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTaa                                                                                             <  1:26700/150‑1 (MQ=255)
                                                            gTCGGCGCCGGACATCATCGTGCGCAACGAACAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCa                                                                                          <  2:2721679/150‑1 (MQ=255)
                                                             tCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCaa                                                                                         <  1:94785/150‑1 (MQ=255)
                                                              cGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAg                                                                                        <  2:4321648/150‑1 (MQ=255)
                                                              cGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAg                                                                                        <  2:4773486/150‑1 (MQ=255)
                                                                 cgcCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTc                                                                                     <  2:4923959/150‑1 (MQ=255)
                                                                 cgcCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTc                                                                                     <  2:3584330/150‑1 (MQ=255)
                                                                  gcCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCg                                                                                    >  2:3496029/1‑150 (MQ=255)
                                                                    ctggcatcatTATCCGTAACGAAAAGCGTATGTTGCAGGAAGCTGTTGACGCGCTGTTCGATAACGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGCGTCCGCTGAAGTCGCTCTCCGACATGCTCAAGGGTAAGCAGGGTCGtt                                                                                  <  1:2800700/146‑1 (MQ=255)
                                                                    cGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGtt                                                                                  <  2:2148811/150‑1 (MQ=255)
                                                                    cGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGtt                                                                                  <  2:1999871/150‑1 (MQ=255)
                                                                     ggACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGttt                                                                                 <  2:4585546/150‑1 (MQ=255)
                                                                       aCATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTcc                                                                               <  1:1237853/150‑1 (MQ=255)
                                                                       aCATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTcc                                                                               <  1:3433155/150‑1 (MQ=255)
                                                                        catcatCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCg                                                                              <  2:2714708/150‑1 (MQ=255)
                                                                         atcatcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGt                                                                             <  1:4936014/150‑1 (MQ=255)
                                                                          tcatcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTc                                                                            <  2:3036284/150‑1 (MQ=255)
                                                                           catcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCa                                                                           >  2:1434157/1‑150 (MQ=255)
                                                                           catcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCa                                                                           <  2:4840305/150‑1 (MQ=255)
                                                                                tGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACt                                                                      >  2:4667236/1‑150 (MQ=255)
                                                                                  cgcAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTg                                                                    <  1:2716809/150‑1 (MQ=255)
                                                                                   gcAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGc                                                                   >  1:3610690/1‑150 (MQ=255)
                                                                                   gcAACGAAAAGCGCATGCTGCAGGAAGTGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCt                                                                  >  2:1375594/1‑150 (MQ=255)
                                                                                    cAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCAGTCAGAACTTGCt                                                                  >  1:688542/1‑150 (MQ=255)
                                                                                     aaCGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTc                                                                 >  2:4414733/1‑150 (MQ=255)
                                                                                     aaCGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTc                                                                 >  1:491806/1‑150 (MQ=255)
                                                                                      aCGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCg                                                                <  2:1511546/150‑1 (MQ=255)
                                                                                       cGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCgg                                                               <  1:1380230/150‑1 (MQ=255)
                                                                                        gAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGt                                                              >  1:3657042/1‑150 (MQ=255)
                                                                                         aaaaGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTa                                                             >  1:2005332/1‑150 (MQ=255)
                                                                                         aaaaGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTa                                                             >  2:2815089/1‑150 (MQ=255)
                                                                                            aGCGTATGTTGCAGGAAGCTGTTGACGCGCTGTTCGATAACGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGCGTCCGCTGAAGTCGCTCTCCGACATGCTCAAGGGTAAGCAGGGTCGTTTCCGTCAGAACCTTCTCGGTAAGc                                                          <  1:2145666/150‑1 (MQ=255)
                                                                                            aGCGTATGTTGCAGGAAGCTGTTGACGCGCTGTTCGATAACGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGCGTCCGCTGAAGTCGCTCTCCGACATGCTCAAGGGTAAGCAGGGTCGTTTCCGTCAGAACCTTCTCGGTAAGc                                                          <  1:1503954/150‑1 (MQ=255)
                                                                                             gcgcATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCg                                                         >  1:3902433/1‑150 (MQ=255)
                                                                                             gcgcATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCg                                                         >  2:3999725/1‑150 (MQ=255)
                                                                                               gcATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCgtg                                                       >  1:2373414/1‑150 (MQ=255)
                                                                                                cATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCgtgt                                                      <  1:1371336/150‑1 (MQ=255)
                                                                                                 aTGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGtt                                                     <  2:229684/150‑1 (MQ=255)
                                                                                                  tgctgcAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTg                                                    >  2:1085082/1‑150 (MQ=255)
                                                                                                  tgctgcAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTg                                                    >  2:3665400/1‑150 (MQ=255)
                                                                                                    ctgcAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGac                                                  <  1:2251857/150‑1 (MQ=255)
                                                                                                        aGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGactact                                              >  1:3241091/1‑150 (MQ=255)
                                                                                                        aGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGactact                                              >  2:190453/1‑150 (MQ=255)
                                                                                                         ggAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGAc                                                                                                           <  2:4734146/88‑1 (MQ=255)
                                                                                                             gctgTTGACGCGCTGTTCGATAACGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGCGTCCGCTGAAGTCGCTCTCCGACATGCTCAAGGGTAAGCAGGGTCGTTTCCGTCAGAACCTTCTCGGTAAGCGCGTCGACTAGTCgggt                                         <  2:3496649/147‑2 (MQ=255)
                                                                                                                gTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGt                                      >  1:1299251/1‑150 (MQ=255)
                                                                                                                  tGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTc                                    <  2:300402/150‑1 (MQ=255)
                                                                                                                   gACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCg                                   <  2:4809203/150‑1 (MQ=255)
                                                                                                                         ctgctgGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAAtt                             <  2:706653/150‑1 (MQ=255)
                                                                                                                          tgctgGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTa                            >  2:4396926/1‑150 (MQ=255)
                                                                                                                              ggACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGt                        >  2:3697558/1‑150 (MQ=255)
                                                                                                                                ataaCGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGCGTCCGCTGAAGTCGCTCTCCGACATGCTCAAGGGTAAGCAGGGTCGTTTCCGTCAGAACCTTCTCGGTAAGCGCGTCGACTATTCGGGTCGTTCGGt                                 <  1:1832578/137‑1 (MQ=255)
                                                                                                                                ataaCGGCCGTCGTGGCCGCGTCATCACCGGTGCTAACAAGCGTCCGCTGAAGTCGCTCTCCGACATGCTCAAGGGTAAGCAGGGTCGTTTCCGTCAGAACCTTCTCGGTAAGCGCGTCGACTATTCGGGTCGTTCGGt                                 >  2:1832578/3‑139 (MQ=255)
                                                                                                                                  aaaggCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGt                    >  2:1077740/4‑150 (MQ=255)
                                                                                                                                  aaCGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGt                    >  2:1586695/1‑150 (MQ=255)
                                                                                                                                     ggCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCg                 >  2:34746/1‑150 (MQ=255)
                                                                                                                                      gCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGGAATTACCGTAGGTCCGa                >  1:53246/1‑150 (MQ=255)
                                                                                                                                            gtGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGc          <  1:1762264/150‑1 (MQ=255)
                                                                                                                                            gtGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGc          <  2:261988/150‑1 (MQ=255)
                                                                                                                                               gTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGCGTc       >  1:4169702/1‑150 (MQ=255)
                                                                                                                                                tCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGCGTCt      <  2:3535329/150‑1 (MQ=255)
                                                                                                                                                  gcgcCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGCGTCTGc    >  1:874812/1‑150 (MQ=255)
                                                                                                                                                    gcCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGCGTCTGCAc  <  2:4119141/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCGTAGGTCCGACCCTGCGTCTGCAC  >  NC_002947/542558‑542855

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