Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,155,549 C→A 34.0% T183T (ACC→ACA argC → N‑acetylglutamylphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,155,5490CA34.0% ‑0.3 / 4.8 12T183T (ACC→ACAargCN‑acetylglutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base A (2/2);  total (4/8)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

GGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACG  >  NC_000913/4155408‑4155676
                                                                                                                                             |                                                                                                                               
ggAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGCTCTTCTTGACGTCAATCAGTCGACGGTGATCAACGccacca                                                                                                                                >  1:212366/1‑143 (MQ=255)
  aGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTAGCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAgc                                                                                                                              <  1:139335/143‑1 (MQ=255)
                       aaggagCGAATGTGATGGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACTTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGt                                                                                                         <  2:78975/143‑1 (MQ=255)
                             cGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCCACGCCACAAGCGGCGTGAGCGGTGCTGGGCGTCAAGCg                                                                                                   >  1:82021/1‑143 (MQ=255)
                             cGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACAAGCGGCGTGAGCGGTGCAGGGCGTAAAGCg                                                                                                   >  1:364344/1‑143 (MQ=255)
                                                             acgacgGCACAGCTGGCGCTGAAACCATTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAACGACCAATTCAAACAGCTGTTGTGAAGTTAGCCt                                                                   >  2:93140/1‑143 (MQ=255)
                                                                                       gTTGATTTATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATGGGGGGCGCCACCCACCAGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATc                                         <  1:361906/143‑1 (MQ=255)
                                                                                        ttgattgaTGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCg                                        >  1:68727/1‑143 (MQ=255)
                                                                                                            tcaCCCCAACCAGCAGCCGGTGATCAACCCCACCAGCGGCGGGAGCGTTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCca                    <  2:133524/141‑1 (MQ=255)
                                                                                                              agctctATCAGTGGCCGGTGATCAACGCCACAAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCcaca                  <  2:82021/141‑1 (MQ=255)
                                                                                                              aCCTCAATCAGTGGCCGGTGATCAACGCCACAAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCcaca                  <  2:364344/143‑1 (MQ=255)
                                                                                                                            ccGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGa    <  2:34364/143‑1 (MQ=255)
                                                                                                                              ggTCATCAGCGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGCTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACg  <  2:196452/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                               
GGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACG  >  NC_000913/4155408‑4155676

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