Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 975,257:1 +A 17.0% coding (358/375 nt) ydeH ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE975,2571.A17.0% 261.8 / 69.2 88coding (358/375 nt)ydeHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (42/31);  new base A (9/6);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGA  >  minE/975189‑975327
                                                                     |                                                                      
gAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTa                                                                       >  1:352161/1‑70 (MQ=255)
gAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAt                                                                      >  1:3532200/1‑71 (MQ=255)
 aaCCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTaat                                                                     >  1:3562574/1‑69 (MQ=255)
      gacgacCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGtt                                                                 <  1:335083/71‑1 (MQ=255)
       acgacCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTaa                                                               >  1:2570394/1‑71 (MQ=255)
          acCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGc                                                            <  1:167766/71‑1 (MQ=255)
            ccGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGcc                                                           >  1:2123206/1‑71 (MQ=255)
              ggCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCg                                                          >  1:2801930/1‑69 (MQ=255)
               gCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATc                                                                     >  1:1163501/1‑57 (MQ=255)
                 aaTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTa                                                     <  1:2510332/71‑1 (MQ=255)
                  aTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTa                                                     >  1:4109807/1‑70 (MQ=255)
                  aTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTa                                                     >  1:3367122/1‑70 (MQ=255)
                   tCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTaaa                                                   <  1:1118899/71‑1 (MQ=255)
                    cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTaa                                                               >  1:1863951/1‑58 (MQ=255)
                    cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAg                                                       >  1:701782/1‑66 (MQ=255)
                    cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTaaa                                                   >  1:2145325/1‑70 (MQ=255)
                     ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTa                                                                       >  1:1937369/1‑49 (MQ=255)
                     ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGt                                                                  >  1:804183/1‑54 (MQ=255)
                     ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGt                                                                  >  1:2943073/1‑54 (MQ=255)
                      cGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaaga                                                <  1:431374/71‑1 (MQ=255)
                      cGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaaga                                                <  1:823662/71‑1 (MQ=255)
                       gTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaaga                                                >  1:82762/1‑70 (MQ=255)
                       gTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaaga                                                >  1:75011/1‑70 (MQ=255)
                        tCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAAt                                                              <  1:3264497/56‑1 (MQ=39)
                        tCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaagaa                                               <  1:407347/71‑1 (MQ=255)
                        tCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaag                                                 <  1:3646041/68‑1 (MQ=255)
                          aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaagaaag                                             >  1:199698/1‑70 (MQ=255)
                          aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGc                                            <  1:805298/71‑1 (MQ=255)
                          aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGc                                            <  1:2016529/71‑1 (MQ=255)
                           aaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGt                                                                         <  1:107772/42‑1 (MQ=255)
                              cATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:3101824/1‑69 (MQ=255)
                               aTCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTaaaa                                                  <  1:3254364/60‑1 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTa                                                                       <  1:3532094/38‑1 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCg                                                                    >  1:1712360/1‑41 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCg                                                                    >  1:2312501/1‑41 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                                 <  1:1479661/44‑1 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTaa                                                               >  1:3247357/1‑46 (MQ=255)
                                   atatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaag                                                 <  1:1180002/57‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                                <  1:2402999/58‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                                <  1:633079/58‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                                 <  1:6352/40‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                                 <  1:282496/40‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                                 <  1:1456259/40‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTaa                                                    <  1:959395/53‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGc                                            <  1:1319666/61‑1 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCcctcct                                  >  1:1873774/1‑70 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCg                                                                    >  1:3693220/1‑35 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTaa                                                               >  1:135282/1‑40 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTa                                                                >  1:3661555/1‑39 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaagaa                                               >  1:1771534/1‑56 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:667926/1‑61 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTg                                 >  1:1773263/1‑70 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                <  1:179374/71‑1 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                <  1:1722203/71‑1 (MQ=255)
                                        gCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATACCGCAGTAAAAGAAAGCAACCCCTCCTGa                                >  1:244399/1‑70 (MQ=255)
                                        gCTACGGATCGTCAGCAAATAAATTTTGTATCGGtt                                                                 >  1:3522991/1‑36 (MQ=255)
                                        gCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAAcc                                        >  1:1918351/1‑61 (MQ=255)
                                         cTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                >  1:3116604/1‑69 (MQ=255)
                                         cTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAaaag                                                 <  1:3311559/51‑1 (MQ=255)
                                          tACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                               >  1:1419218/1‑37 (MQ=37)
                                            cGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGc                           <  1:364552/70‑1 (MQ=255)
                                            cGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCATCCTGa                                <  1:1659/65‑1 (MQ=255)
                                               aTCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAATCaa                                          >  1:3745843/1‑52 (MQ=255)
                                                tCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAgg                            >  1:2157865/1‑65 (MQ=255)
                                                tCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                      <  1:2213885/71‑1 (MQ=255)
                                                    cAGCAAATAAATTTTTTATCGGTTAATGCCGCAGTaaaaca                                                <  1:3923443/41‑3 (MQ=38)
                                                     aGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATg                  <  1:2728124/70‑1 (MQ=255)
                                                       cAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                <  1:357755/55‑1 (MQ=255)
                                                       cAAATAAATTTTGTATCGGTTAATGCCGCAGTaaaa                                                  <  1:3330947/36‑1 (MQ=255)
                                                       cAAATAAATTTTGTATCGGTTAATGCCGCAGTaaaa                                                  >  1:3928041/1‑36 (MQ=255)
                                                       cAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:3998983/1‑44 (MQ=255)
                                                       cAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc               >  1:3484274/1‑71 (MQ=255)
                                                       cAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc               >  1:839929/1‑71 (MQ=255)
                                                         aataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGt             >  1:2669953/1‑71 (MQ=255)
                                                          ataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCcct                                     <  1:1201145/46‑1 (MQ=255)
                                                          ataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGt                         >  1:2950555/1‑58 (MQ=255)
                                                          ataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa                    <  1:2371072/63‑1 (MQ=255)
                                                          ataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa                    <  1:1599181/63‑1 (MQ=255)
                                                          ataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGt             <  1:1567351/70‑1 (MQ=255)
                                                          ataaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTc            >  1:2607781/1‑71 (MQ=255)
                                                           taaatTTTGTATCGGTTAATGCCGCAGTAaaagaa                                               >  1:1827088/1‑35 (MQ=255)
                                                           taaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTcc           >  1:1140163/1‑71 (MQ=255)
                                                           taaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTcc           >  1:2082436/1‑71 (MQ=255)
                                                           taaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTcc           >  1:3894385/1‑71 (MQ=255)
                                                            aaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:420448/1‑40 (MQ=38)
                                                            aaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:294780/1‑40 (MQ=38)
                                                            aaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTc            >  1:3621409/1‑70 (MQ=255)
                                                            aaatTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTc            >  1:3387661/1‑70 (MQ=255)
                                                            aaatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAAccc                                       >  1:1181111/1‑42 (MQ=255)
                                                             aatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAAccc                                       >  1:2200019/1‑41 (MQ=255)
                                                             aatTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCt          <  1:2540613/70‑1 (MQ=255)
                                                                 ttGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa                    <  1:3972727/56‑1 (MQ=255)
                                                                 ttGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAg     >  1:2390880/1‑71 (MQ=255)
                                                                    tATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                      <  1:3098948/51‑1 (MQ=255)
                                                                    tATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGa  <  1:892342/71‑1 (MQ=255)
                                                                     |                                                                      
GAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGA  >  minE/975189‑975327

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: