Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,910,165 T→C D174G (GAC→GGC)  mazG ← nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,910,1650TC95.5% 66.7 / ‑4.5 22D174G (GAC→GGC) mazGnucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (12/9);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAT  >  NC_000913/2910031‑2910300
                                                                                                                                      |                                                                                                                                       
cTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCCCCTCGCCGAtc                                                                                                                                    <  1:588040/140‑1 (MQ=255)
cTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGAtc                                                                                                                                    >  1:615675/1‑140 (MQ=255)
     tGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTc                                                                                                                               <  2:484107/140‑1 (MQ=255)
        aaTGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTa                                                                                                                            <  1:631658/140‑1 (MQ=255)
               tttCTGCTTTCGTCCCTAAATGGCGAGCCAGTTTAACCGGGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTa                                                                                                                     <  2:294436/140‑1 (MQ=255)
               tttCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAACCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTa                                                                                                                     <  2:834110/140‑1 (MQ=255)
                                       gAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCgtcgtc                                                                                             <  2:86074/140‑1 (MQ=255)
                                          ccAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCaa                                                                                          <  2:192386/140‑1 (MQ=255)
                                                  aCCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTATTCGTAGAATTTATCGGCTACCGGACCAAGCGTCGTCCAATCGAAGcc                                                                                  >  2:345943/1‑140 (MQ=255)
                                                        gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGtt                                                                            >  2:480826/1‑140 (MQ=255)
                                                                        ccATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCTAACGTTAGCGCAACGTTTCTgg                                                            >  1:416046/1‑140 (MQ=255)
                                                                                 cctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTgc                                                   >  2:593794/1‑140 (MQ=255)
                                                                                 cctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTATTCGTAGACTTTATCGACTACCGGACCAAGAGTCGTCCAGTCGAAGCCAACGTTGGAGCAATGTTTCTGGATTTTTTgc                                                   >  2:321633/1‑140 (MQ=255)
                                                                                    cgcgTTTACCGTGGTCGCCACTAGTCTGCCGCTTTGCGAACATCACCTCGCCGCTCTCTACGTAGACTTTATCTGTTTGCGGACCAATCTCCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTTGAATTTTTGcgca                                                <  2:535173/137‑1 (MQ=255)
                                                                                    ccAGTTTAGCCTGGTGGACAACCGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGTCTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGcgca                                                <  1:345943/140‑1 (MQ=255)
                                                                                                ggTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGcc                                    <  2:615675/140‑1 (MQ=255)
                                                                                                 gTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCg                                   >  1:698690/1‑140 (MQ=255)
                                                                                                        acaGCCTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACt                            <  2:708175/140‑1 (MQ=255)
                                                                                                            ccTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGa                        >  1:447490/1‑140 (MQ=255)
                                                                                                            ccTGCCGCGCTTCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGa                        >  1:442064/1‑140 (MQ=255)
                                                                                                                      ttCGTACATCACCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGt              >  1:160164/1‑140 (MQ=255)
                                                                                                                                aCCTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGa    >  1:392108/1‑140 (MQ=255)
                                                                                                                                  cTCGCCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAt  >  2:715467/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAT  >  NC_000913/2910031‑2910300

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