Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,021 G→A 5.0% E86E (GAG→GAA 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,0210GA 5.0% 1089.6 / 28.1 338E86E (GAG→GAArpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (177/144);  new base A (10/7);  total (187/151)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 344 total aligned reads displayed.

GAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGA  >  NC_002947/541872‑542170
                                                                                                                                                     |                                                                                                                                                     
gAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCgag                                                                                                                                                       <  2:1960458/150‑1 (MQ=255)
    ttAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGt                                                                                                                                                   >  1:3757850/1‑150 (MQ=255)
     tAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTg                                                                                                                                                  <  1:5580391/150‑1 (MQ=255)
      aagaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgc                                                                                                                                                 <  1:6267414/150‑1 (MQ=255)
       agaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGTGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcg                                                                                                                                                >  2:9024233/1‑150 (MQ=255)
        gaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcgg                                                                                                                                               <  1:1739076/150‑1 (MQ=255)
        gaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcgg                                                                                                                                               >  2:3066080/1‑150 (MQ=255)
            ccGGAGACGATCAACTATCGTACGTTCAAGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGtt                                                                                                                                           <  1:9080514/150‑1 (MQ=255)
              ggAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGa                                                                                                                                         >  1:2390148/1‑150 (MQ=255)
               gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCGTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                                        <  1:3140645/150‑1 (MQ=255)
               gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                                        <  2:2400151/150‑1 (MQ=255)
                aaaCCATCAACTACCGTACGGTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAg                                                                                                                                       <  2:4673927/150‑1 (MQ=255)
                    cATCAACTACCGTACGTTCACGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                                   <  2:6023780/150‑1 (MQ=255)
                     aTCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGcc                                                                                                                                  <  1:8030553/150‑1 (MQ=255)
                      tCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGccc                                                                                                                                 >  2:625829/1‑150 (MQ=255)
                      tCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGccc                                                                                                                                 >  2:7527952/1‑150 (MQ=255)
                        aaCTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTg                                                                                                                               >  1:1308446/1‑150 (MQ=255)
                         aCTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTgg                                                                                                                              >  1:642533/1‑150 (MQ=255)
                          cTACCGTACGTTCCAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTGGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGc                                                                                                                             <  2:6692795/150‑1 (MQ=255)
                          cTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGc                                                                                                                             >  2:3023623/1‑150 (MQ=255)
                           tACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCg                                                                                                                            >  2:8642294/1‑150 (MQ=255)
                           tACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCg                                                                                                                            >  2:6852330/1‑150 (MQ=255)
                                tACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGt                                                                                                                       >  1:4106134/1‑150 (MQ=255)
                                  cGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtc                                                                                                                     <  2:160913/150‑1 (MQ=255)
                                   gTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcg                                                                                                                    >  1:2592435/1‑150 (MQ=255)
                                   gTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcg                                                                                                                    <  2:2075322/150‑1 (MQ=255)
                                     tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  <  1:275795/150‑1 (MQ=255)
                                     tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  >  1:7270904/1‑150 (MQ=255)
                                     tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  >  2:7091528/1‑150 (MQ=255)
                                       aaGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTcgtcgt                                                                                                                <  2:4495908/150‑1 (MQ=255)
                                       aaGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTcgtcgt                                                                                                                >  2:2796812/1‑150 (MQ=255)
                                         gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                                              >  2:4038808/1‑150 (MQ=255)
                                          ccTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCTTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAg                                                                                                             >  1:3740235/1‑150 (MQ=255)
                                             gAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgcgc                                                                                                          >  1:3634823/1‑150 (MQ=255)
                                                cGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATg                                                                                                       <  2:1436780/150‑1 (MQ=255)
                                                  tGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGgc                                                                                                     >  2:5799/1‑150 (MQ=255)
                                                   gACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGgcg                                                                                                    <  1:8485659/150‑1 (MQ=255)
                                                    aCGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGgcgc                                                                                                   <  2:7789098/150‑1 (MQ=255)
                                                      ggCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGcac                                                                                                 >  1:6566830/1‑150 (MQ=255)
                                                          tGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCg                                                                                             >  2:122764/1‑150 (MQ=255)
                                                            ttCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAg                                                                                           <  1:5738925/150‑1 (MQ=255)
                                                              tgggccaaaGATCTTGGCGCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                                        >  2:5931238/8‑150 (MQ=255)
                                                              tgggccaaaGATCTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                                        >  2:1703068/8‑150 (MQ=255)
                                                               tGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                                        <  2:2522734/150‑1 (MQ=255)
                                                               tGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                                        <  1:1066477/150‑1 (MQ=255)
                                                                gcgcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTgg                                                                                       <  2:5974892/150‑1 (MQ=255)
                                                                  gcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCTCCTCAAGCACCGAGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCACTGGCGAAGGTTCGTCGTGAGAGCATGGCGCACATCGAGCTGGcc                                                                                     >  2:5159400/1‑150 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                                   >  2:52543/1‑150 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                                   >  2:1007712/1‑150 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                                   >  1:4783064/1‑150 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                                   >  2:4113034/1‑150 (MQ=255)
                                                                     aaGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCg                                                                                  >  1:8379556/1‑150 (MQ=255)
                                                                     aaGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCg                                                                                  >  1:1024267/1‑150 (MQ=255)
                                                                     aaGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCg                                                                                  >  1:1481218/1‑150 (MQ=255)
                                                                      aGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGc                                                                                 <  2:3357010/150‑1 (MQ=255)
                                                                         tCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCgg                                                                              <  2:5650843/150‑1 (MQ=255)
                                                                          cTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGt                                                                             <  2:8887865/150‑1 (MQ=255)
                                                                           tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                                            <  1:7738079/150‑1 (MQ=255)
                                                                            ttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTg                                                                           <  2:3757841/150‑1 (MQ=255)
                                                                              ggCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGcc                                                                         <  2:7324636/150‑1 (MQ=255)
                                                                              ggCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGcc                                                                         <  2:3785151/150‑1 (MQ=255)
                                                                               gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGccc                                                                        >  1:1629131/1‑150 (MQ=255)
                                                                                cccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCca                                                                       >  1:7106717/1‑150 (MQ=255)
                                                                                cccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCca                                                                       >  1:2800531/1‑150 (MQ=255)
                                                                                cccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCca                                                                       >  1:7553885/1‑150 (MQ=255)
                                                                                 ccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCcac                                                                      >  1:5007420/1‑150 (MQ=255)
                                                                                 ccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCcac                                                                      >  1:7473557/1‑150 (MQ=255)
                                                                                   aGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAGGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCGGGTTGccctccg                                                                    >  1:7006829/1‑146 (MQ=255)
                                                                                    gTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATc                                                                   <  2:2689349/150‑1 (MQ=255)
                                                                                        aGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGt                                                               >  1:9105470/1‑150 (MQ=255)
                                                                                          gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                             <  1:5812703/150‑1 (MQ=255)
                                                                                                     ccTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCt                                                  <  2:3114334/150‑1 (MQ=255)
                                                                                                        gtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGcc                                               >  1:7445976/1‑150 (MQ=255)
                                                                                                             gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                        >  2:1761930/1‑136 (MQ=255)
                                                                                                             gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                          <  2:8410031/150‑1 (MQ=255)
                                                                                                                aGTACTAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATAGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCTTCCCGTa                                       <  1:3013232/150‑1 (MQ=255)
                                                                                                                 gTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTAt                                      >  2:1308376/1‑150 (MQ=255)
                                                                                                                     aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                                  >  1:3238549/1‑150 (MQ=255)
                                                                                                                     aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGGTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                                  <  1:6867309/150‑1 (MQ=255)
                                                                                                                      aGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCt                                 <  1:6928030/150‑1 (MQ=255)
                                                                                                                            tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTc                                            >  2:2316948/1‑133 (MQ=255)
                                                                                                                             cAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGAt                          <  1:1964023/150‑1 (MQ=255)
                                                                                                                             cAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGAt                          <  1:2279665/150‑1 (MQ=255)
                                                                                                                             cAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGAt                          <  1:3065943/150‑1 (MQ=255)
                                                                                                                             cAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGAt                          <  2:6549770/150‑1 (MQ=255)
                                                                                                                              aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                         >  1:5166564/1‑150 (MQ=255)
                                                                                                                                gCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGa                       <  1:7795173/150‑1 (MQ=255)
                                                                                                                                 cACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGAc                      <  2:3853871/150‑1 (MQ=255)
                                                                                                                                  aCCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACa                     >  2:5421535/1‑150 (MQ=255)
                                                                                                                                   ccGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACAt                    <  1:3955990/150‑1 (MQ=255)
                                                                                                                                    cgcgGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATg                   >  1:4894108/1‑150 (MQ=255)
                                                                                                                                     gcggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa                  <  1:8649552/150‑1 (MQ=255)
                                                                                                                                     gcggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa                  <  1:5982706/150‑1 (MQ=255)
                                                                                                                                          gTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTg             <  1:466586/150‑1 (MQ=255)
                                                                                                                                           tAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGc            <  2:7490306/150‑1 (MQ=255)
                                                                                                                                             aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCtcaag          <  1:6651336/150‑5 (MQ=255)
                                                                                                                                             aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCtcaag          <  1:617254/150‑5 (MQ=255)
                                                                                                                                                tGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGat       >  2:7236192/1‑150 (MQ=255)
                                                                                                                                                tGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCtcaaggat       >  2:7988880/1‑143 (MQ=255)
                                                                                                                                                   gagaAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATc    <  2:8189437/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGA  >  NC_002947/541872‑542170

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: