Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,231,607 A→G 64.9% T383T (ACT→ACC ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,231,6070AG64.9% 80.6 / 65.0 111T383T (ACT→ACCptrBprotease II
Reads supporting (aligned to +/- strand):  ref base A (10/29);  new base G (46/26);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 1.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 113 total aligned reads displayed.

ACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCA  >  minE/1231539‑1231674
                                                                    |                                                                   
aCTTCCGTTTGTTTTATTACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTa                                                                     <  1:2165917/69‑1 (MQ=255)
  ttCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGt                                                                   <  1:2686963/69‑1 (MQ=255)
         ttttttAAAACACGACTCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGt                                                                   <  1:3221606/60‑1 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:2630992/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:2056091/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:287352/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:2142282/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:1729237/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:167747/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:1636839/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:3230329/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:379347/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:1361891/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:431297/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:873488/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:2787583/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:124149/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                                 >  1:2471787/1‑59 (MQ=255)
              tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCAGGTGTGGTCa                                                                 >  1:1930172/1‑59 (MQ=255)
                 tACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAc                                                     <  1:1958744/68‑1 (MQ=255)
                 tACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAc                                                     <  1:1792476/68‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAAt                                                        <  1:3006602/63‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:2832489/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:2932690/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:1797883/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:2538510/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:2822346/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:2476746/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:1420457/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:1000790/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  <  1:409629/69‑1 (MQ=255)
                   cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTCTAGTCATGGAAGAATAACCAt                                                  <  1:302048/69‑1 (MQ=255)
                        cgCTCACCTGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:245800/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGc                                             <  1:694025/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2810485/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2224888/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2166342/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2572572/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2431351/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2122444/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1807778/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2897260/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1030048/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:818912/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:799241/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:113628/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1173882/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:735819/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:723876/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1262928/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1392668/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:3073758/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1720725/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2965517/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:2897861/69‑1 (MQ=255)
                        cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                                             <  1:1940247/69‑1 (MQ=255)
                                 gTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGTATAACCATAACGc                                             <  1:656886/60‑1 (MQ=255)
                                     ccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  >  1:244202/1‑51 (MQ=255)
                                     ccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  >  1:2844628/1‑51 (MQ=255)
                                     ccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  >  1:2430752/1‑51 (MQ=255)
                                     ccatatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                                  >  1:1252544/1‑51 (MQ=255)
                                      catatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTc                               >  1:202803/1‑69 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1439544/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:3120904/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1281355/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:52244/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:610076/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1274960/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1840200/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1926184/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1239595/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1074232/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2093629/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2101085/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1040066/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2538617/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2753325/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2735063/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:1035623/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2681148/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2638253/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                                           >  1:2581698/1‑51 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCa                                            >  1:1285526/1‑50 (MQ=255)
                                            caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCa                                            >  1:1767251/1‑50 (MQ=255)
                                                 caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    >  1:1501107/1‑69 (MQ=255)
                                                 caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    <  1:2620113/69‑1 (MQ=255)
                                                 caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    >  1:301609/1‑69 (MQ=255)
                                                 caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    <  1:1808860/69‑1 (MQ=255)
                                                 caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    >  1:2915366/1‑69 (MQ=255)
                                                 caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    >  1:921425/1‑69 (MQ=255)
                                                  aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    <  1:308751/68‑1 (MQ=255)
                                                  aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                    <  1:532123/68‑1 (MQ=255)
                                                     caaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa                <  1:1933873/69‑1 (MQ=255)
                                                     caaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa                <  1:1207647/69‑1 (MQ=255)
                                                         gtgtCTGGTGTAGTCATGGAAGAATAACCATAACGc                                             <  1:2417478/36‑1 (MQ=255)
                                                            tCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATc         <  1:1565096/69‑1 (MQ=255)
                                                                   tAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTc                               <  1:2551362/40‑1 (MQ=255)
                                                                   tAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTca  <  1:430626/69‑1 (MQ=255)
                                                                   tAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTca  <  1:1681507/69‑1 (MQ=255)
                                                                   tAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTca  <  1:1990840/69‑1 (MQ=255)
                                                                    |                                                                   
ACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCA  >  minE/1231539‑1231674

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: