Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,383 A→G 22.2% intergenic (+189/+40) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3830AG22.2% 54.0 / 11.6 27intergenic (+189/+40)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base G (3/3);  total (14/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 = 9.79e-01

AGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTG  >  NC_000913/3423264‑3423413
                                                                                                                       |                              
agtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt    >  3:47975‑M1/113‑148 (MQ=255)
              cccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacttcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc       >  4:1038098‑M1/99‑131 (MQ=255)
              cccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaCTAAAACAAAAGGCCCAGTCTTCCGACTGAGcc       <  3:1038098‑M1/33‑1 (MQ=255)
                      cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                    >  3:661421/1‑110 (MQ=40)
                      cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCAGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                    <  4:661421/110‑1 (MQ=255)
                          catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGa          >  1:491071‑M1/87‑116 (MQ=255)
                          catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGa          <  2:491071‑M1/30‑1 (MQ=255)
                             cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg         <  4:71976/114‑1 (MQ=40)
                             cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg         >  3:71976/1‑114 (MQ=40)
                                 gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              <  4:713140/105‑1 (MQ=255)
                                 gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              >  3:713140/1‑105 (MQ=40)
                                       gaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTg           <  4:881736‑M1/29‑1 (MQ=255)
                                       gaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTg           >  3:881736‑M1/74‑102 (MQ=255)
                                                cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt      <  2:358880‑M1/34‑1 (MQ=255)
                                                cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt      >  1:358880‑M1/65‑98 (MQ=255)
                                                     tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                        <  3:629843/75‑1 (MQ=255)
                                                     tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                        >  4:629843/1‑75 (MQ=255)
                                                        aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      <  2:791787/90‑1 (MQ=255)
                                                        aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      >  1:691611/1‑90 (MQ=255)
                                                        aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      >  1:791787/1‑90 (MQ=255)
                                                        aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      <  2:691611/90‑1 (MQ=255)
                                                                cccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg  >  4:149945/1‑86 (MQ=255)
                                                                cccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg  <  3:149945/86‑1 (MQ=255)
                                                                              cctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc       >  3:926911‑M1/35‑67 (MQ=255)
                                                                              cctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc       <  4:926911‑M1/33‑1 (MQ=255)
                                                                                      ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              <  4:474855/52‑1 (MQ=255)
                                                                                      ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              >  3:474855/1‑52 (MQ=255)
                                                                                                                       |                              
AGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTG  >  NC_000913/3423264‑3423413

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