Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,335,681 T→C 8.0% intergenic (+223/‑150) yciX → / → acnA uncharacterized protein/aconitate hydratase 1; aconitase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,6810TC 8.0% 164.4 / 4.4 50intergenic (+223/‑150)yciX/acnAuncharacterized protein/aconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (2/2);  total (25/25)
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.17e-01

CCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTAT  >  NC_000913/1335554‑1335771
                                                                                                                               |                                                                                          
ccGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGcggcgg                                                                                     <  4:468850/135‑1 (MQ=255)
ccGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGcggcgg                                                                                     >  3:468850/1‑135 (MQ=255)
                           acaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCcaca                                                               >  3:135463/1‑130 (MQ=255)
                           acaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCcaca                                                               <  4:135463/130‑1 (MQ=255)
                            caAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAgg                                                            <  4:1044127/132‑1 (MQ=255)
                            caAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAgg                                                            >  3:1044127/1‑132 (MQ=255)
                              aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGcggcgg                                                                                     <  3:841215/105‑1 (MQ=255)
                              aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGcggcgg                                                                                     >  4:841215/1‑105 (MQ=255)
                                                        cgcgTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAtt                                        >  4:424416/1‑124 (MQ=255)
                                                        cgcgTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAtt                                        <  3:424416/124‑1 (MQ=255)
                                                             tttCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                      >  4:878593/1‑73 (MQ=255)
                                                             tttCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                      <  3:878593/73‑1 (MQ=255)
                                                                ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACaa                                                              <  2:421798/94‑1 (MQ=255)
                                                                ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACaa                                                              >  1:421798/1‑94 (MQ=255)
                                                                    gTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTg                                      <  2:961169/114‑1 (MQ=255)
                                                                    gTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTg                                      >  1:961169/1‑114 (MQ=255)
                                                                     tCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctttg                                                                    >  4:278961/1‑83 (MQ=255)
                                                                     tCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctttg                                                                    <  3:278961/83‑1 (MQ=255)
                                                                        aTGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCaa                          >  3:907711/1‑122 (MQ=255)
                                                                        aTGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCaa                          <  4:907711/122‑1 (MQ=255)
                                                                         tGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgct                                                                       >  2:864620/1‑76 (MQ=255)
                                                                         tGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgct                                                                       <  1:864620/76‑1 (MQ=255)
                                                                         tGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGcagcgg                                                                                     <  3:444845/62‑1 (MQ=255)
                                                                         tGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGcagcgg                                                                                     >  4:444845/1‑62 (MQ=255)
                                                                          gCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCttt                                                <  3:414538/98‑1 (MQ=255)
                                                                          gCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCttt                                                >  4:414538/1‑98 (MQ=255)
                                                                             cttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgt       >  3:996911/1‑136 (MQ=255)
                                                                             cttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgt       <  4:996911/136‑1 (MQ=255)
                                                                               tctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGcagcgg                                                                                     >  3:329757/1‑56 (MQ=255)
                                                                               tctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGcagcgg                                                                                     <  4:329757/56‑1 (MQ=255)
                                                                                               ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGTTTGTTATCaa                          <  3:652187/99‑1 (MQ=255)
                                                                                               ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCaa                          >  3:197202/1‑99 (MQ=255)
                                                                                               ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCaa                          <  4:197202/99‑1 (MQ=255)
                                                                                               ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCaa                          >  4:652187/1‑99 (MQ=255)
                                                                                                       ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTctcctc                                                   >  3:803554/1‑66 (MQ=255)
                                                                                                       ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTctcctc                                                   <  4:803554/66‑1 (MQ=255)
                                                                                                        tACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgttgt                                <  3:22156/84‑1 (MQ=255)
                                                                                                        tACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgttgt                                >  4:22156/1‑84 (MQ=255)
                                                                                                            cACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGtt                    <  2:829511/92‑1 (MQ=255)
                                                                                                            cACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGtt                    >  1:829511/1‑92 (MQ=255)
                                                                                                                gATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt     <  2:652672/103‑1 (MQ=255)
                                                                                                                gATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt     >  1:652672/1‑103 (MQ=255)
                                                                                                                  tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAtt                                        <  2:529092/66‑1 (MQ=255)
                                                                                                                  tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAtt                                        >  1:529092/1‑66 (MQ=255)
                                                                                                                   ccTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGtt          >  2:26233/1‑95 (MQ=255)
                                                                                                                   ccTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGtt          <  1:26233/95‑1 (MQ=255)
                                                                                                                      tGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAAt                        <  4:549016/78‑1 (MQ=255)
                                                                                                                      tGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAAt                        >  3:549016/1‑78 (MQ=255)
                                                                                                                           ccGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTAt  >  2:173240/1‑95 (MQ=255)
                                                                                                                           ccGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTAt  <  1:173240/95‑1 (MQ=255)
                                                                                                                               |                                                                                          
CCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTAT  >  NC_000913/1335554‑1335771

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