Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,124,793 C→T 6.4% Q70Q (CAG→CAA yhbP ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,124,7930CT 6.4% 202.5 / 3.7 62Q70Q (CAG→CAAyhbPconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/23);  new base T (2/2);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGTTTTTTCT  >  minE/2124725‑2124859
                                                                    |                                                                  
tCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCt                                                                   >  1:470578/1‑70 (MQ=255)
tCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCt                                                                   >  1:519810/1‑70 (MQ=255)
      ccTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            <  1:1085927/71‑1 (MQ=255)
      ccTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            <  1:651874/71‑1 (MQ=255)
      ccTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            <  1:665424/71‑1 (MQ=255)
       cTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGTTGACCGtt                                                            >  1:1278609/1‑70 (MQ=255)
       cTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            >  1:3772731/1‑70 (MQ=255)
       cTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            >  1:2393027/1‑70 (MQ=255)
       cTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             >  1:2271971/1‑69 (MQ=255)
          ccAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTA‑‑GCTACCGTTTTCGGTTGACCGTTTACCGt                                                      >  1:913345/1‑71 (MQ=255)
                  cGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCgg                                                 >  1:1233423/1‑70 (MQ=255)
                   ggATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGc                                                >  1:1386188/1‑70 (MQ=255)
                   ggATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGc                                                >  1:519489/1‑70 (MQ=255)
                   ggATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCa                                               <  1:3788212/71‑1 (MQ=255)
                     aTCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             <  1:4245434/55‑1 (MQ=255)
                     aTCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             <  1:1999298/55‑1 (MQ=255)
                     aTCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             <  1:425798/55‑1 (MQ=255)
                     aTCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             <  1:918424/55‑1 (MQ=255)
                      tctcACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAAcc                                            >  1:652232/1‑71 (MQ=255)
                      tctcACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAAcc                                            >  1:1275830/1‑71 (MQ=255)
                       ctcACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCg                                           >  1:2165237/1‑71 (MQ=255)
                       ctcACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCg                                           >  1:2978978/1‑71 (MQ=255)
                         cACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCg                                           >  1:731510/1‑69 (MQ=255)
                         cACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCg                                           >  1:301996/1‑69 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGTCTGACCGtt                                                            >  1:688054/1‑51 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGTTGACCGTTTACCGTTCCGGCAACCGCAg                                        <  1:2547781/71‑1 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGTTGACCGTTTACCGTTCCGGCAACCGCAg                                        <  1:378144/71‑1 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            >  1:465674/1‑51 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGtt                                                            >  1:1114873/1‑51 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             >  1:2819316/1‑50 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             >  1:3299760/1‑50 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             >  1:3848130/1‑50 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGt                                                             >  1:1678486/1‑50 (MQ=255)
                          aCCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCa                                               <  1:1971442/64‑1 (MQ=255)
                               tAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCaa                                              >  1:4249625/1‑60 (MQ=255)
                               tAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCaa                                              >  1:1059211/1‑60 (MQ=255)
                                        cacCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTc                          <  1:3484463/71‑1 (MQ=255)
                                          ccGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGc                                                <  1:3498283/47‑1 (MQ=255)
                                             cgAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTg                     <  1:1177067/71‑1 (MQ=255)
                                             cgAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTg                     >  1:3022791/1‑71 (MQ=255)
                                              gAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTgc                    <  1:3752954/71‑1 (MQ=255)
                                                aTTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTc                          <  1:4199154/63‑1 (MQ=255)
                                                aTTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTgcgc                  <  1:1339700/71‑1 (MQ=255)
                                                aTTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTgcgc                  <  1:452558/71‑1 (MQ=255)
                                                 ttAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGc                                                <  1:3647475/40‑1 (MQ=255)
                                                   aaCGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCaa                                              >  1:1214297/1‑40 (MQ=255)
                                                   aaCGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCaa                                              >  1:2139687/1‑40 (MQ=255)
                                                      gCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCACCCGAAGCCTGCGGCCCACTCATCTGCGCGTGgcg            >  1:3865806/1‑71 (MQ=255)
                                                      gCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGgcg            >  1:3138212/1‑71 (MQ=255)
                                                      gCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGAAGCCTGCGGCCCACTCATCTGCGCGTGgcg            >  1:3861132/1‑71 (MQ=255)
                                                          ccGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCa                             >  1:334751/1‑50 (MQ=255)
                                                          ccGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCa                             >  1:1863859/1‑50 (MQ=255)
                                                          ccGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCa                             >  1:4061844/1‑50 (MQ=255)
                                                          ccGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGt         <  1:4250802/70‑1 (MQ=255)
                                                            gTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGtttt      >  1:2942537/1‑71 (MQ=255)
                                                            gTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGtttt      <  1:2894605/71‑1 (MQ=255)
                                                            gTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGtttt      <  1:1401409/71‑1 (MQ=255)
                                                            gTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGtttt      >  1:3867549/1‑71 (MQ=255)
                                                               ttCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGtttttt    >  1:1787616/1‑70 (MQ=255)
                                                                tCGGCTGACCGTTTACCGTTCCGGCAACCCCAGCCTGCGGCCCACTCATCTg                     <  1:326420/52‑1 (MQ=255)
                                                                 cGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGTTTTttct  <  1:771674/70‑1 (MQ=255)
                                                                 cGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGTTTTttct  <  1:96647/70‑1 (MQ=255)
                                                                    |                                                                  
TCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTCCGGCAACCGCAGCCTGCGGCCCACTCATCTGCGCGTGGCGCGTTTTTTCT  >  minE/2124725‑2124859

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: