Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 86,384 C→A intergenic (‑51/+42) 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
*REL60686,3840CA91.3% 68.9 / NA 23intergenic (‑51/+42)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (12/9);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACC  >  REL606/86311‑86455
                                                                         |                                                                       
aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAacc                                                                       >  2:84959/1‑76 (MQ=255)
 aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACcg                                                                      >  2:725475/1‑76 (MQ=255)
        gANTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACTGCGCANTg                                                               <  2:1095901/76‑1 (MQ=255)
        gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATg                                                               >  1:805419/1‑76 (MQ=255)
            tGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGgcgc                                                           <  2:654531/76‑1 (MQ=255)
                 gCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGtt                                                      <  2:1048835/76‑1 (MQ=255)
                  cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGttt                                                     <  2:1255923/76‑1 (MQ=255)
                  cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGttt                                                     >  2:413898/1‑76 (MQ=255)
                  cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGttt                                                     <  2:261139/76‑1 (MQ=255)
                   tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGgttt                                                    >  1:893685/1‑76 (MQ=255)
                        gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGGTTtttgtt                                               >  1:1055589/1‑76 (MQ=255)
                           tGGGTCCTTCTCTTTTTTAGAGCCCCTCGCTGCGGGCAAAAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGa                                            <  1:867197/76‑1 (MQ=37)
                             ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACt                                          >  2:563609/1‑76 (MQ=255)
                              gTCCTTGTCTCTTTAAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGACTg                                         <  2:748751/76‑1 (MQ=255)
                               tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGc                                        >  2:1158807/1‑76 (MQ=255)
                                ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAACCCGCCAAGGGGCGGGGTTTTTTTTTTTCTGCg                                       >  1:1215499/1‑76 (MQ=255)
                                   tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAACACCGCGCCCTTGCGCGGGTGGTTTGTTTGTCTGCGTGc                                    >  1:753191/1‑76 (MQ=255)
                                    gTCTCTTTTATAGCGCCCCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCt                                   >  1:113511/1‑76 (MQ=39)
                                       tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGc                                >  2:557548/1‑76 (MQ=255)
                                                gcgcCTCGCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTGTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                       <  1:1270169/76‑1 (MQ=255)
                                                   ccTCGCTTCGTGCATAAAAAAAACCGCGCAATGGCGCGGGGTTTTTGTTTCACTACGTGGTGGCTTTATGCTGGAt                    >  2:598694/1‑76 (MQ=255)
                                                       gCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg                >  2:648321/1‑76 (MQ=255)
                                                       gCTTCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTGTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg                <  2:258325/76‑1 (MQ=255)
                                                        cTTCGGGCATAAGAACAACAGCGCAATGGTGCGTGTTTTTTGTATGAGTGCGTGCTGGCTTAATGCTGGATGCCGc               >  1:1296211/1‑76 (MQ=255)
                                                        cTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGc               <  1:245700/76‑1 (MQ=255)
                                                         ttCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTGTTTTGTTTGACTGTGTGCTGGCTTAATGCTGGATGCCGCt              <  2:1263564/76‑1 (MQ=255)
                                                          tCGGGCATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCNCTc             <  1:1201941/76‑1 (MQ=255)
                                                               cATAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCg        >  2:451918/1‑76 (MQ=255)
                                                                aTAAAAAAAACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGt       >  2:1209941/1‑76 (MQ=255)
                                                                  aaaaaaaaCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCt     >  2:143977/1‑76 (MQ=255)
                                                                   aaaaagACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTa    <  2:49454/76‑1 (MQ=255)
                                                                    aaagaCCCGCGCAGTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAc   <  1:437915/76‑1 (MQ=255)
                                                                     aaaaACCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGTTGGCTTAATGCTGGATGCCGCTCACTCGTCTAcc  >  2:1279381/1‑76 (MQ=255)
                                                                         |                                                                       
AAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACC  >  REL606/86311‑86455

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