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→GTA) ghrAglyoxylate/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

CGAACCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGGAAT  >  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)
                    gggCGTGGAAAGCTGTGCCGGACGGGAAGAACTGGCTGCATTGCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCGCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAAc                                                                                                                                                                                                                <  2:259416/250‑1 (MQ=255)
                       cGTGCAAAGATTTGCCGGACGGGACGAACTGTCTGCATTTCTGAGCCAATGTCGGGGATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAAAAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGATGTTCTAGTTGGGGAAGATGACCTGCTCTCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGGTGTTTTTAATCG‑TGAACCCt                                                                                                                                                                                                             >  2:397909/1‑251 (MQ=255)
                                         aCGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGGTGA‑CCCTTACCGCCTGAAAGTCCGc                                                                                                                                                                                           >  2:179732/1‑251 (MQ=255)
                                            ggAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAAGTTGTTACCGAAAACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCGTTCCGCCTGAAAGTCTGctct                                                                                                                                                                                        <  2:364420/251‑1 (MQ=255)
                                                                            gggTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACAt                                                                                                                                                        <  1:221850/251‑1 (MQ=255)
                                                                                               aCCGAATACCCCTGAAACCGGCGGCATTATTAATCAACAATTACTCGAAAAATTCCCGGATGGCGCGTATCTCCTCAACCGGGCGCGTGGGGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGTCGATAACACCACATGTCGCCGCGATTACCCGTc                                                                                                                                     <  2:323341/251‑1 (MQ=255)
                                                                                                           tGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTAGCGCGAGGTGTTCATGTTGTGGAAGATGAACTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACACTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGCAGCtg                                                                                                                         >  1:178040/1‑251 (MQ=255)
                                                                                                                     ggCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACAt                                                                                                               >  1:313534/1‑251 (MQ=255)
                                                                                                                                                aaaTTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGaaaa                                                                                    <  1:179732/251‑1 (MQ=255)
                                                                                                                                                 aaTTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACCCCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTTGCACCATTGCCCAGCTCGaaaa                                                                                    >  2:268727/1‑250 (MQ=255)
                                                                                                                                                    tACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGCTAACACCACATGTCGCCGCGATTACCCGTCCCGCTGACGCTGTGGAGTCCATTTCTCGCACCATTGCCCAGCCCGAAAAgggg                                                                                >  2:137869/1‑251 (MQ=255)
                                                                                                                                                                   aTCTCCTAAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGATCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCggg                                                                 >  2:38699/1‑251 (MQ=255)
                                                                                                                                                                                   gcgtgtttttCATGGTGTGGAAGATGACCTGCTCGCGGCGCTGGGTGGCGGCAAAGTAAAAGGAGCAATTTTGGTGGTTGTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCAc                                                 <  2:144277/243‑1 (MQ=255)
                                                                                                                                                                                                                         cgcTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGGCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCtt           <  2:383320/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGGAAt  >  1:35994/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aGTTAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATGACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCtt           >  2:190756/1‑235 (MQ=255)
                                                                                                                                                                                                                                         aGTTAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATGACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCtt           <  1:190756/235‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
CGAACCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCG‑TGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGGAAT  >  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: ‑