Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,595 G→A 23.1% intergenic (‑890/+21) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,5950GA23.1% 52.0 / 12.2 26intergenic (‑890/+21)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (14/6);  new base A (3/3);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

GCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATA  >  NC_000913/3708482‑3708743
                                                                                                                 |                                                                                                                                                    
gCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGa                                                                                                                                                  <  1:1041828/118‑1 (MQ=255)
gCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGa                                                                                                                                                  >  2:1041828/1‑118 (MQ=255)
        gcctgataagacgcgtcagcgtcgcatcaggcctacgTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGc                                                                                                                                                >  2:887340‑M2/38‑112 (MQ=255)
        gcctgataagacgcgtcagcgtcgcatcaggcctacgTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGc                                                                                                                                                <  1:887340‑M2/75‑1 (MQ=255)
                         gcAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGc                                                                                                                                            <  3:485555/99‑1 (MQ=17)
                         gcAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGc                                                                                                                                            >  4:485555/1‑99 (MQ=17)
                                   caggcctacgTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAgg                                                                                                                                             <  1:1762091‑M2/78‑1 (MQ=255)
                                   caggcctacgTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAgg                                                                                                                                             >  2:1762091‑M2/11‑88 (MQ=255)
                                    tcgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                    <  1:1075174/80‑1 (MQ=255)
                                    tcgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                    >  2:1075174/1‑80 (MQ=255)
                                    tcgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAgg                                                                                                                                             <  4:1752605/87‑1 (MQ=255)
                                    tcgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAgg                                                                                                                                             >  3:1752605/1‑87 (MQ=255)
                                           cgTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTg                                                                                                                         <  2:116548‑M2/98‑1 (MQ=255)
                                           cgTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTg                                                                                                                         >  1:116548‑M2/3‑100 (MQ=255)
                                             tCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCa                                                                      >  1:1525045/1‑149 (MQ=255)
                                                                                                       aaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAAt                                                          <  2:1312452/103‑1 (MQ=255)
                                                                                                       aaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAAt                                                          >  1:1312452/1‑103 (MQ=255)
                                                                                                       aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCAcc                                                      >  4:1715686/1‑107 (MQ=255)
                                                                                                       aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCAcc                                                      <  3:1715686/107‑1 (MQ=255)
                                                                                                            tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCt       >  2:68672/1‑149 (MQ=255)
                                                                                                            tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCt       >  3:199057/1‑149 (MQ=255)
                                                                                                            tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCt       >  3:1860822/1‑149 (MQ=255)
                                                                                                            tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCt       >  2:270907/1‑149 (MQ=255)
                                                                                                             gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCtt      >  3:1446403/1‑149 (MQ=255)
                                                                                                              cTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTa     >  4:1280219/1‑149 (MQ=255)
                                                                                                               tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTaa    >  2:916291/1‑149 (MQ=255)
                                                                                                                 gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCaa            >  3:776133/1‑139 (MQ=255)
                                                                                                                 gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCaa            >  3:679816/1‑139 (MQ=255)
                                                                                                                 gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCaa            <  4:679816/139‑1 (MQ=255)
                                                                                                                 gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCaa            <  4:776133/139‑1 (MQ=255)
                                                                                                                 gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATa  >  4:1253341/1‑149 (MQ=255)
                                                                                                                 |                                                                                                                                                    
GCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATA  >  NC_000913/3708482‑3708743

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