Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,975,345 G→T 20.0% intergenic (‑109/+285) fliC ← / ← PP_4379 flagellin/3‑oxoacyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,975,3450GT20.0% 49.1 / 2.5 20intergenic (‑109/+285)fliC/PP_4379flagellin/3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base T (2/2);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTTGGTGTTAACGGTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTG  >  NC_002947/4975213‑4975477
                                                                                                                                    |                                                                                                                                    
gTTGGTGTTAACGGTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGtt                                                                                                                                   >  1:114824/1‑136 (MQ=255)
            ggTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCaaaaaaa                                                                                                                       >  1:519472/1‑136 (MQ=255)
                aaaGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAAt                                                                                                                   <  1:120646/136‑1 (MQ=255)
                                      tGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGGTTGATGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTgg                                                                                             >  1:71777/1‑135 (MQ=255)
                                      tGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTga                                                                                             >  1:545520/1‑136 (MQ=255)
                                           ggATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAAt                                                                                        <  1:178243/136‑1 (MQ=255)
                                                     gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTGGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATg                                                                                                          <  2:1041/108‑1 (MQ=255)
                                                     gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTa                                                                              >  1:211171/1‑136 (MQ=255)
                                                     gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTa                                                                              >  1:473862/1‑136 (MQ=255)
                                                     gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATg                                                                                                          >  1:1041/1‑108 (MQ=255)
                                                                                   ccTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGttt                                                >  1:311625/1‑136 (MQ=255)
                                                                                     aagatAACCTTAGTCATGGTTTTCGTCGGGGGCGCGGGTTGCTTGTGGGTTATTTCAAATAAAATATTCTGGCATGGGGCCTGGCCTGAGAAATCAAGGGGTTTATGTCGGTTTATTGGCGTCATAAACTTTTTGa                                              <  2:405706/135‑1 (MQ=255)
                                                                                                     tGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTTGGCTGAATGCCt                              <  2:114824/136‑1 (MQ=255)
                                                                                                          aTCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCa                                                          >  2:511836/1‑103 (MQ=255)
                                                                                                          aTCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCa                                                          <  1:511836/103‑1 (MQ=255)
                                                                                                          aTCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATcc                         >  1:138878/1‑136 (MQ=255)
                                                                                                                        gggTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCgg           <  2:472492/136‑1 (MQ=255)
                                                                                                                             gCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAg      <  1:447023/136‑1 (MQ=255)
                                                                                                                                 gAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTg  <  2:211171/136‑1 (MQ=255)
                                                                                                                                 gAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTg  <  2:473862/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GTTGGTGTTAACGGTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTG  >  NC_002947/4975213‑4975477

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