Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,047,174 G→A 24.1% I40I (ATC→ATT ydhH ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,047,1740GA24.1% 113.3 / 34.9 58I40I (ATC→ATTydhHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/16);  new base A (8/6);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACA  >  minE/1047107‑1047243
                                                                   |                                                                     
cAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATc                                                                    >  1:1780844/1‑71 (MQ=255)
   aCTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGc                                                                 >  1:1326641/1‑71 (MQ=255)
    cTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCgg                                                                  >  1:276077/1‑69 (MQ=255)
    cTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGc                                                                 >  1:2007939/1‑70 (MQ=255)
    cTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGcc                                                                <  1:1214565/71‑1 (MQ=255)
          aaGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCa                                                          <  1:2842943/71‑1 (MQ=255)
          aaGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCa                                                          <  1:3460200/71‑1 (MQ=255)
             tgtAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCggg                                                                       <  1:3907244/55‑1 (MQ=255)
             tgtAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAAc                                                       <  1:2813595/71‑1 (MQ=255)
             tgtAAGCTGCTGGCCCTGACAAATATCCAGTACAGCCTGTTTCAGAGATACCggg                                                                       <  1:1530173/55‑1 (MQ=255)
              gtAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCgg                                                                  <  1:1982997/59‑1 (MQ=255)
              gtAAGCTGCTGGCCCTGGCAAATATCCAGTACAGACTGTTTCAGAGATACCGGGATCgg                                                                  <  1:2511669/59‑1 (MQ=255)
               tAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACtt                                                     <  1:794279/71‑1 (MQ=255)
                aaGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGcc                                                                >  1:341170/1‑59 (MQ=255)
                 aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACt                                                      <  1:2711617/68‑1 (MQ=255)
                 aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGc                                                   >  1:3291858/1‑71 (MQ=255)
                 aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGc                                                   >  1:27738/1‑71 (MQ=255)
                 aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGc                                                   >  1:2288798/1‑71 (MQ=255)
                 aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGc                                                   >  1:818605/1‑71 (MQ=255)
                  gctgctGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGc                                                   >  1:2464201/1‑70 (MQ=255)
                  gctgctGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGc                                                   >  1:1881296/1‑70 (MQ=255)
                      ctgGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCaa                                                              >  1:2268431/1‑55 (MQ=255)
                      ctgGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAg                                                <  1:158209/69‑1 (MQ=255)
                      ctgGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAg                                                <  1:1203/69‑1 (MQ=255)
                       tgGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCaaa                                                        >  1:2471228/1‑60 (MQ=255)
                             tGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCggg                                                                       <  1:2826918/39‑1 (MQ=255)
                             tGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCggg                                                                       <  1:3868437/39‑1 (MQ=255)
                               gCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGc                                                   >  1:892638/1‑57 (MQ=255)
                               gCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGc                                                   >  1:2856720/1‑57 (MQ=255)
                                 aaaTATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCaaa                                                        >  1:3053385/1‑50 (MQ=255)
                                  aaTATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgt                                  <  1:3245072/71‑1 (MQ=255)
                                  aaTATCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgt                                  <  1:2589934/71‑1 (MQ=255)
                                   atatCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAAc                                                       >  1:2008863/1‑49 (MQ=255)
                                   atatCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtg                                 <  1:1330399/71‑1 (MQ=255)
                                   atatCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtg                                 >  1:1553009/1‑71 (MQ=255)
                                   atatCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgt                                  >  1:2137634/1‑70 (MQ=255)
                                    tatCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtgt                                <  1:1946532/71‑1 (MQ=255)
                                    tatCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtgt                                <  1:2361683/71‑1 (MQ=255)
                                    tatCCAGTACAGCCTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtgt                                <  1:2263040/71‑1 (MQ=255)
                                      tCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGtt                               >  1:3055477/1‑70 (MQ=255)
                                      tCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGtt                               >  1:1410306/1‑70 (MQ=255)
                                       ccAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtgt                                >  1:1030699/1‑68 (MQ=255)
                                       ccAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtca                             >  1:2546778/1‑71 (MQ=255)
                                       ccAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtca                             >  1:2715944/1‑71 (MQ=255)
                                        cAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtcat                            >  1:2784734/1‑71 (MQ=255)
                                        cAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtcat                            >  1:1134923/1‑71 (MQ=255)
                                          gTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtcatca                          <  1:3228269/71‑1 (MQ=255)
                                           tACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCaaa                                                        <  1:3676942/40‑1 (MQ=255)
                                             cAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAACTGTGCGACCCGGTGTTCATCAAtt                       >  1:1116639/1‑71 (MQ=255)
                                              tttcTGTTTCAGAGATACCGGAATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTg                      <  1:2864966/68‑1 (MQ=255)
                                              aGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTg                      >  1:2554519/1‑71 (MQ=255)
                                              aGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTg                      >  1:2029668/1‑71 (MQ=255)
                                              aGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCAAGCTGTGCGACCCGGTGTTAATCAATTg                      >  1:3597045/1‑71 (MQ=255)
                                                 cTGTTTCAGAGATACCGGGATCGGCCAACTCAAACtt                                                     <  1:2397033/37‑1 (MQ=255)
                                                 cTGTTTCAGAGATACCGGGATCGGCCAACTCAAACtt                                                     <  1:3352316/37‑1 (MQ=255)
                                                 cTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCg                   >  1:883677/1‑71 (MQ=255)
                                                  tGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGc                  >  1:1280520/1‑71 (MQ=255)
                                                       cAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca             >  1:1341258/1‑71 (MQ=255)
                                                       cAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca             >  1:3895271/1‑71 (MQ=255)
                                                        agagATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca             >  1:1512651/1‑70 (MQ=255)
                                                            aTACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacacaaca        >  1:3708977/1‑71 (MQ=255)
                                                              aCCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacacaa          <  1:816834/67‑1 (MQ=255)
                                                                 gggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCaa    >  1:340242/1‑70 (MQ=255)
                                                                 gggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCaa    >  1:2337352/1‑70 (MQ=255)
                                                                 gggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCaa    >  1:83919/1‑70 (MQ=255)
                                                                 gggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCaa    >  1:18254/1‑70 (MQ=255)
                                                                 gggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGCTGTTCATCAATTGTCGCCAACACAACATCAac   <  1:2789035/71‑1 (MQ=255)
                                                                  ggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacacaa          >  1:117311/1‑63 (MQ=255)
                                                                  ggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAaca  <  1:3889224/71‑1 (MQ=255)
                                                                  ggATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAaca  >  1:427661/1‑71 (MQ=255)
                                                                   |                                                                     
CAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACA  >  minE/1047107‑1047243

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: