Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,908,319 C→G D450H (GAC→CAC)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,908,3190CG84.6% 16.9 / ‑5.0 13D450H (GAC→CAC) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (4/7);  minor base T (0/2);  total (4/9)
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.39e-01

AAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCG  >  NC_000913/2908072‑2908500
                                                                                                                                                                                                                                                       |                                                                                                                                                                                     
aaGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTgaac                                                                                                                                                                                    <  1:169249/251‑1 (MQ=255)
                      cATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGGAATTAGGAACTTCGATGATCTCGACCAACTGATCATCCCAGGAACAGCCCGCAACGCGCAGACATGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGAAGCTCAACAATTCTCGGCGCATTGTACAGCTGGAGAACAAGGCTATCGTGAACCAACTGGCACTGCTGTGCGCCg                                                                                                                                                              >  1:139207/1‑251 (MQ=255)
                                      aGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGGTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTGAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGTAccgc                                                                                                                                              <  1:228728/251‑1 (MQ=255)
                                                        cTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTGAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGc                                                                                                                            >  1:5656/1‑251 (MQ=255)
                                                                      cccatgggTGATTCGGAACTTCGATGATCTCGACCAACTGTTCATCCCCGGAACGGCCCTCAACGCGCAGCCCTCCTTCTTCAATCTTTTTCAACAGCTTTTTTTTGACTTCGTCACGGTGCCGTTTCCGCTCACCAATTTCCGGCGCATTCCCCATCTGGCGAACCCGGCTATCGTGAACCATCTGGCCCTGCTGTGTGCCGATCCGCATGGTTCCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACg                                                                                                              <  1:277807/246‑1 (MQ=255)
                                                                                aTTCGGAACTTCGATGATCTCGACCAACTGATCAGCCCCGGAACGGCCCGCACCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGGCGGCGCATTGTACAGCTGGCTAACCAGCCTATCGTGAACCAACTGGCACTGCTGTGCGCCGAGACGCAGGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTgn                                                                                                          <  1:50436/245‑2 (MQ=255)
                                                                                          gcGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTTTTTCAACAGCATGTTTTTGACTTCGTAACGGTGACGAGGACTCTCAACAATGGTCGGGGCATTGTGCATCGGGCGAACCAGGCTATCGTTAACCAACTGGCACTTCTGTGCTCCGGGACGCATGGTACCGCCGAGATCGCTCTTGTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGcc                                                                                           <  1:133684/249‑1 (MQ=255)
                                                                                            natgatCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTGAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGTACCGCCGAGATCGc                                                                                                                                     >  1:25367/2‑206 (MQ=255)
                                                                                                                          aCGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGCCTATCGTGAACCAACAGGCACTTCTGTGCGCCGAGACGCATTTGACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAgcg                                                                            <  1:189828/233‑1 (MQ=255)
                                                                                                                                                      tCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTTAACCAACTGGCACTGCGGTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTgg                               <  1:254022/250‑1 (MQ=255)
                                                                                                                                                             gTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTGAACCAACTGTCACTGCTGTGCGCCGAGACGCAAGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGTTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTtctc                        <  1:283558/250‑1 (MQ=255)
                                                                                                                                                                                  nttCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTGAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCg  >  1:182326/2‑251 (MQ=255)
                                                                                                                                                                                                      cGCTCAACAATTGTCGGCGCATTGTACAGCTGTCGAACCAGGCTATCTTGAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGt                                                                                                                                                   <  1:271339/86‑1 (MQ=255)
                                                                                                                                                                                                                            tNNACAGCTNNCGAACCAGGCTATCGTGAACAAACTGGCACTGCTGTGCGCCGGGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCaa                                                                                                                <  2:5656/99‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                     
AAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCACGGGTGATTCGGAACTTCGATGATCTCGACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAACCAGGCTATCGTCAACCAACTGGCACTGCTGTGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCG  >  NC_000913/2908072‑2908500

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