Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,714,472 G→C 32.0% A320A (GCC→GCG yphG ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,714,4720GC32.0% 114.8 / 73.4 103A320A (GCC→GCGyphGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base C (22/11);  total (57/46)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 103 total aligned reads displayed.

TACTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCA  >  minE/1714413‑1714521
                                                           |                                                 
tACTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGTTTATCGGCaaaaa                                            <  1:2214775/67‑1 (MQ=255)
tACTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCaaaaa                                            <  1:791374/67‑1 (MQ=255)
tACTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCaaaaa                                            <  1:1359580/67‑1 (MQ=255)
 aCTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATAc                                         >  1:958396/1‑69 (MQ=255)
 aCTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATAc                                         >  1:89668/1‑69 (MQ=255)
 aCTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATAc                                         >  1:427593/1‑69 (MQ=255)
 aCTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATAc                                         >  1:114479/1‑69 (MQ=255)
 aCTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATAc                                         >  1:598713/1‑69 (MQ=255)
 aCTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGAGCTGGTTATCGGCAAAAATAc                                         >  1:1156137/1‑69 (MQ=255)
       cAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGt                                     <  1:2470035/67‑1 (MQ=255)
       cAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGt                                     <  1:1576970/67‑1 (MQ=255)
       cAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGt                                     <  1:1236682/67‑1 (MQ=255)
           aCGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAg                               >  1:1456754/1‑69 (MQ=255)
           aCGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAg                               >  1:345126/1‑69 (MQ=255)
           aCGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAg                               >  1:664840/1‑69 (MQ=255)
           aCGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAg                               >  1:969390/1‑69 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:68990/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:106402/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:2508572/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:1000291/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:3230018/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:3346326/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:904669/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:653889/69‑1 (MQ=255)
            cGCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:1389009/69‑1 (MQ=255)
             gCTTCTCTTCATAAGCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:1628954/68‑1 (MQ=255)
                tctcTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGa                          <  1:30150/69‑1 (MQ=255)
                tctcTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGa                          <  1:2131507/69‑1 (MQ=255)
                tctcTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGa                          <  1:2359782/69‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:813179/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:78212/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:2948397/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:293229/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:2829794/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:694656/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:2655909/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:2339830/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:1948379/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:1110059/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:2145972/68‑1 (MQ=255)
                   ctTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:2085733/68‑1 (MQ=255)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1382836/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:2986087/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:2999099/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1837242/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:3106790/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1756369/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1745223/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1710080/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:538607/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1568447/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:2306156/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1226374/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:702331/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1208382/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:875375/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:1097920/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:93545/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:969029/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:2895806/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:2012901/1‑39 (MQ=38)
                           gCATCAAGCCAGGTAAAATCGGGCTGATTATCCGCaaaa                                             >  1:2774205/1‑39 (MQ=38)
                               cAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTAt           >  1:2495712/1‑69 (MQ=255)
                               cAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTAt           >  1:35712/1‑69 (MQ=255)
                                     aGGTAAAATCGGGCTGATTATCCGCAAAAATGCCGGTCATCAGt                              <  1:644028/44‑1 (MQ=38)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTc   >  1:2778877/1‑68 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:1044904/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:2183025/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:1629877/1‑69 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:371300/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:347830/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:2696645/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:3378969/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:2347840/69‑1 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:3350809/1‑69 (MQ=255)
                                        tAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:2195878/1‑69 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:2320531/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:82969/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:104687/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:1143350/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:2439403/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:738495/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:2289827/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:2057332/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:2526185/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:1398686/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:14500/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:455092/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcatc                                 >  1:3341423/1‑37 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGtcat                                  >  1:3155131/1‑36 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:2477495/68‑1 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:2560445/1‑68 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:565103/1‑68 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:2969252/1‑68 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:3376224/68‑1 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  <  1:1074049/68‑1 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:2939356/1‑68 (MQ=255)
                                         aaaaTCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCa  >  1:2973069/1‑68 (MQ=255)
                                          aaaTCGGTCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGata                        <  1:527497/45‑1 (MQ=255)
                                                 gCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGa                          <  1:1871640/36‑1 (MQ=255)
                                                           |                                                 
TACTGCTCAAAACGCTTCTCTTCGTAAGCATCAAGCCAGGTAAAATCGGGCTGGTTATCGGCAAAAATACCGGTCATCAGTTCGATATACGGGCCGTTATTGTCGGTCA  >  minE/1714413‑1714521

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