Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,098,620 T→A V245V (GTT→GTA ghrA → glyoxylate/hydroxypyruvate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,098,6200TA90.0% 47.3 / ‑0.5 21V245V (GTT→GTAghrAglyoxylate/hydroxypyruvate reductase A
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (9/9);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CGAACCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGGAAT  >  NC_000913/1098384‑1098860
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
cGAACCCGTAAATCGTGGCCTGTCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCGTAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATTACCTGCTCTCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGt                                                                                                                                                                                                                                      <  2:412658/250‑1 (MQ=255)
  aaCCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGGGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTCACAGGCGCAATGTTgg                                                                                                                                                                                                                                   >  1:386540/1‑251 (MQ=255)
         aaagcGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTAATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGtttt                                                                                                                                                                                                                            <  2:204259/247‑1 (MQ=255)
         aaaTCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACACCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGtttt                                                                                                                                                                                                                            >  1:383320/1‑251 (MQ=255)
                    gggCGTGGAAAGCTGTGCCGGACGGGAAGAACTGGCTGCATTGCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCGCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAAc                                                                                                                                                                                                                <  2:259416/250‑1 (MQ=255)
                       cGTGCAAAGATTTGCCGGACGGGACGAACTGTCTGCATTTCTGAGCCAATGTCGGGGATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAAAAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGATGTTCTAGTTGGGGAAGATGACCTGCTCTCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGGTGTTTTTAATCGTGAACCCt                                                                                                                                                                                                             >  2:397909/1‑251 (MQ=255)
                                         aCGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGGTGACCCTTACCGCCTGAAAGTCCGc                                                                                                                                                                                           >  2:179732/1‑251 (MQ=255)
                                            ggAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAAGTTGTTACCGAAAACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCGTTCCGCCTGAAAGTCTGctct                                                                                                                                                                                        <  2:364420/251‑1 (MQ=255)
                                                                            gggTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACAt                                                                                                                                                        <  1:221850/251‑1 (MQ=255)
                                                                                               aCCGAATACCCCTGAAACCGGCGGCATTATTAATCAACAATTACTCGAAAAATTCCCGGATGGCGCGTATCTCCTCAACCGGGCGCGTGGGGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGTCGATAACACCACATGTCGCCGCGATTACCCGTc                                                                                                                                     <  2:323341/251‑1 (MQ=255)
                                                                                                           tGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTAGCGCGAGGTGTTCATGTTGTGGAAGATGAACTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACACTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGCAGCtg                                                                                                                         >  1:178040/1‑251 (MQ=255)
                                                                                                                     ggCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACAt                                                                                                               >  1:313534/1‑251 (MQ=255)
                                                                                                                                                aaaTTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGaaaa                                                                                    <  1:179732/251‑1 (MQ=255)
                                                                                                                                                 aaTTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACCCCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTTGCACCATTGCCCAGCTCGaaaa                                                                                    >  2:268727/1‑250 (MQ=255)
                                                                                                                                                    tACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGCTAACACCACATGTCGCCGCGATTACCCGTCCCGCTGACGCTGTGGAGTCCATTTCTCGCACCATTGCCCAGCCCGAAAAgggg                                                                                >  2:137869/1‑251 (MQ=255)
                                                                                                                                                                   aTCTCCTAAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGATCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCggg                                                                 >  2:38699/1‑251 (MQ=255)
                                                                                                                                                                                   gcgtgtttttCATGGTGTGGAAGATGACCTGCTCGCGGCGCTGGGTGGCGGCAAAGTAAAAGGAGCAATTTTGGTGGTTGTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCAc                                                 <  2:144277/243‑1 (MQ=255)
                                                                                                                                                                                                                         cgcTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGGCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCtt           <  2:383320/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGGAAt  >  1:35994/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATGACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCtt           >  2:190756/1‑235 (MQ=255)
                                                                                                                                                                                                                                         aGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATGACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCtt           <  1:190756/235‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
CGAACCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGGAAT  >  NC_000913/1098384‑1098860

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