Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,472,103 G→A T98I (ACT→ATT)  pyrB ← aspartate carbamoyltransferase catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,1030GA88.9% 17.7 / ‑5.8 9T98I (ACT→ATT) pyrBaspartate carbamoyltransferase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (4/4);  total (4/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAG  >  NC_000913/4471889‑4472240
                                                                                                                                                                                                                       |                                                                                                                                         
ggTGCGGCCATATTGCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAAGAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGGTGGTTGGAGCCATCACCGGCATTCAGGACCTGTACATTGCCGGAAAACTCGGTGGCCAGTCGCGCCGCACCTTCCTGCTTATGACGCATCTCTATCGCATCGACGTAAAGTCTGATAACCGAAATGGTATCGGCCAGCGTTTCGc                                                                                                       <  1:227338/251‑1 (MQ=255)
        cATATTTCAGGTCACCCACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGAAGTTGGCTGGAGCCATCCCCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAATGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCttttt                                                                                                >  1:83942/1‑250 (MQ=255)
                                                                         ntGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAATGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGt                              >  1:63892/2‑251 (MQ=255)
                                                                                              tGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCGCTATCGCATGGACGTAAATGCGGATATCCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTg                                                                                <  1:99891/180‑1 (MQ=255)
                                                                                                 gTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGAGCGCCGCACCTTCCTGCGGATGACGCATCACTATCCAATCGACGTAAATGCTg                                                                                                                                      >  1:128944/1‑123 (MQ=255)
                                                                                                     gATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGTCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAATGCTGATAACCGAAATGGTATCGGCCAGCGTGTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCAGTGCATAGATGTTTCGAAAGAGAGTCGGTTACGGGTagng  <  1:197995/251‑1 (MQ=255)
                                                                                                                                                                gCCAGGCGCGCCGCACCTTCCTGCNGNTGACGCATCACTATCGCATCGACGTAAATGCTGATAACCGAAATGGTANCGGCCAGCGTNTCGNNCTTTNTNCNCAGTGATGTATTGGCGCTGTCGGAGAAGCNCACCACGCTGGCCCCCAGGCGGTNNATNNNTGNNTCGaa                        <  2:191304/170‑1 (MQ=255)
                                                                                                                                                                gCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAATGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACAACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGaa                        >  1:191304/1‑170 (MQ=255)
                                                                                                                                                                                                cATCACTATCGCATCGACGTAAATGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAg                                                                 <  1:117086/97‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                         
GGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAG  >  NC_000913/4471889‑4472240

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