Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,730,046 C→T 36.0% T177T (ACG→ACA tadA ← tRNA‑specific adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,730,0460CT36.0% 44.7 / 38.6 70T177T (ACG→ACAtadAtRNA‑specific adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (11/33);  new base T (16/9);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 2.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGT  >  minE/1729985‑1730104
                                                             |                                                          
tGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAg                                                      <  1:900039/68‑1 (MQ=255)
tGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAg                                                      <  1:841922/68‑1 (MQ=255)
tGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAg                                                      <  1:537641/68‑1 (MQ=255)
tGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAg                                                      <  1:1678125/68‑1 (MQ=255)
tGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAg                                                      <  1:2729991/68‑1 (MQ=255)
                tctGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTg                                     >  1:3336022/1‑69 (MQ=255)
                tctGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTg                                     >  1:1811530/1‑69 (MQ=255)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:3153489/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:211994/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:2253356/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:1930744/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:670330/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:699510/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:857474/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:980540/45‑1 (MQ=38)
                    ctgTTTTCCAGAAAAGGGAGTGAAGAGAAACAAAATTAATCTGTc                                                         <  1:1051654/45‑1 (MQ=38)
                      gTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCttt                               <  1:2679933/69‑1 (MQ=255)
                         ttCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCttt                               <  1:261681/66‑1 (MQ=255)
                         ttCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCAGTCGAGGATTGCGCTTTTTTCTGCGCTTTAAt                            >  1:2436451/1‑69 (MQ=255)
                           ccAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAAtt                           <  1:3232396/68‑1 (MQ=255)
                           ccAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAAtt                           <  1:2832756/68‑1 (MQ=255)
                           ccAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAAtt                           <  1:1025617/68‑1 (MQ=255)
                           ccAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAAtt                           <  1:3342482/68‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:39728/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:2414909/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:462008/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:2455229/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:506663/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:606085/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:1647907/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:2737661/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                <  1:70005/69‑1 (MQ=255)
                                     gAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                >  1:320792/1‑69 (MQ=255)
                                      aGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAAtt                           <  1:456379/57‑1 (MQ=255)
                                       gTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                >  1:3302745/1‑67 (MQ=255)
                                       gTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                >  1:650743/1‑67 (MQ=255)
                                       gTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                >  1:2884853/1‑67 (MQ=255)
                                       gTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                >  1:1264884/1‑67 (MQ=255)
                                            gagaAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCg          <  1:1051742/68‑1 (MQ=255)
                                            gagaAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCg          <  1:2651953/68‑1 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:514009/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:998077/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:951599/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:1224532/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:722956/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:12949/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:1692346/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:524817/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:3074446/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:1777042/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:471490/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:1955336/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:1015563/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:3402374/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGc           >  1:2274894/1‑66 (MQ=255)
                                             agaAACAAAATTAATCTGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGa         >  1:1938/1‑68 (MQ=255)
                                              gaAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGaa        >  1:2006972/1‑68 (MQ=255)
                                              gaAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGaa        >  1:534882/1‑68 (MQ=255)
                                              gaAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGaa        >  1:2083496/1‑68 (MQ=255)
                                              gaAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGa         >  1:3397006/1‑67 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAATGAAGt  <  1:1392962/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGGAGt  <  1:1731192/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:1501227/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:322766/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:1090311/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:2636531/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:1643687/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:613094/69‑1 (MQ=255)
                                                   aaaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:3415783/69‑1 (MQ=255)
                                                    aaaTTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGt  <  1:3192883/68‑1 (MQ=255)
                                                             |                                                          
TGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGT  >  minE/1729985‑1730104

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