Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 67,890 +TG coding (448/1503 nt) araA ← L‑arabinose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091367,8871.G92.9% 68.0 / 0.0 28coding (451/1503 nt)araAL‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (14/12);  total (15/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 = 8.72e-01
*NC_00091367,8872.T92.6% 73.0 / ‑0.2 28coding (451/1503 nt)araAL‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (14/11);  total (15/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.48e-01

TTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGA  >  NC_000913/67753‑68018
                                                                                                                                       ||                                                                                                                                     
tttATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑gtgacc                                                                                                                             >  1:648278/1‑138 (MQ=255)
         aTCGGTGACCGCCACTTCACGCATTTTATCTCCAAATCGACAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCAc                                                                                                                      <  1:393801/143‑1 (MQ=255)
              tGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCCCGGCAt                                                                                                                 >  1:472690/1‑143 (MQ=255)
                 tctctaCTTCACGCTTGTTATCGTCATATCGGTAGACTTTCATATGACGGGTATTCTGTTTAGTGTCCGCTTGATGTATCTAGTAGCCGATACGCTCATGGGCTTGTTTATCCTGTCAGTGT‑‑GTCCGGTAACCACGGCATGtt                                                                                                              <  2:605099/138‑1 (MQ=255)
                    ccACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACGACGAAATGATGCt                                                                                                           >  2:25/1‑143 (MQ=255)
                     cACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGCTGCCGGGTATCCTGTATAGCGACCGCCTGACGCATCCAGGGGCCGATACGCTCATGGGCTTGTTTAACCTGCCAGTGT‑‑GACCGGTAACCACGCCATGctgctg                                                                                                          >  2:319381/1‑143 (MQ=255)
                      gctgcACGCATGGTATCGCCAAATCGGCAGACTTGCAGGTGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGa                                                                                                         <  2:611711/139‑1 (MQ=255)
                         tCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCTATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAAGCACGGCATGTTGCTGACGc                                                                                                      >  1:361284/1‑143 (MQ=255)
                                 gTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATAcgcgc                                                                                              <  1:606392/143‑1 (MQ=255)
                                                gCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGaa                                                                               >  2:101658/1‑143 (MQ=255)
                                                                   gTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCCCCAttttca                                                            >  1:718083/1‑138 (MQ=255)
                                                                      tccgttttAGGGACCGCCGGACTCATCCGGTACCCGTTACGCTCATGGGCTTGTTTATCCTGCCAGGGT‑‑GACCGGTAACCACGGCATGTTGCTGCCGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTc                                                         <  2:76710/138‑1 (MQ=255)
                                                                         tGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCCGCCTgg                                                      >  2:305453/1‑143 (MQ=255)
                                                                           gttggacaCCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggtt                                                    <  2:427988/136‑1 (MQ=255)
                                                                              agagACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCCATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcac                                                 >  1:447444/1‑142 (MQ=255)
                                                                                cgcccGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggt                                               <  2:611397/140‑1 (MQ=255)
                                                                                agaCCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGCCCGCCGTGTGCAGTCTggttaagcg                                               >  1:125869/1‑141 (MQ=255)
                                                                                     gCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATa                                          <  2:19388/143‑1 (MQ=255)
                                                                                          aCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTc                                     <  1:25/143‑1 (MQ=255)
                                                                                            gCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATACTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCCATGAAGCCGAACTCGCGACCGCCATGTGCATCCTGGTTTAAGTTCTTAAAGTCCa                                   >  1:464828/1‑143 (MQ=255)
                                                                                                    gAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAg                                       >  2:232093/1‑131 (MQ=255)
                                                                                                    gAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAg                                       <  1:232093/131‑1 (MQ=255)
                                                                                                         gATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCCATGAAGCCGAACTCGCGCCCGCCATGTGCCGTCTGGTTCTGCTCCTTATCGTCTATATAGATaatata                      >  1:311862/1‑138 (MQ=255)
                                                                                                            atgCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCCCGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCa                   <  1:320620/141‑1 (MQ=255)
                                                                                                              gCTCATGGGCTTGTTTATCCTGCCA‑‑GTGGGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACg                 <  1:614640/143‑1 (MQ=255)
                                                                                                                     ggCCTGTTCATCCTGCGAGTGT‑‑GGCCGGTAACCGCGGCATGTTGCTGCCGCATGCGCGCGCCAATGAAGCCGCACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAgcgc          <  2:658222/143‑1 (MQ=255)
                                                                                                                          gTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACGGTCCCACGGCAGCGCCCCGt     >  1:380385/1‑143 (MQ=255)
                                                                                                                             tATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGa  >  1:27684/1‑143 (MQ=255)
                                                                                                                                       ||                                                                                                                                     
TTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGA  >  NC_000913/67753‑68018

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