Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,146,384 A→G V71A (GTA→GCA)  yceF ← m(7)GTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,3840AG96.3% 70.8 / ‑4.3 27V71A (GTA→GCA) yceFm(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (12/14);  total (12/15)
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.52e-01

CCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTG  >  NC_000913/1146150‑1146623
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
cccGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGCTATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGACCCaat                                                                                                                                                                                                                                 >  2:143955/1‑251 (MQ=255)
  cgcgCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGAGTTCGCCATCAAGAACACATGCCTGGTCTGAACCaataa                                                                                                                                                                                                                               <  1:332050/251‑1 (MQ=255)
       aGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTaa                                                                                                                                                                                                                           >  2:133720/1‑250 (MQ=255)
       aGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTACAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAAGAAATGCCTGGTCTGAACCAATAATTaa                                                                                                                                                                                                                           >  1:340911/1‑250 (MQ=255)
               gggtgatctTTACTCTCGTAATTATAGATCTCCTCCTCGGTCAGTGGAATGAAAGGGACTTCATAAGGCTCCACTTCTGGTTGCGGATGCTCATGCGCAGAATTAACCAACGCCAGGCTAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCATACGGGCATTTTCTTCCGTTAACGGTTTACCAGTGATTTCGCCAGCAAGAACACAGGCCTGGTCTGAGCCAATAACTAAATGTTCCgg                                                                                                                                                                                                                  <  2:149274/245‑1 (MQ=255)
                                                        cAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCCAGAACACATGCCTGGTCTGAACCaat                                                                                                                                                                                                                                 >  2:103690/1‑195 (MQ=255)
                                                        cAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCaat                                                                                                                                                                                                                                 <  1:103690/195‑1 (MQ=255)
                                                                                aGGCTCCACGTATGTTTTCAGATGCCCATCCTACGAATTAAACAACGTCAGTCCAGTCTAGAATGTGTCGATATTGCCGCGTGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCAGCAAGAACACATGCCTGGTCTTATCCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCTCGTTCTCTTGTGCCAGTCGAAGCACCAACTGGCGTGGTGAtt                                                                                                                                                 <  2:298363/251‑1 (MQ=255)
                                                                                    tCCACTTCTGTTTGCAGCTACCCATTCGCCGAATTCAACAACGCCAGTCCAGTATAGAAGGTGACGCTATTGCCGCTGGCTTCGCGTAATTCCAGACGGGCATTTTCTTCCGTAAACGGTTTACCGGTGATTTAGCCATCAAGAACACATGCCTGGACTGAACCACTAATTAAATGATCCGGATAACGTAACGCCagagtt                                                                                                                                                                                               >  2:322506/1‑199 (MQ=255)
                                                                                    tCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGACCACATGCCTGGTCTTAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGAt                                                                                                                                                                                               <  1:322506/201‑1 (MQ=255)
                                                                                     ccACTTCTGTTAGCAGATGCCCATTCGCCGCATTAAACGACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTACTTGCAACCGGGCATTTTCTTCCGTTAACGGTTTAACGGCGATTTCGCCATCAAGGACACGTGCCAGGTCTGCACCAATAATTAAATGATCCGGATAACGGGACGCCCACGATTGCGCTTTTACTTGTGCCAGTAGAAGCACCAACTGGCGGGGTGAtttgt                                                                                                                                              >  2:57156/1‑246 (MQ=255)
                                                                                             gTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCg                                                                                                                                       <  1:22517/248‑1 (MQ=255)
                                                                                                                             cgccatgccagtATTGAAGGTGACGATATTTCCGCTGGATTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCTCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCaa                                                                                                    <  2:47404/244‑1 (MQ=255)
                                                                                                                                         atataaGGTTACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGTGTTTCCGGTTATGTCGCCATCAAGAACGCATGCCTGGTCTGAACCAATAATTTAATAATCCGGAGAACGTGACGCCAGAGATTTCGCTTTTTCTTGTGCCATGCGAAGCACCATCGGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGCCCTCTGGTGCTGCACATTCAAAAGAGATTTGc                                                                                         <  2:384669/246‑1 (MQ=255)
                                                                                                                                                           gCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGCTGCTTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAgcgc                                                                      >  1:38974/1‑251 (MQ=255)
                                                                                                                                                            ccGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGcgcg                                                                     <  1:195650/251‑1 (MQ=255)
                                                                                                                                                                         gTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAg                                                        <  1:57156/251‑1 (MQ=255)
                                                                                                                                                                            tttGCAGACGGGCATTATCTTCCGTAAACGGGTTACCGGAGCTTTCGCCATAAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATACGGATAACGTGGCGCCAGAGATTGCGCTTTTTATTGTGCCAGTCGAAGCACCAACTGTCGCGGAGATTAGTCGATGCTTGGGGTTTCGTCGTCATCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCAAGAAAAGCGCGACGCTATGGAGAGGTg                                                     >  2:24238/2‑251 (MQ=255)
                                                                                                                                                                              tGCAGACGGACATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCCGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAAc                                                                                                                                                              >  2:373896/1‑144 (MQ=255)
                                                                                                                                                                                 aGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGCAGAGCGCGACGCCATGGCGAGGTGGAGGc                                                >  1:274405/1‑251 (MQ=255)
                                                                                                                                                                                  gACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGACAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGGGGCt                                               <  1:112670/251‑1 (MQ=255)
                                                                                                                                                                                     gggCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTaaa                                            <  1:133720/251‑1 (MQ=255)
                                                                                                                                                                                         aTTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATaa                                        <  1:143955/251‑1 (MQ=255)
                                                                                                                                                                                                tcCGTTAACTGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAATAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAg                                  <  1:327733/250‑1 (MQ=255)
                                                                                                                                                                                                    ttAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCatc                             >  1:349604/1‑250 (MQ=255)
                                                                                                                                                                                                                               tCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAATTAAGTTTAGGCTTATTTTTTTCCATCAGATATAGCGTAtt   >  2:79652/1‑251 (MQ=255)
                                                                                                                                                                                                                                cAAGAACTCATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATtg  <  1:112423/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
CCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTG  >  NC_000913/1146150‑1146623

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