Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,971,099 A→G 13.8% noncoding (85/86 nt) PP_mr25 → gyrA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,971,0990AG13.8% 66.0 / 3.6 29noncoding (85/86 nt)PP_mr25gyrA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (8/17);  new base G (2/2);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTTCCAAGACCGCCCATCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGTGACGGCTTGAAGCCTGTG  >  NC_002947/1970954‑1971246
                                                                                                                                                 |                                                                                                                                                   
tttCCAAGACCGCCCATCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGaa                                                                                                                                                 <  1:264027/150‑1 (MQ=255)
 ttCCAAGACCGCCCATCACGGCGGCCTTCAATGCGCAGATCAATAGAAAAAAAAAATGATTTGTAGTAAGTCGTAGCACCCTTATGCACTGCGGCGGCGAGCTACGTTCGGAAATTGATTGAGCTGCTAAATTTCAACAAAAAAAGGAAt                                                                                                                                                <  2:352253/150‑1 (MQ=255)
  tCCAAGACCGCCCATCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGAGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATc                                                                                                                                               <  1:213739/150‑1 (MQ=255)
         ccGCCCATCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTtc                                                                                                                                        <  1:104170/150‑1 (MQ=255)
              cATCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAGGGAATCAGGCTTCTCATg                                                                                                                                   >  1:759811/1‑150 (MQ=255)
               aTCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGGCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATgg                                                                                                                                  <  2:319843/150‑1 (MQ=255)
                tCACGGCGGCATTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTCCGTCTGATAACAAGATAGCGCTGCAACGGATAGCTTCGTTCGGCCCTTGATGGCGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATggg                                                                                                                                 <  1:390852/150‑1 (MQ=255)
                       ggcCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTg                                                                                                                          <  2:750509/150‑1 (MQ=255)
                        gcCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTgg                                                                                                                         >  2:948612/1‑150 (MQ=255)
                              aTTGCCCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAg                                                                                                                   <  1:356478/150‑1 (MQ=255)
                                   gcAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAAtc                                                                                                              >  1:531075/1‑150 (MQ=255)
                                       aTCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAAtgctcc                                                                                                          <  2:715181/150‑1 (MQ=255)
                                           aTGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGt                                                                                                      >  2:1083511/1‑150 (MQ=255)
                                            tGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTTACCAGAAAAGGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTc                                                                                                     >  2:194856/1‑150 (MQ=255)
                                            tGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATg                                                                                                                                   >  1:46360/1‑120 (MQ=255)
                                            tGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATg                                                                                                                                   <  2:46360/120‑1 (MQ=255)
                                                acaACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAata                                                                                                 <  1:166243/150‑1 (MQ=255)
                                                           aTTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGaa                                                                                      <  2:292111/150‑1 (MQ=255)
                                                                   tAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCtt                                                                                                                                         >  1:11141/1‑91 (MQ=255)
                                                                   tAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCtt                                                                                                                                         <  2:11141/91‑1 (MQ=255)
                                                                             cACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTc                                                                    <  1:1083511/150‑1 (MQ=255)
                                                                                                         cgttcgGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAGGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGc                                        <  2:759811/150‑1 (MQ=255)
                                                                                                                cccTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTc                                                                    >  2:414425/1‑115 (MQ=255)
                                                                                                                cccTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTc                                                                    <  1:414425/115‑1 (MQ=255)
                                                                                                                       tGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGAt                          <  1:948612/150‑1 (MQ=255)
                                                                                                                             tGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGt                    <  2:531075/150‑1 (MQ=255)
                                                                                                                                    tttCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGTGACGGCt             >  2:1045862/1‑150 (MQ=255)
                                                                                                                                    tttCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGTGACGGCt             >  1:1127147/1‑150 (MQ=255)
                                                                                                                                     ttCACCAGAAAAGGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGTGACGGCtt            <  1:194856/150‑1 (MQ=255)
                                                                                                                                             aaaaaGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGt                    >  1:65502/1‑134 (MQ=255)
                                                                                                                                             aaaaaGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGt                    <  2:65502/134‑1 (MQ=255)
                                                                                                                                               aaaGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGTGACGACTTGAAGCCtgtg  >  1:848989/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTTCCAAGACCGCCCATCACGGCGGCCTTCATTGCGCAGATCCATGGCACAACTCATTGATTTGTCGTAAGTCGTCGCACCCTTATGCGCTGCGGCGGCGAGCTTCGTTCGGCCCTTGATGGAGCTGCGAAGTTTCACCAGAAAAAGGAATCAGGCTTCTCATGGGCGAACTGGCCAAAGAAATCCTCCCGGTCAATATCGAAGACGAACTGAGACAGTCTTACCTCGACTACGCGATGAGCGTGATTGTCGGGCGAGCGCTGCCCGATGCGCGTGACGGCTTGAAGCCTGTG  >  NC_002947/1970954‑1971246

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