Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 3,142,272 C→T 73.7% A74V (GCG→GTG)  PP_2758 → putative Ribose ABC transporter, periplasmic ribose‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154513,142,2720CT73.7% 29.1 / 9.8 19A74V (GCG→GTG) PP_2758putative Ribose ABC transporter, periplasmic ribose‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (6/8);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCACCATCGAGGTGCCGTCGCAGTACGCGGTCAATACGA  >  AE015451/3142066‑3142502
                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
cTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCCTCTACTTCGACAGCCGGGCCTTCAAGGCGCCGGTCGAg                                                                                                                                                                                                                              >  1:32859/1‑217 (MQ=255)
        tcATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCg                                                                                                                                                                                                   <  1:278037/235‑1 (MQ=255)
                               ggCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCaa                                                                                                                                                                            <  1:777056/236‑1 (MQ=255)
                                                                          tGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTg                                                                                                                                          <  2:789263/227‑1 (MQ=255)
                                                                                    cATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATg                                                                                                                       >  1:318020/1‑236 (MQ=255)
                                                                                               gAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTca                                                                                                             <  1:1228856/235‑1 (MQ=255)
                                                                                                                   gCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGcc                                                                                                                            <  1:57896/200‑1 (MQ=255)
                                                                                                                                    ttCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATg                                                                                                                       <  1:201673/188‑1 (MQ=255)
                                                                                                                                     tCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGACGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACa                                                                                                                                                                       >  2:474089/1‑139 (MQ=255)
                                                                                                                                                   gcgcATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACg                                                                                                                                                                                                      <  1:676436/94‑1 (MQ=255)
                                                                                                                                                   gcgcATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACg                                                                                                                                                                                                      >  2:676431/1‑94 (MQ=255)
                                                                                                                                                                      gTGGGCACCAGCATCTACGTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCg                                                                                                                                                                                                   <  2:74094/78‑1 (MQ=255)
                                                                                                                                                                        gggCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCAcc                                    >  1:85407/1‑235 (MQ=255)
                                                                                                                                                                         ggCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGcatc                                                                                       <  1:917047/183‑1 (MQ=255)
                                                                                                                                                                         ggCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGcatc                                                                                       >  2:917042/1‑183 (MQ=255)
                                                                                                                                                                          gCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCAcc                                    >  2:1055394/1‑233 (MQ=255)
                                                                                                                                                                                 cATCTACTTCGACAGCCGTGCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCACCATCGAGGTg                           <  1:608161/235‑1 (MQ=255)
                                                                                                                                                                                                  gggCCTTCAAGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCACCATCGAGGTGCCGTCGCAGTACGCGGt          >  1:602630/1‑235 (MQ=255)
                                                                                                                                                                                                           aGGTGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCACCATCGAGGTGCCGTCGCAGTACGCGGGCCATACgc  >  1:428615/1‑233 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
CTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAACGACAAGGCCATGGTGACCCAGTTGCAGAATCTGGTGACGCTGAAGCCCGATGCAGTTATTCACACCCTCGGCACGCTGAGCATCATCGAGCCCTGGTTCAAGCGCATCAGCGCCGCCGGCATCCCGCTGTTCACCATCGAGGTGCCGTCGCAGTACGCGGTCAATACGA  >  AE015451/3142066‑3142502

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑