Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,069,743 C→T A33T (GCA→ACA)  rutG ← pyrimidine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,069,7430CT96.3% 71.1 / NA 28A33T (GCA→ACA) rutGpyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/10);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGC  >  NC_000913/1069512‑1069978
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
tGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGCCGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGCATCGGCATTAACACCGTCTCGCCAAACACCGCCACCGCGTGTTGAACGCCCCTTACTGTCTTCTGCGCACACGGCAAt                                                                                                                                                                                                                          >  2:461896/1‑251 (MQ=255)
                   cTGACAGTTAAAGCCCGTCGCGGCGATCACCACGCCGACACAAGCCGCGCTGGAGCCGCGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACCTTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCACGACAAACATCGCCACCGCTTGTTGAACGCCCATTACTGTCTGCTGAGCAAACGGCCATCGTTCATCCGCCGCGAccc                                                                                                                                                                                                       >  2:481477/1‑250 (MQ=255)
                                            aTCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGAttt                                                                                                                                                                              >  1:494788/1‑251 (MQ=255)
                                                      cGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACCGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGcc                                                                                                                                                                      >  2:162749/1‑249 (MQ=255)
                                                      cGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGcc                                                                                                                                                                      >  1:396704/1‑249 (MQ=255)
                                                              gCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCACGACATTAAAATGGAAAGATTGGGATCAAGCCCCATCAGCATCGGCATTAACACCCTCGCGCCAAACATCGCCACCGCGCGTTGAACGCCCATTACTGTCGTCTGAGCACACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGaa                                                                                                                                                             >  1:215783/1‑250 (MQ=255)
                                                                                       ggccccgcccgccGGTAATAAAGAAAAACAGTAACGTCCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGACGCGCCAAACATCGCCACCGCTTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCACCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtctc                                                                                                                                   <  1:359294/246‑1 (MQ=255)
                                                                                                          aaaGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCAAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGc                                                                                                                >  1:262349/1‑251 (MQ=255)
                                                                                                            aaaaaaaCAGTACCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCAATACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCTCGTCCACGCCTCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAg                                                                                                              <  1:527573/249‑1 (MQ=255)
                                                                                                                     gTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                      >  1:291485/1‑250 (MQ=255)
                                                                                                                      aaaCGTGCCGATCCCCGACATAAAAATGGAAAGATTGGGAGCCAGCCCCAGCAGAATCGGCATTAACCCCGCCGCGCAAAACATCGCCCCCGCGTGTTGAACGCCCATTACGGTCGTCGGAGCAAAGGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAg                                                                                                    <  2:262523/250‑1 (MQ=255)
                                                                                                                        aCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGt                                                                                                                                      >  2:46624/1‑215 (MQ=255)
                                                                                                                        aCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGt                                                                                                                                      <  1:46624/215‑1 (MQ=255)
                                                                                                                                              aagggAAAGATTGGGACCCAGCCCCATCAGAACCGGCATTAACACCGTCGCGCCAAACATCGCCCCCGCTTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACca                                                                            <  1:408081/248‑1 (MQ=255)
                                                                                                                                                 ggAAAGATGGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCacac                                                                         <  1:240893/251‑1 (MQ=255)
                                                                                                                                                      gATTGGGATCCAGCCACATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCAATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGCTGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGCAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCTAACCACACCAGcc                                                                    >  2:63421/1‑251 (MQ=255)
                                                                                                                                                                    cccATCAGAATCGGCATTAACCCCGTCGCGCAAAAAACCGCCCCCGCGTGTTGAACGCCGATTATGGCCGTCCGAGCAAACGGCAATCTTTCATCCGGCGCGACAACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtg                                                      <  1:63421/251‑1 (MQ=255)
                                                                                                                                                                    cccATCAGAATCGGCATTAACACCGTCCCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtg                                                      >  1:457642/1‑251 (MQ=255)
                                                                                                                                                                       ataaGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTGGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGa                                                   <  1:481477/248‑1 (MQ=255)
                                                                                                                                                                                     tAACCCCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCACCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCtt                                      <  1:199870/250‑1 (MQ=255)
                                                                                                                                                                                     tAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                     >  1:497291/1‑251 (MQ=255)
                                                                                                                                                                                     aaaCACAGTCGCTACAAACCTCGCCCCTGGGTGTGGAACGCCCATGACTGACTTCTGAGAAAAGGGCAACCGTTAATCCGGCGCTGCCACGCCGCTCTCGGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGGAGGGTGTGGTGTGACGATGAATCGCTTc                                     <  2:525706/250‑1 (MQ=255)
                                                                                                                                                                                       acacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGCGAGGAACACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATAAGCGCGTGATAAATGCGATCGCACCACACCAGCCCGTAGGGTGCGGTGTGAAGATGAATCGCTTCg                                    >  2:115849/1‑250 (MQ=255)
                                                                                                                                                                                            gTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGCGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGACATCAGCTGTCTCCTTGAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACAAGCACGTAGGGTGGGGTGTGACGATGAATCGCTTCGATGgcg                              >  1:164018/1‑251 (MQ=255)
                                                                                                                                                                                                    aaactgcGCCACCGCGGGTTGAACGCCCTTTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAg                      <  1:162749/245‑1 (MQ=255)
                                                                                                                                                                                                                        gAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGc  >  1:390210/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGCCCATTACTGTCGTCTGAGCAAACGGCAATCTTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACAGTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCa                                                                                                               <  1:161787/139‑1 (MQ=255)
                                                                                                                                                                                                                           cGCCCATTACTGTCGTCTGAGCAAAAGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCa                                                                                                               >  2:161787/1‑139 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
TGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGC  >  NC_000913/1069512‑1069978

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