Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,590 G→T 12.9% intergenic (‑61/+32) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5900GT12.9% 86.5 / 3.5 38intergenic (‑61/+32)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base T (2/3);  total (18/20)
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

CACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGA  >  NC_000913/83455‑83726
                                                                                                                                       |                                                                                                                                        
cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCggg                                                                                                                                   <  2:157952/143‑1 (MQ=255)
    cTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCGTTTAGAGCGCCTCGCTTCGGGCATCAAAACACCCGCGCAATGGCGCGGGtttt                                                                                                                               >  1:211413/1‑143 (MQ=255)
     tccggCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTAGCGAGAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCGTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttttt                                                                                                                              <  2:603772/139‑1 (MQ=255)
     tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCGATGGCGCGGGttttt                                                                                                                              >  2:273865/1‑143 (MQ=255)
     tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttttt                                                                                                                              >  2:266441/1‑143 (MQ=255)
      cAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGtttttt                                                                                                                             >  1:643246/1‑143 (MQ=255)
      cAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGtttttt                                                                                                                             <  1:612213/143‑1 (MQ=255)
        ggcgggCCTGGAACGCGTGTGCACCGTCGCGCAGTGTGGTATCGAAAATAATGACTTGGTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgt                                                                                                                           <  1:269505/138‑1 (MQ=255)
          cTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgttt                                                                                                                         >  1:461247/1‑143 (MQ=255)
             gCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGAc                                                                                                                      >  1:269708/1‑143 (MQ=255)
                            tCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTaa                                                                                                       <  1:557655/143‑1 (MQ=255)
                                     cgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATAGTTTGGGTCCTGGTCTCTTTTAGAGCGCCTAGCTTCGGGGATAAAAAAACCCGCGCAATGGTGCGGGGTGTTTGTGTGGCTGCGTGGTGGCTGAATGATggtgg                                                                                              >  1:275512/1‑140 (MQ=255)
                                                 tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGt                                                                                  <  1:266441/143‑1 (MQ=255)
                                                 tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCGTAATGCTGGATGCCGCTCACTCGt                                                                                  >  1:309374/1‑143 (MQ=255)
                                                              cTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaag                                                                     >  2:221479/1‑143 (MQ=255)
                                                               ttGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATACCGCTCACTCGTCTACCGCGCAaaga                                                                    <  1:511544/143‑1 (MQ=255)
                                                               ttGCTGGCTCATGGTTTGGGCCCTTGTCTCTTTTAGAGCGCCTCGCTTCGCGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaaga                                                                    <  1:511581/143‑1 (MQ=255)
                                                                tGCTGGCGCATGGTGTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                                                                   <  1:267777/143‑1 (MQ=255)
                                                                   tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGGTGCCGCTCACTCGTCTACCGCGCAACGACGAt                                                                >  2:448694/1‑143 (MQ=255)
                                                                     ggtCATGGTTTGGGTCCTTGTCCCTTTCAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                              <  1:274921/141‑1 (MQ=255)
                                                                        catgcgttGGGGCCTTGGCTCTTTGAGAGCGCCGCGCTTGTGGCATAGAAAGACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                           <  2:158101/137‑1 (MQ=255)
                                                                           ggTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAGAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                                                        <  2:4350/143‑1 (MQ=255)
                                                                           ggTTTGGGTCCTTGGCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAGGCCCGCGCGATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATTCTGGATGCCGCTCACTCGTCGACCGCGCAAAGAAGATGCGTTtag                                                        <  1:81405/143‑1 (MQ=255)
                                                                                   tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagta                                                >  1:439394/1‑143 (MQ=255)
                                                                                      ttGTCTCTTTTAGAGCGCCTAGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAGAGAAGATGCGTTTAGTAGTAGTAGAc                                             >  2:160741/1‑143 (MQ=255)
                                                                                       ggtggctTTTAGGGCGCCGCGCTTCGGGCATAAAAAGAGCCGCGCAATTGCGCGGGTTTTTTGTTTGCCTCCCTTCTTGCTTACTGCTGGGTGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAcc                                            <  2:232550/138‑1 (MQ=255)
                                                                                         tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGCCCGa                                          >  1:39177/1‑143 (MQ=255)
                                                                                           tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATa                                        <  2:82068/143‑1 (MQ=255)
                                                                                           tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATCCCGCTCCCTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGCTa                                        >  1:24628/1‑143 (MQ=255)
                                                                                             gtttAGAGCGCCTCGCGTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaaa                                      <  1:40630/142‑1 (MQ=255)
                                                                                                      gcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                             >  1:576834/1‑143 (MQ=255)
                                                                                                                    aTAAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt               <  1:400394/143‑1 (MQ=255)
                                                                                                                     tAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGGTGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATc              >  2:550456/1‑143 (MQ=255)
                                                                                                                       aaaaaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCACTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg            >  1:119549/1‑143 (MQ=255)
                                                                                                                         aagaCCCGCGCAAGGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCt          <  1:159835/140‑1 (MQ=255)
                                                                                                                            aCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa       >  2:178477/1‑143 (MQ=255)
                                                                                                                               cgcgCAATTGCGCGGGTTTTTCGTTTGACTGCGTGCGGGCTTAATACTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGAt    <  1:211073/143‑1 (MQ=255)
                                                                                                                                 cgcAATCGCGCGGGTCTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGa  <  1:171930/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGA  >  NC_000913/83455‑83726

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