Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,150,418 A→G intergenic (‑153/+29) eptC ← / ← ppc phosphoethanolamine transferase EptC/phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,150,4180AG87.5% 29.9 / ‑2.4 17intergenic (‑153/+29)eptC/ppcphosphoethanolamine transferase EptC/phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (8/6);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

CCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGA  >  NC_000913/4150191‑4150653
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
ccAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAACCAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGCAGTAT‑TTAAGAAAACCCTCGCGCAGACGCACGAGGGCTTGC‑agaag                                                                                                                                                                                                                        >  2:112593/1‑249 (MQ=255)
  aaaaGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAgg                                                                                                                                                                                                                     >  1:76301/1‑250 (MQ=255)
           gCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGcc                                                                                                                                                                                                           <  1:160223/251‑1 (MQ=255)
                                                                                      gTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTTAAGGTATATTCAGAATTTGAAAAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCt                                                                                                                                 <  1:201010/250‑1 (MQ=255)
                                                                                                         gAATGTAACGACCATTCCTTAAGGATATCTGAAGGTATATTAAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGATTGA‑AGAAGAGGAAGATTAGCCGATATTACGCATACCTGCCGCAATCCCGGCAATAGTGACAATTAAAGCTTGTTCGACGCGAGGATCCGATTCCTGGCCTTCTTTTTCTGCCTgg                                                                                                             >  1:237174/1‑251 (MQ=255)
                                                                                                         gAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTgg                                                                                                             >  1:274628/1‑251 (MQ=255)
                                                                                                               aactaataTTCCTTAAGTAAATCTGTAGGTATATTCAGAATTTGAATAAAATGCAGTGAAAAATATATATAAAACG‑AGGGGTTTAGAACAGGAGTAT‑TGCAAAAAACCCGCGCGCAGAAGCACGAGG‑TTTGCCAGATGAGGAAGATTAGCCGTTATTACGCATACCTGCCGCAATCCCGGCAATAGTGTCCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAg                                                                                                       <  2:57727/244‑1 (MQ=255)
                                                                                                                           ttAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttct                                                                                                                         >  1:10868/1‑221 (MQ=255)
                                                                                                                             aaGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGGTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGTGCGTTGCAGCACCTc                                                                                         >  2:128520/1‑251 (MQ=255)
                                                                                                                               ggATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCgg                                                                                       >  1:123458/1‑251 (MQ=255)
                                                                                                                                      tGAATGTATATTCAAAAGTTGAAAAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAATTAT‑TTCAGAAACCCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTGTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCa                                                                                 <  1:75733/250‑1 (MQ=255)
                                                                                                                                        aaGGTATATGCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATa                                                                              <  1:242222/251‑1 (MQ=255)
                                                                                                                                                                      gacagaAATAAATTGAAAAGG‑AGGTTGTAAAAACAGAA‑TATATGCAAAAAACCCACCGGTAAAAGCACGGGGTTTTGC‑AAAAGATGAAAATTACCCGTTATTACGCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTGCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTa                                                <  2:274628/251‑1 (MQ=255)
                                                                                                                                                                          aaaaTTTATTGAAA‑CGTAGGGTGTTAGAACAGAATTAT‑TTCAAAAAACCCTCGCGCATAAGCAGGAGTTTTTGC‑AAAAGAGGAAGATTAGCCGGTATTACGCATACCGGCCGCAAGCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTTCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATCCGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTg                                            <  1:191502/250‑1 (MQ=255)
                                                                                                                                                                                    gAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAACCCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGActct                                  <  1:139874/251‑1 (MQ=255)
                                                                                                                                                                                        aCG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCAGACCTGCCGCAATACCGGCCATAGTGACCATTAACGCTTGTGCGACGCGAGGCTCCGGTTCCTGGCCTTCTTTTTCTTCCTGGCGGGAACGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGGAAATATTCCGTCACTGAATAGCATCTGCa                               >  1:220171/1‑250 (MQ=255)
                                                                                                                                                                                                                      gAAAACCCTCGCGCAGAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAACTCTTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGGTCGGCCATAAGATgggg  >  2:293763/1‑250 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
CCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACG‑AGGGTGTTAGAACAGAAGTAT‑TTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGC‑AGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGA  >  NC_000913/4150191‑4150653

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