Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 249,414 A→C 12.6% intergenic (+35/‑41) brnQ → / → proY predicted branched chain amino acid transporter/predicted cryptic proline transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE249,4140AC12.6% 66.5 / 3.1 31intergenic (+35/‑41)brnQ/proYpredicted branched chain amino acid transporter/predicted cryptic proline transporter
Reads supporting (aligned to +/- strand):  ref base A (10/17);  new base C (2/2);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGG  >  minE/249348‑249484
                                                                  |                                                                      
cAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGcc                                                                    <  1:3687704/71‑1 (MQ=255)
      gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTTTTTTAACCACGGGGCATCGCAGCCCCGTg                                                               <  1:4007120/70‑1 (MQ=25)
      gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                               <  1:1648202/70‑1 (MQ=37)
                    aGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:1055242/71‑1 (MQ=255)
                    aGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:1235186/71‑1 (MQ=255)
                    aGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:1421521/71‑1 (MQ=255)
                    aGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:480523/71‑1 (MQ=255)
                     gcgcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:102235/70‑1 (MQ=255)
                        cTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:1556641/67‑1 (MQ=255)
                           aCTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGt                                         >  1:1155508/1‑71 (MQ=255)
                            cTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:800536/63‑1 (MQ=255)
                            cTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGtt                                                <  1:2222894/63‑1 (MQ=255)
                               aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGtttttt                                                        >  1:3737439/1‑52 (MQ=255)
                               aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAgg                                     <  1:3641500/71‑1 (MQ=255)
                               aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAg                                      >  1:2257211/1‑70 (MQ=255)
                               aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAg                                      >  1:3188291/1‑70 (MQ=255)
                               aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAg                                      >  1:1704588/1‑70 (MQ=255)
                               aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGAAGGGTTAg                                      >  1:1709453/1‑70 (MQ=255)
                                  cACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAt                                  >  1:4063729/1‑71 (MQ=255)
                                    cTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTg                                <  1:758387/71‑1 (MQ=255)
                                       aaCATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                               <  1:806386/69‑1 (MQ=255)
                                        aCATTTGTTTTAACCACGGGGCTGCGCTGCCCCGTGGTTTTTTATtgtg                                                  >  1:596464/1‑49 (MQ=255)
                                        aCATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGttt                                                           >  1:3677541/1‑40 (MQ=255)
                                           tttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaaa                         <  1:3603529/71‑1 (MQ=255)
                                            ttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAg                        >  1:3507603/1‑71 (MQ=255)
                                            ttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAg                        >  1:1826248/1‑71 (MQ=255)
                                             tGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAg                        <  1:562965/70‑1 (MQ=255)
                                              gTTTTAACCACGGGGCTGCGCTGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGaa                                   >  1:3526216/1‑58 (MQ=255)
                                                       aCGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGa                               <  1:1651436/53‑1 (MQ=255)
                                                          gggCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAtt                                 <  1:803254/48‑1 (MQ=255)
                                                             cTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGc       <  1:4249942/71‑1 (MQ=255)
                                                                  aTGCCCCGTGGTTTTTTATTGTGTTGATGGTTTAgg                                     <  1:2446801/36‑1 (MQ=255)
                                                                  aTGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAgg                                     <  1:2851897/36‑1 (MQ=255)
                                                                  aTGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAgg                                     <  1:1933706/36‑1 (MQ=255)
                                                                  aTGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAg              <  1:607190/59‑1 (MQ=255)
                                                                  aTGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTggg  >  1:2024805/1‑71 (MQ=255)
                                                                  |                                                                      
CAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGG  >  minE/249348‑249484

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: