Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,021 G→A 6.1% 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 6.1% 658.3 / 21.0 212E86E (GAG→GAArpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (112/87);  new base A (6/7);  total (118/94)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 216 total aligned reads displayed.

GAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGA  >  NC_002947/541872‑542170
                                                                                                                                                     |                                                                                                                                                     
gAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCgag                                                                                                                                                       >  2:1320561/1‑150 (MQ=255)
gAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCgag                                                                                                                                                       >  2:2003284/1‑150 (MQ=255)
    ttAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGt                                                                                                                                                   >  2:1154073/1‑150 (MQ=255)
             cGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTg                                                                                                                                          <  1:445828/150‑1 (MQ=255)
              ggAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGa                                                                                                                                         <  2:2071839/150‑1 (MQ=255)
               gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                                        <  2:1815823/150‑1 (MQ=255)
               gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTGGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                                        <  1:1070848/150‑1 (MQ=255)
                    gaTCAACTATCGTACGTTCAAGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGtcac                                                                                                                                   <  2:445522/149‑4 (MQ=255)
                    gaTCAACTATCGTACGTTCAAGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGACATCATCTGCGAAAAGTGCGGCGTTGAAGtcac                                                                                                                                   <  1:642205/149‑4 (MQ=255)
                    cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                                   >  1:1685013/1‑150 (MQ=255)
                         aCTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTgg                                                                                                                              >  1:3603219/1‑150 (MQ=255)
                           tACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCg                                                                                                                            >  1:2467225/1‑150 (MQ=255)
                             ccGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGaa                                                                                                                          <  1:2441665/150‑1 (MQ=255)
                             ccGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGaa                                                                                                                          <  1:766503/150‑1 (MQ=255)
                              cGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAg                                                                                                                         <  2:4657163/150‑1 (MQ=255)
                              cGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAg                                                                                                                         >  1:4950282/1‑150 (MQ=255)
                              cGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAg                                                                                                                         >  2:2327387/1‑150 (MQ=255)
                                  cGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtc                                                                                                                     <  1:4766341/150‑1 (MQ=255)
                                  cGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtc                                                                                                                     <  2:4509746/150‑1 (MQ=255)
                                     tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  >  2:4406910/1‑150 (MQ=255)
                                     tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTGGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  <  2:2902761/150‑1 (MQ=255)
                                      cAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCCCCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtcg                                                                                                                 >  1:440127/1‑150 (MQ=255)
                                       aaGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTcgtcgt                                                                                                                <  2:3603219/150‑1 (MQ=255)
                                         gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                                              <  1:1944652/150‑1 (MQ=255)
                                         gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                                              >  1:1883286/1‑150 (MQ=255)
                                           cTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgc                                                                                                            >  1:1977140/1‑150 (MQ=255)
                                            tGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  <  2:2065393/143‑1 (MQ=255)
                                               gCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgc                                                                                                            >  1:3668437/1‑146 (MQ=255)
                                                cGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATg                                                                                                       <  2:5353870/150‑1 (MQ=255)
                                                  tGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                                              <  2:229120/141‑1 (MQ=255)
                                                      ggCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGCGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGcac                                                                                                 >  1:4716649/1‑150 (MQ=255)
                                                       gCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGcaca                                                                                                >  2:1198517/1‑150 (MQ=255)
                                                         cTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATc                                                                                              >  1:4146163/1‑150 (MQ=255)
                                                         cTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATc                                                                                              <  1:3413944/150‑1 (MQ=255)
                                                              cTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCt                                                                                         <  2:1757597/150‑1 (MQ=255)
                                                               tGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                                        <  1:4322093/150‑1 (MQ=255)
                                                                  gcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGcc                                                                                     >  2:857443/1‑150 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  >  2:2634616/1‑119 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                                  <  1:2634616/119‑1 (MQ=255)
                                                                    cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                                   >  2:4325665/1‑150 (MQ=255)
                                                                       gATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGcc                                                                                <  1:2380258/150‑1 (MQ=255)
                                                                       gATCTTGGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGACGTGAGCGCATGGCGCACATCGAGCTGGCCTCGcc                                                                                <  2:3039400/150‑1 (MQ=255)
                                                                         tCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGACGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCgg                                                                              >  2:2528983/1‑150 (MQ=255)
                                                                           tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                                            <  2:547713/150‑1 (MQ=255)
                                                                            ttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTg                                                                           <  2:2566159/150‑1 (MQ=255)
                                                                             tGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCt                                                                                                                                >  2:3929220/1‑96 (MQ=255)
                                                                               gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGccc                                                                        >  2:4203598/1‑150 (MQ=255)
                                                                               gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGccc                                                                        >  1:3027314/1‑150 (MQ=255)
                                                                               gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGc                                                                          <  2:3426912/148‑1 (MQ=255)
                                                                                   aGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAt                                                                    <  2:1406721/150‑1 (MQ=255)
                                                                                      cAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTg                                                                 <  2:1839402/150‑1 (MQ=255)
                                                                                       aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACCTCTgg                                                                >  1:2723120/1‑150 (MQ=255)
                                                                                       aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                >  1:1409332/1‑150 (MQ=255)
                                                                                       aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                <  2:3244329/150‑1 (MQ=255)
                                                                                       aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                >  1:4271731/1‑150 (MQ=255)
                                                                                        aGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGt                                                               >  2:1377240/1‑150 (MQ=255)
                                                                                          gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                             <  2:1907620/150‑1 (MQ=255)
                                                                                          gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                             <  2:5258731/150‑1 (MQ=255)
                                                                                          gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                             >  2:3829967/1‑150 (MQ=255)
                                                                                          gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                             <  1:119190/150‑1 (MQ=255)
                                                                                           aCTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTcc                                                            <  2:1006815/150‑1 (MQ=255)
                                                                                             tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                                          >  1:953933/1‑150 (MQ=255)
                                                                                             tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                                          >  1:2047578/1‑150 (MQ=255)
                                                                                               tgaatGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGaa                                                        >  1:1778378/5‑150 (MQ=255)
                                                                                               cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                        >  2:1267421/1‑150 (MQ=255)
                                                                                               cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                        <  1:4829244/150‑1 (MQ=255)
                                                                                                gaatGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAg                                                       >  1:918944/4‑150 (MQ=255)
                                                                                                    gcttgtgCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGc                                                   >  1:598492/4‑150 (MQ=255)
                                                                                                    gCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGc                                                   <  1:3247122/150‑1 (MQ=255)
                                                                                                      cTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTg                                                 <  1:4406633/150‑1 (MQ=255)
                                                                                                      cTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTg                                                 <  1:1246051/150‑1 (MQ=255)
                                                                                                       tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCt                                                  <  2:1458392/148‑1 (MQ=255)
                                                                                                       tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCt                                                  >  1:1458392/1‑148 (MQ=255)
                                                                                                       tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGc                                                <  2:3772530/150‑1 (MQ=255)
                                                                                                       tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGc                                                >  1:4658095/1‑150 (MQ=255)
                                                                                                            ggCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTcc                                           <  1:4921872/150‑1 (MQ=255)
                                                                                                            ggCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCAGCCGTcc                                           >  2:3165712/1‑150 (MQ=255)
                                                                                                             gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCATATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                          >  1:1268839/1‑150 (MQ=255)
                                                                                                             gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                          >  1:4915232/1‑150 (MQ=255)
                                                                                                               aaGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGt                                        >  2:4714640/1‑150 (MQ=255)
                                                                                                                   aCAAGCGCCTCAAGCACCGCGACGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCg                                    >  1:5416991/1‑150 (MQ=255)
                                                                                                                     aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGc                                                                                          <  2:668984/94‑1 (MQ=255)
                                                                                                                     aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCACTGGCGAAGGTTAGTCGTGAGCGCATGGCGCACATCGAGc                                                                                          >  1:668984/1‑94 (MQ=255)
                                                                                                                      aGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCCAAGGTTCGTCGTGAGCGCATCGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCt                                 >  1:4430963/1‑150 (MQ=255)
                                                                                                                       gcgcCTCAAGCACCGCGGCGTACTCTGCGAGAAGTGCGGCGTTGAAGTGGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCtt                                <  2:1496743/150‑1 (MQ=255)
                                                                                                                         gcCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGc                              >  2:2481673/1‑150 (MQ=255)
                                                                                                                          ccTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCt                             >  1:3857621/1‑150 (MQ=255)
                                                                                                                            tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa                           >  2:3923379/1‑150 (MQ=255)
                                                                                                                            tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAGCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa                           >  2:3303981/1‑150 (MQ=255)
                                                                                                                              aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                         >  2:3390688/1‑150 (MQ=255)
                                                                                                                              aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                         <  1:1822044/150‑1 (MQ=255)
                                                                                                                                gCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGa                       >  1:5084779/1‑150 (MQ=255)
                                                                                                                                 cACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGAc                      <  2:4045201/150‑1 (MQ=255)
                                                                                                                                 cACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGAc                      <  1:4190741/150‑1 (MQ=255)
                                                                                                                                       ggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGAcc                <  1:654012/150‑1 (MQ=255)
                                                                                                                                             aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATcgc                                          <  1:2128623/118‑3 (MQ=255)
                                                                                                                                             aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCtcaag          <  1:3035398/150‑5 (MQ=255)
                                                                                                                                              tCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATAGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTg         >  2:59152/1‑150 (MQ=255)
                                                                                                                                                tGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGat       >  1:1355422/1‑150 (MQ=255)
                                                                                                                                                     gaAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGa  <  1:2219681/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: