Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,485,318 C→T A753V (GCG→GTG)  hrpA → ATP‑dependent 3'‑>5' RNA helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,485,3180CT90.9% 15.9 / ‑3.4 12A753V (GCG→GTG) hrpAATP‑dependent 3'‑>5' RNA helicase HrpA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (5/5);  total (6/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

AGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGCGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTT  >  NC_000913/1485082‑1485549
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
aGAATCCAGCCGCCTGTGGGGGCGCATTGCTGCTCGTATCTACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCATGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTGTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGTGCTTGTGGAAGGt                                                                                                                                                                                                                            <  1:68600/250‑1 (MQ=255)
                            nctgcGCGTATCGACCCTGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGTGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCg                                                                                                                                                                                               >  1:84367/2‑251 (MQ=255)
                                                   ggTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCAAGATCGATCCGGCGTGAGGTCGTGAACTCTTTATTCGCCACGTGCTGGTGGAAgg                                                                                                                                                                                                                             <  1:280141/198‑1 (MQ=255)
                                                           cAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCCGGGCGCGGTGATGGCAACGGAAAAAGTAACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGCGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAAc                                                                                                                                                                                 >  1:310770/1‑234 (MQ=255)
                                                                                                   gCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGTGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTtg                                                                                                                        <  1:265445/251‑1 (MQ=255)
                                                                                                                gcgcAGGGCGCGGTGATGGCAACGGAACAAGTCACTGCTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGTGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGaaa                                                                                                                                                                                          >  1:267964/1‑172 (MQ=255)
                                                                                                                         gcgGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGTGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGtt                                                                                                  >  1:240233/1‑251 (MQ=255)
                                                                                                                                                                                             cAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATGCGCCACGTGCTGGGGGAAGGTGACTGGCAGACGCGGCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTgg                                                                                                                                                         <  1:145094/128‑1 (MQ=255)
                                                                                                                                                                                                    atcgatcCGGCGTTATGTCGTGAACTCTTTATGCGCCACGGGCTGGTGGAAGGTGACTGGCAGACGCTTCACGCTTTCTTCCGTGAAAGCCTGAAACTAGGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTcgccgc                                                                                                                                     <  1:151637/141‑1 (MQ=255)
                                                                                                                                                                                                    atcgancCNGNGTNNTNCCGTNAACTCTTTATTCGCCACGTGCTGGTGGAAGGTGACTggn                                                                                                                                                                                                                     >  2:151637/1‑60 (MQ=255)
                                                                                                                                                                                                                         gAACTCTTTATTCTCCACGTGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGATCACAAAGCACATCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCTCATCAGCCACGATTTAATCTCCGCTCGCCACGTCGACAGCTGTTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACtt  <  1:159703/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccACGTGCTGGTGGAAGGTGAGTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACCGCTGGTGGa                                       >  1:260180/1‑200 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
AGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTTGCCGCGCGCAAGGTCAACTACAGCCAGATCGATCCGGCGTTATGTCGTGAACTCTTTATTCGCCACGCGCTGGTGGAAGGTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTT  >  NC_000913/1485082‑1485549

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