Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,231,598 G→A 64.7% D386D (GAC→GAT ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,231,5980GA64.7% 74.9 / 68.3 113D386D (GAC→GATptrBprotease II
Reads supporting (aligned to +/- strand):  ref base G (10/30);  new base A (47/26);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 7.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

AGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTA  >  minE/1231535‑1231660
                                                               |                                                              
aGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTgg                                                           <  1:879701/69‑1 (MQ=255)
aGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTgg                                                           <  1:2035021/69‑1 (MQ=255)
aGGAACCTCGGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTgg                                                           <  1:1884022/69‑1 (MQ=255)
    aCTTCCGTTTGTTTTATTACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTa                                                       <  1:2165917/69‑1 (MQ=255)
    aCTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTa                                                       <  1:670834/69‑1 (MQ=255)
      ttCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGt                                                     <  1:2686963/69‑1 (MQ=255)
             ttttttAAAACACGACTCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGt                                                     <  1:3221606/60‑1 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGgt                                                          >  1:2141271/1‑52 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:1366698/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:1825974/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:1729237/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:2787583/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:1701607/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:1636839/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:2056091/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:2630992/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:3230329/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:3358299/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:431297/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:124149/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:2142282/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:873488/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCa                                                   >  1:2471787/1‑59 (MQ=255)
                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCAGGTGTGGTCa                                                   >  1:1930172/1‑59 (MQ=255)
                     tACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAc                                       <  1:1792476/68‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAAt                                          <  1:3006602/63‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:1000790/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:2538510/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:2476746/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:2822346/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:2832489/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:1797883/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:2932690/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:1420457/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCAt                                    <  1:409629/69‑1 (MQ=255)
                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTCTAGTCATGGAAGAATAACCAt                                    <  1:302048/69‑1 (MQ=255)
                          gacgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTgg                                                           <  1:1553286/43‑1 (MQ=255)
                            cgCTCACCTGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:245800/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2655924/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2431351/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:735819/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:723876/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:3073758/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2897260/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2965517/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2810485/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2897861/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:1392668/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:1807778/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:113628/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:1720725/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:1173882/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:1940247/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2028650/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2224888/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2122444/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:2166342/69‑1 (MQ=255)
                            cgCTCACCGGTATCCATATCCAGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGc                               <  1:1262928/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:1252544/1‑51 (MQ=255)
                                          catatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTc                 >  1:202803/1‑69 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1239595/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1074232/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1274960/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1439544/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:52244/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:3120904/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:610076/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1281355/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1840200/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2844129/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1926184/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1040066/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2753325/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2681148/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2093629/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2638253/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2101085/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:1035623/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2581698/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2538617/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2479384/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCaa                             >  1:2358148/1‑51 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCa                              >  1:1767251/1‑50 (MQ=255)
                                                caGTTCAAACAAAGTATCTGGTGTGGTCATGGACGAATAACCATAACGCa                              >  1:1285526/1‑50 (MQ=255)
                                                    tCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTgg        >  1:2989885/1‑68 (MQ=255)
                                                     caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt      >  1:1501107/1‑69 (MQ=255)
                                                     caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt      >  1:2915366/1‑69 (MQ=255)
                                                     caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt      <  1:2620113/69‑1 (MQ=255)
                                                     caaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt      >  1:921425/1‑69 (MQ=255)
                                                      aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt      <  1:532123/68‑1 (MQ=255)
                                                      aaacaaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt      <  1:308751/68‑1 (MQ=255)
                                                         caaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa  <  1:1207647/69‑1 (MQ=255)
                                                         caaaGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa  <  1:1933873/69‑1 (MQ=255)
                                                             gtgtCTGGTGTAGTCATGGAAGAATAACCATAACGc                               <  1:2417478/36‑1 (MQ=255)
                                                               |                                                              
AGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTA  >  minE/1231535‑1231660

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: