Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,952,505 (G)8→9 noncoding (72/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,952,4971.G81.2% 36.0 / 4.9 16noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (6/7);  total (7/9)
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.91e-01

GCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACC  >  minE/2952428‑2952559
                                                                       |                                                              
gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               <  1:3349012/71‑1 (MQ=255)
gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               <  1:1726933/71‑1 (MQ=255)
gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               <  1:2972384/71‑1 (MQ=255)
 cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               <  1:3918885/70‑1 (MQ=255)
       tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                 >  1:94068/1‑63 (MQ=255)
       tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                 >  1:1773986/1‑63 (MQ=255)
       tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                 >  1:4238408/1‑63 (MQ=255)
       tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                 >  1:2180775/1‑63 (MQ=255)
         ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                         <  1:1230459/68‑1 (MQ=255)
         ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                         <  1:247419/68‑1 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:879678/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:392835/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:3963061/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:3026063/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:2716454/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:72576/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:2056311/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:1939735/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:883631/1‑57 (MQ=255)
              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                               >  1:1029015/1‑57 (MQ=255)
               ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGaa                                                 >  1:1362171/1‑71 (MQ=255)
                       ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                          <  1:2395265/70‑1 (MQ=255)
                                  aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                          <  1:3559298/59‑1 (MQ=37)
                                  aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                 >  1:1093642/1‑68 (MQ=255)
                                  aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                               >  1:3613259/1‑70 (MQ=255)
                                  aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                               >  1:271549/1‑70 (MQ=255)
                                  aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGa                                >  1:3094204/1‑69 (MQ=255)
                                  aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg                                                        >  1:1185032/1‑44 (MQ=255)
                                    cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                          <  1:3729535/57‑1 (MQ=25)
                                     ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                          <  1:4448720/56‑1 (MQ=25)
                                                      gCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg          >  1:111233/1‑71 (MQ=255)
                                                      gCCGCGTTTGGTGCGTGGGGGGGGGACTTGAACCCCCACGTCCTTAAGaa                               <  1:628071/50‑1 (MQ=14)
                                                           gTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                 <  1:35436/43‑1 (MQ=11)
                                                           gTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                               <  1:426931/45‑1 (MQ=11)
                                                           gTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAac                        <  1:2410704/52‑1 (MQ=14)
                                                             ttGGTGCGTGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAac                        <  1:3062683/50‑1 (MQ=11)
                                                              tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc  >  1:4322517/1‑70 (MQ=35)
                                                                       |                                                              
GCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACC  >  minE/2952428‑2952559

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