Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,130,210 G→A N218N (AAC→AAT flgA ← flagellar basal body P‑ring formation protein FlgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,2100GA96.2% 65.4 / ‑6.9 26N218N (AAC→AATflgAflagellar basal body P‑ring formation protein FlgA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/8);  minor base T (0/1);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAAC  >  NC_000913/1130038‑1130455
                                                                                                                                                                            |                                                                                                                                                                                                                                                     
ccgccCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGGGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCtgt                                                                                                                                                                         <  1:24014/251‑1 (MQ=255)
           ttCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCACGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGc                                                                                                                                                              >  2:7217/1‑251 (MQ=255)
                      cgcgTCAGTGGTTTCGCGCGTGTGAACGGTGATTACAGGCTTCAGAGTCGAAGTGCGATCAATACTCTTGGTTTATTCCGCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGTTAATCTTTAACAGCTTACAGATTTAAAAGAATATGCCCACCTGCATCAACAACGCCGCTGCCTACCTGTCCCGATACCATGCGCACCCGCGCATTTTGTGCGGCGGCTGCTTTGTTCAgct                                                                                                                                                   <  2:87523/251‑2 (MQ=255)
                        cgTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTg                                                                                                                                                               <  1:224723/237‑1 (MQ=255)
                        cgTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTg                                                                                                                                                               >  2:224723/1‑237 (MQ=255)
                             gTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTCCAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCCTTGTTCAGCCCCTgccc                                                                                                                                            >  2:219455/1‑248 (MQ=255)
                                                       tACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCa                                                                                                                   <  1:200786/250‑1 (MQ=255)
                                                        aCAGGCTTCAGAGGCGAAGTGCGAGCAATACTCATGGTTTATTCCTCATTGAGGGCGATTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCAcc                                                                                                                 <  2:137963/251‑1 (MQ=255)
                                                                          gTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGAc                                                                                               >  2:345659/1‑251 (MQ=255)
                                                                                 cAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCCTGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGt                                                                                        >  2:49749/1‑251 (MQ=255)
                                                                                  aaTACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTc                                                                                       >  1:264893/1‑251 (MQ=255)
                                                                                    tACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTc                                                                                                                                         <  1:331935/199‑1 (MQ=255)
                                                                                    tACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTc                                                                                                                                         >  2:331935/1‑199 (MQ=255)
                                                                                      cTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGc                                                                                   >  1:292021/1‑251 (MQ=255)
                                                                                         aTGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGCTCAAATTGACGCGTTGTCCCGCttt                                                                                >  1:211215/1‑251 (MQ=255)
                                                                                                    cTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCagg                                                                     >  2:91306/1‑249 (MQ=255)
                                                                                                              gcgcTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTACTCTTTAACAGCTTACAGATTTATAAGAATAATCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGaa                                                            >  2:216453/1‑250 (MQ=255)
                                                                                                                      atcattttttGCTTATTTATCGGCAAGGGACGGGTAATTTTTAACAGCTTAAAGATTTAAAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTTGGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCACCCCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGtt                                                    <  1:405710/242‑1 (MQ=255)
                                                                                                                          tgttttGCTTATTTATCGGCAAGGGACGGTAAATCTTAAACAGCTTACAGATTTATAAGAATTTCCCCACCTGCATCAACACCGCCGTTGCCTCCCTGCCCCGATACCATGCGCCCCCGCGCATGTTGTGCGACGGCTGCATTGTTCAGCGCTTGCCCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGAGGGTTTGTCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAAc                                                 <  1:186632/245‑1 (MQ=255)
                                                                                                                           gtgtTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTTACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTgg                                              >  2:80676/1‑251 (MQ=255)
                                                                                                                           gtgtTGCTTATTTACCGGCAAGGGACGGGTAATCTTTACCACTTTACAGTTTTAAAAGAATATTCCCACCGGCACCACCAACGCCGCTGACTCCCTGCCCCGATACCATGCGCCCCCGCGCATTCTGTGCGACGGCGGCATTGTTCAGCGCCTGCCCTTCTGCGTGGGCGCTAACCCCATCCCCGCTGGCGATCACATTGACGCGTTGTCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTgg                                              <  1:216453/251‑1 (MQ=255)
                                                                                                                                           ggCAAGGGACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTg                              >  1:192076/1‑251 (MQ=255)
                                                                                                                                                 ggACGGGTAATCTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGAt                            >  1:273659/1‑247 (MQ=255)
                                                                                                                                                            cTTTAACAGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGACGCTGACTACCTGTCCCGATACCCTGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGCACCTGGGTTAACTGGATAGGTTGATCGGGTGCTAGATCGCGCCGGCt              >  2:13353/1‑250 (MQ=255)
                                                                                                                                                                   aGCTTACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCAt      >  1:391295/1‑251 (MQ=255)
                                                                                                                                                                       tACAGATTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaac  >  1:32928/1‑251 (MQ=255)
                                                                                                                                                                            |                                                                                                                                                                                                                                                     
CCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAAC  >  NC_000913/1130038‑1130455

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