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

GAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGA  >  minE/975189‑975327
                                                                     |                                                                      
gAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑a                                                                       >  1:352161/1‑70 (MQ=255)
gAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑At                                                                      >  1:3532200/1‑71 (MQ=255)
 aaCCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑aat                                                                     >  1:3562574/1‑69 (MQ=255)
      gacgacCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGtt                                                                 <  1:335083/71‑1 (MQ=255)
       acgacCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTaa                                                               >  1:2570394/1‑71 (MQ=255)
          acCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGc                                                            <  1:167766/71‑1 (MQ=255)
            ccGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGcc                                                           >  1:2123206/1‑71 (MQ=255)
              ggCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCg                                                          >  1:2801930/1‑69 (MQ=255)
               gCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATc                                                                     >  1:1163501/1‑57 (MQ=255)
                 aaTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTa                                                     <  1:2510332/71‑1 (MQ=255)
                  aTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTa                                                     >  1:4109807/1‑70 (MQ=255)
                  aTCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTa                                                     >  1:3367122/1‑70 (MQ=255)
                   tCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTaaa                                                   <  1:1118899/71‑1 (MQ=255)
                    cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTaa                                                               >  1:1863951/1‑58 (MQ=255)
                    cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAg                                                       >  1:701782/1‑66 (MQ=255)
                    cccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTaaa                                                   >  1:2145325/1‑70 (MQ=255)
                     ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑a                                                                       >  1:1937369/1‑49 (MQ=255)
                     ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGt                                                                  >  1:804183/1‑54 (MQ=255)
                     ccGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGt                                                                  >  1:2943073/1‑54 (MQ=255)
                      cGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaaga                                                <  1:431374/71‑1 (MQ=255)
                      cGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaaga                                                <  1:823662/71‑1 (MQ=255)
                       gTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaaga                                                >  1:82762/1‑70 (MQ=255)
                       gTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaaga                                                >  1:75011/1‑70 (MQ=255)
                        tCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAAt                                                              <  1:3264497/56‑1 (MQ=39)
                        tCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaagaa                                               <  1:407347/71‑1 (MQ=255)
                        tCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaag                                                 <  1:3646041/68‑1 (MQ=255)
                          aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaagaaag                                             >  1:199698/1‑70 (MQ=255)
                          aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGc                                            <  1:805298/71‑1 (MQ=255)
                          aaaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGc                                            <  1:2016529/71‑1 (MQ=255)
                           aaaCATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGt                                                                         <  1:107772/42‑1 (MQ=255)
                              cATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:3101824/1‑69 (MQ=255)
                               aTCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTaaaa                                                  <  1:3254364/60‑1 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑a                                                                       <  1:3532094/38‑1 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCg                                                                    >  1:1712360/1‑41 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCg                                                                    >  1:2312501/1‑41 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGtt                                                                 <  1:1479661/44‑1 (MQ=255)
                                tCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTaa                                                               >  1:3247357/1‑46 (MQ=255)
                                   atatTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaag                                                 <  1:1180002/57‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                                <  1:2402999/58‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAaaaga                                                <  1:633079/58‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGtt                                                                 <  1:6352/40‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGtt                                                                 <  1:282496/40‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGtt                                                                 <  1:1456259/40‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTaa                                                    <  1:959395/53‑1 (MQ=255)
                                    tatTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGc                                            <  1:1319666/61‑1 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCcctcct                                  >  1:1873774/1‑70 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCg                                                                    >  1:3693220/1‑35 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTaa                                                               >  1:135282/1‑40 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTa                                                                >  1:3661555/1‑39 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaagaa                                               >  1:1771534/1‑56 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:667926/1‑61 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTg                                 >  1:1773263/1‑70 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                <  1:179374/71‑1 (MQ=255)
                                      ttGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                <  1:1722203/71‑1 (MQ=255)
                                        gCTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATACCGCAGTAAAAGAAAGCAACCCCTCCTGa                                >  1:244399/1‑70 (MQ=255)
                                        gCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGtt                                                                 >  1:3522991/1‑36 (MQ=255)
                                        gCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAAcc                                        >  1:1918351/1‑61 (MQ=255)
                                         cTACGGATCGTCAGCAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                >  1:3116604/1‑69 (MQ=255)
                                         cTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAaaag                                                 <  1:3311559/51‑1 (MQ=255)
                                          tACGGATCGTCAGCAAATAAATTTTGTAATCGGTTaa                                                               >  1:1419218/1‑37 (MQ=37)
                                            cGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGc                           <  1:364552/70‑1 (MQ=255)
                                            cGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCATCCTGa                                <  1:1659/65‑1 (MQ=255)
                                               aTCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAATCaa                                          >  1:3745843/1‑52 (MQ=255)
                                                tCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAgg                            >  1:2157865/1‑65 (MQ=255)
                                                tCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                      <  1:2213885/71‑1 (MQ=255)
                                                    cAGCAAATAAATTTTTT‑ATCGGTTAATGCCGCAGTaaaaca                                                <  1:3923443/41‑3 (MQ=38)
                                                     aGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATg                  <  1:2728124/70‑1 (MQ=255)
                                                       cAAATAAATTTTGTAATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGa                                <  1:357755/55‑1 (MQ=255)
                                                       cAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTaaaa                                                  <  1:3330947/36‑1 (MQ=255)
                                                       cAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTaaaa                                                  >  1:3928041/1‑36 (MQ=255)
                                                       cAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCaa                                          >  1:3998983/1‑44 (MQ=255)
                                                       cAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc               >  1:3484274/1‑71 (MQ=255)
                                                       cAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATgcgc               >  1:839929/1‑71 (MQ=255)
                                                         aataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGt             >  1:2669953/1‑71 (MQ=255)
                                                          ataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCcct                                     <  1:1201145/46‑1 (MQ=255)
                                                          ataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGt                         >  1:2950555/1‑58 (MQ=255)
                                                          ataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa                    <  1:2371072/63‑1 (MQ=255)
                                                          ataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa                    <  1:1599181/63‑1 (MQ=255)
                                                          ataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGt             <  1:1567351/70‑1 (MQ=255)
                                                          ataaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTc            >  1:2607781/1‑71 (MQ=255)
                                                           taaatTTTGT‑ATCGGTTAATGCCGCAGTAaaagaa                                               >  1:1827088/1‑35 (MQ=255)
                                                           taaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTcc           >  1:1140163/1‑71 (MQ=255)
                                                           taaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTcc           >  1:2082436/1‑71 (MQ=255)
                                                           taaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTcc           >  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)
                                                            aaatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAAccc                                       >  1:1181111/1‑42 (MQ=255)
                                                             aatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAAccc                                       >  1:2200019/1‑41 (MQ=255)
                                                             aatTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCt          <  1:2540613/70‑1 (MQ=255)
                                                                 ttGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGaaa                    <  1:3972727/56‑1 (MQ=255)
                                                                 ttGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAg     >  1:2390880/1‑71 (MQ=255)
                                                                    t‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGa                      <  1:3098948/51‑1 (MQ=255)
                                                                    t‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGa  <  1:892342/71‑1 (MQ=255)
                                                                     |                                                                      
GAACCCGACGACCCGGCAATCCCGTCAAAACATCCATATTGCTACGGATCGTCAGCAAATAAATTTTGT‑ATCGGTTAATGCCGCAGTAAAAGAAAGCAACCCCTCCTGAAAGGCGTCGAAATGCGCGTCCTGCCAGTGA  >  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: ‑