Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,392,725 C→T intergenic (+67/‑167) glyY → / → yjeV tRNA‑Gly/uncharacterized protein YjeV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,392,7250CT83.3% 43.7 / 6.3 24intergenic (+67/‑167)glyY/yjeVtRNA‑Gly/uncharacterized protein YjeV
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (9/11);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

ACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCAC  >  NC_000913/4392497‑4392938
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                     
aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGAAAGGCACAGAGCACCCAAGCGGGAAAAGCTCAGTTGGTAGAGCACGAACTTGACACGGTCGGGGACGCGAGTTCGAGTCTGGATTCCCGCACCAAAATCTTATCGCAATAAAATACACACAGCGGCATGATGTGTTATTGTTGGTTTTTGTTGTCTCTGATAAAATTTTGTAAACaggtgttt                                                                                                                                                                                                 >  2:127660/1‑245 (MQ=9)
       tgcaaaGGTCGGGGTCGCGATTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCcaca                                                                                                                                                                                            <  1:40573/245‑1 (MQ=255)
             ggTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGCGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAgg                                                                                                                                                                                    >  1:289180/1‑251 (MQ=255)
               tCGGGGTCGCTAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACa                                                                                                                                                                                                        <  2:74791/229‑1 (MQ=255)
               tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACa                                                                                                                                                                                                        >  1:74791/1‑229 (MQ=255)
                        cgAGTTCGAGTCTCTTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTa                                                                                                                                                                         <  2:289415/251‑1 (MQ=255)
                              cGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTc                                                                                                                                                                    >  1:313165/1‑250 (MQ=255)
                                                acAAAATTTGAAAGTGCTGTAAGACACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTg                                                                                                                                                 <  2:13390/250‑1 (MQ=255)
                                                 cAAAATTTGAAATTACGGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCCCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGTTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTGCGCACGGTGa                                                                                                                                                <  2:30967/251‑1 (MQ=37)
                                                  aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGa                                                                                                                                                >  1:253817/1‑250 (MQ=255)
                                                  aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCACGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGa                                                                                                                                                >  2:292613/1‑250 (MQ=255)
                                                        tgaaattgcgggaatgcccaGAACACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTa                                                                                                                                          <  2:289180/232‑1 (MQ=38)
                                                                    aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCt                                                                                                                              <  2:253817/250‑1 (MQ=255)
                                                                     aagcaaaGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGACAAGGTCGGGGGCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCAAAGCGACGCGATGCGTTAGTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATTAAAAGCCGCAAGCTGTAATCTTTATTTCACAGAAACTTCGCCCGAAGAATAGCATTTTTTTAACCTAt                                                                                                                            >  2:188096/7‑251 (MQ=25)
                                                                      ggCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTAtt                                                                                                                           <  2:319690/251‑1 (MQ=255)
                                                                                        gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAAt                                                                                                         <  2:98906/251‑1 (MQ=255)
                                                                                           aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGGCGCGAGCTCACGTCACGTTTCACGCTCCAAATTCTTCTCTCAATAAAATATCCAAAGCGACGCGATGCGTTATTGCTGGTTCTTGTTGAGTCTGATTAACTCTTGTAAACAGAGATATCCACCGCATCAGGCTGGAATCTTAATATCAAAGAAACTTAGCACGGTGAATAGTATTTATTTAACCTTTTGATAGTGAAGTTAACAATTa                                                                                                       >  2:74760/1‑250 (MQ=255)
                                                                                                                      tGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCTAGTCGATCACt                                                                            <  1:103532/250‑1 (MQ=255)
                                                                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAACATATCCACAGCGCCGCGATGCGTTATTGCTGGTTTTTGTTGTATCTGATAAACTCTTGTAACCAGAGTTATCCACAGCCTAAGGATGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGTGTCGATCACTTGAc                                                                        >  2:100294/1‑251 (MQ=255)
                                                                                                                                                     gTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCg                                            <  2:221736/251‑1 (MQ=255)
                                                                                                                                                                cAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGtttttattttt                                 <  1:168098/251‑1 (MQ=255)
                                                                                                                                                                          tctcAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGATAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTGGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGc                       >  2:174122/1‑251 (MQ=255)
                                                                                                                                                                                                                              ctctGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCAc  >  2:86629/1‑220 (MQ=255)
                                                                                                                                                                                                                              ctctGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTAGTTTTATTCACAAGCTGTGGATGAATCAGGCGTCAc  <  1:86629/220‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                     
ACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCAC  >  NC_000913/4392497‑4392938

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