Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 669,952 G→A N207N (AAC→AAT nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913669,9520GA86.7% 25.1 / ‑2.4 15N207N (AAC→AATnadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (7/6);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGGTTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGCTGAAGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTG  >  NC_000913/669720‑670193
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
cgcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACAACATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGATTAATGTAAGTCAGTAcc                                                                                                                                                                                                                                 >  1:98285/1‑251 (MQ=255)
 gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACTCCATGTAGCAGCGTATGCAGATTTTGCGCTGGCGCAATACCGCGCGCGGTCAGGGGGGTGTGCGCATGACCGCTGTAAAGACCATCGATATGCTCTTGCGTGTCACCATGACTAAGTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCATCGATACAAGCCTTGTTGATTAATGTAAGTCAGTACCg                                                                                                                                                                                                                                <  1:241800/251‑1 (MQ=255)
                 ccGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACTCCAGGTAGCGGCGTATGCAGATTTTGCGCTTGCTCAATGCCGCTCGGGGTCAGGGGGGTGGGCGCATGACTGCTGTAAAGACGATCGATAGGCGCGTGCGTTGCACCATGACGAGGTAACCACAGTCGCAAGATGCCCTCCGTAACGACAGGTATCAGCTATACAAGCCTTGTGGATTAATGTAAGTCAGTACCGGTTCCGGCAATAAATc                                                                                                                                                                                                                <  1:187346/251‑1 (MQ=255)
                                                                                                     cgcgGTCAGGGTGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGGTTCACCATTACGAAGTAACCACAGTCGCATCAGGCGCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGATTAATGTAATTCAGTACCg                                                                                                                                                                                                                                <  1:241809/151‑1 (MQ=255)
                                                                                                     cgcgGNNAGNNNGGNNGGCGCATGACCGCTGTAAAGACCNTCGATATTCGCTTGCGTTTCAACATGACGNANTNACCANNGTCNCATCATGCCCNCCGTAACGACAGGTATCAGCGATCCAAGCCTTGTTGATTCATGTAAGTNANTCCCg                                                                                                                                                                                                                                >  2:241809/1‑151 (MQ=255)
                                                                                                              ggggggggGCGCATGACCGCTGGAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGGCGCATCATGCCCGCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGATTAATGTAAGTCAGTCCCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGTTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAgg                                                                                                                   <  1:27811/251‑1 (MQ=255)
                                                                                                               ngggTGGGCGCATGACCGCTGTAAAGACCATCGCTATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGGTTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGc                                                                                                                  >  1:200287/2‑251 (MQ=255)
                                                                                                               ggggTGGGCGCATGACCGCTGTACAGACCATCGATATTCGCTTGCGTTTCCCCCTGACGAATTAACCCCCGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGATTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTc                                                                                                                                                                                                    >  1:135683/1‑169 (MQ=255)
                                                                                                                       cgcATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACCGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGCTACAAGCCTTGTTGATTAATGTAAGTCAGTACCGGTTCCGGCAATa                                                                                                                                                                                                                    >  1:146773/1‑145 (MQ=255)
                                                                                                                                             tCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGATTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGt                                                                                                                                                                                               <  1:252035/144‑1 (MQ=255)
                                                                                                                                              cGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGATTAATGTAAGACAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGCTGAAGGTGAAGATCTTCCGGGTTATGTGTc                                                                                    <  1:318260/250‑1 (MQ=255)
                                                                                                                                                                  ccATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGGTTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGCTGAAGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGc                                                               >  1:88070/1‑251 (MQ=255)
                                                                                                                                                                                               gCCCTCCGTAACGACAGGTATCAGCGATCCAAGCCTTGTTGATTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGCTGAAGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTgg                                  >  1:273906/1‑251 (MQ=255)
                                                                                                                                                                                                               ggTATCAGCGATACAAGCCTTGTTGATTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGCCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGCTGAAGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCg                   >  1:296342/1‑250 (MQ=255)
                                                                                                                                                                                                                               gCCTTGTTGATTAATGTAAGTCAGTACCGGTTACGGCACTAAATCCTCACATGATTTACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCCGATAAATTTTACCGGCAGGCTGACGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGGGGGTAACCTGGACGCCGACAGACGATCAAAtgtg  >  1:312794/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggTTANTGTAANTNNGTCCCGGTTCTGGCAATAANTCCTCACATGATTCANCGTTNTGCAANCGnt                                                                                                                                                                                   >  2:131099/1‑66 (MQ=255)
                                                                                                                                                                                                                                       ggTTAATGTAAGTCAGTACCGGTTCCGGCAATAAAGCCTCACAGGATTCACCGTTTTGCAAACGTTCGCGGATGAGGGTCGCCGAGATGTTAAACCACGGCGTTTCAGcc                                                                                                                                       <  1:131099/110‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
CGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTATCAGCGATACAAGCCTTGTTGGTTAATGTAAGTCAGTACCGGTTCCGGCAATAAATCCTCACATGATTCACCGTTTTGCAAACGTTCGCGGATGATGGTCGCCGAGATGTTAAACCACGGCGTTTCAGCCAGATAAATTTTACCGGCAGGCTGAAGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTG  >  NC_000913/669720‑670193

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