Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,875,975 G→A L119L (CTG→TTG)  cysD ← sulfate adenylyltransferase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,875,9750GA93.8% 40.2 / NA 16L119L (CTG→TTG) cysDsulfate adenylyltransferase subunit 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/10);  minor base C (1/0);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCAGCCAGATGTATTGCCAGATATCCTGCTCGGTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATA  >  NC_000913/2875729‑2876176
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                         
ccagccagATGTATTGCCAGATATCCTGCTCGGTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTTGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACAGTATTTGTTCAGCGCCTGTT‑TCAAGcc                                                                                                                                                                                                        <  1:4051/250‑1 (MQ=255)
                                gTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCATTTGTTAATTTGCCCGTTGTAGTTGCGACACAGCTCCGGGGGCTGATTTTTCGGATCACAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTAACGGCGCGCACAACCGACGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATAGCGGTATGTTTCGCGCt                                                                                                                                                                       >  1:217947/1‑251 (MQ=255)
                                                                                tGCCCGTTGTAGTTGTGCCACAGCTCGGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATGCGCTCTGTAGCGCGGGATTGCTCTTCGTCACGGCGCGGGCCAGCGAATGCGGCATCAAAACCGTATTTGTGCAGCGCC‑GTTCTCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCGCAAAGTGATTAATCCCCATCGCCACGCCTTCCGGGtttt                                                                                                                        <  1:196686/250‑1 (MQ=255)
                                                                                              gtgCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTAACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCCTCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATcc                                                                                                                                                >  1:227724/1‑212 (MQ=255)
                                                                                                gCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCg                                                                                                       <  1:400544/251‑1 (MQ=255)
                                                                                                           gggCGCTGATTTTTCGGATCCCAGCGATGGAAACGGCCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCt                                                                                                                                                                       >  1:159773/1‑176 (MQ=255)
                                                                                                                cTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGGTTAATCCCCATCGcc                                                                                                                                        <  1:171425/202‑1 (MQ=255)
                                                                                                                         cGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTGCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATATGTTCAGCGCCCGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTcgcg                                                                                                                                                                         <  1:157449/160‑1 (MQ=255)
                                                                                                                                                    naaaGAGTAAATTCGATCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTcgcg                                                                                                                                                                         >  1:131774/2‑133 (MQ=255)
                                                                                                                                                      aaGAGTAAATTCGATCTGTAGCGCGGGATTTCTCTTCGTCACGGCGCGCCCCACCGAAGGCGGCAGCAAAACCGTATTTGTTCAGCGCCGGTT‑TCAAGCCTTCGGTTTTCATAATATCGGTATGTTTCGCTCTGCCGTGCACGAATTGATTAATCCCCATCGCCACGCCTTCCGGGttn                                                                                                                         <  1:57683/179‑2 (MQ=255)
                                                                                                                                                                         aGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGcc                                                                                                                                                                                                        <  1:177589/81‑1 (MQ=255)
                                                                                                                                                                            gcgGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAgn                                                                                                   <  1:96445/179‑2 (MQ=255)
                                                                                                                                                                                                ggCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCa                                                                                                     <  1:35545/157‑1 (MQ=255)
                                                                                                                                                                                                      caccacCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCACGCCTTCAGTTTTCCTAATATCGGTCTGTTTCGCGCTGCCCTGCACGACTGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTACGATCGCGGAACCCATACATCTCGCGGCATTTCCAGCCGGTATCGACATGCCGCAACGGGCAAGGCAGCGTACCTGGATa  >  1:245004/1‑250 (MQ=255)
                                                                                                                                                                                                      caccacCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAAGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAgg                 >  1:401126/1‑235 (MQ=255)
                                                                                                                                                                                                                    gcATCAAAACCGTATTTGGTCAGCGCCTGTG‑TCAAGCCTTCAGTTGTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCggn                                                                                                                            <  1:58424/114‑2 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                         
CCAGCCAGATGTATTGCCAGATATCCTGCTCGGTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTT‑TCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATA  >  NC_000913/2875729‑2876176

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