Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 854,861 C→T 19.7% A43A (GCG→GCA yciG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE854,8610CT19.7% 143.7 / 27.9 71A43A (GCG→GCAyciGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base T (6/8);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCT  >  minE/854792‑854927
                                                                     |                                                                  
ttAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCg                                                                   >  1:345959/1‑71 (MQ=255)
  aataatTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCt                                                                 >  1:3677709/1‑71 (MQ=255)
      atTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:4757041/70‑1 (MQ=255)
      atTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:4980622/70‑1 (MQ=255)
      atTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAg                                                             <  1:257891/71‑1 (MQ=255)
      atTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCACCTTTTTTACCCGCTTCa                                                              <  1:2275423/70‑1 (MQ=255)
              tGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:3061500/62‑1 (MQ=255)
                aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                                   >  1:5177215/1‑71 (MQ=255)
                aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                                   <  1:2085009/71‑1 (MQ=255)
                aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                                   >  1:1928582/1‑71 (MQ=255)
                  tCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCg                                                 <  1:4809467/71‑1 (MQ=255)
                   cAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:4094745/57‑1 (MQ=255)
                   cAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:2758171/57‑1 (MQ=255)
                   cAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:4607183/57‑1 (MQ=255)
                   cAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:4616181/57‑1 (MQ=255)
                     ggATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                            >  1:204844/1‑57 (MQ=255)
                      gATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                                   <  1:5578218/65‑1 (MQ=255)
                        tttCCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAt                                           <  1:706094/71‑1 (MQ=255)
                          tGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                                         >  1:1673781/1‑71 (MQ=255)
                            ccTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCg                                                 <  1:1303607/61‑1 (MQ=255)
                              tGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:1515083/46‑1 (MQ=255)
                              tGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              <  1:5152521/46‑1 (MQ=255)
                              tGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCgg                                                >  1:1222567/1‑60 (MQ=255)
                              tGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                      >  1:1199558/1‑70 (MQ=255)
                                aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                     <  1:4841389/69‑1 (MQ=255)
                                   tattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                              >  1:785152/1‑41 (MQ=255)
                                   tattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                >  1:4708965/1‑71 (MQ=255)
                                     ttaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                               <  1:4865645/70‑1 (MQ=255)
                                        ccaccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc                            <  1:1450435/70‑1 (MQ=255)
                                          accGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                         <  1:4981776/71‑1 (MQ=255)
                                           ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                            >  1:5148255/1‑35 (MQ=255)
                                           ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                            >  1:3582631/1‑35 (MQ=255)
                                           ccGCTTTGTTGACCGCCTTTTTTACACGCTTCAGa                                                            >  1:3891033/1‑35 (MQ=255)
                                            cGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                       <  1:2394448/71‑1 (MQ=255)
                                              cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATg                     >  1:1814136/1‑71 (MQ=255)
                                                 tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCt                   >  1:599170/1‑70 (MQ=255)
                                                    ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTAccccc                           >  1:2537712/3‑59 (MQ=255)
                                                     taccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca              >  1:5201179/2‑70 (MQ=255)
                                                     gACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                       <  1:262660/62‑1 (MQ=255)
                                                      accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCAtt                                          >  1:4085450/1‑42 (MQ=25)
                                                      accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTgtc                <  1:2902977/68‑3 (MQ=255)
                                                      accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTccac             >  1:2807254/1‑69 (MQ=255)
                                                      accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca              <  1:1115050/70‑2 (MQ=255)
                                                      accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca              <  1:2408153/70‑2 (MQ=255)
                                                      accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca              <  1:2414555/70‑2 (MQ=255)
                                                      aCCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGAcctc             >  1:574382/1‑69 (MQ=255)
                                                       ccGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                              >  1:2708169/1‑37 (MQ=255)
                                                         gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGTATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt          <  1:3742697/71‑1 (MQ=255)
                                                         gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt          <  1:3456577/71‑1 (MQ=255)
                                                         gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt          <  1:49739/71‑1 (MQ=255)
                                                          ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt         >  1:504233/1‑71 (MQ=255)
                                                           cTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCttt         <  1:1081837/70‑1 (MQ=255)
                                                            ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGc             <  1:3452642/65‑1 (MQ=255)
                                                            ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTAc       >  1:226188/1‑71 (MQ=255)
                                                             tttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                  <  1:3374548/43‑1 (MQ=255)
                                                              ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTa                                >  1:223896/1‑44 (MQ=255)
                                                              ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt       <  1:1200742/69‑3 (MQ=255)
                                                              ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt       <  1:2860583/69‑3 (MQ=255)
                                                              ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt       <  1:3709759/69‑3 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                   >  1:2976150/1‑41 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                      >  1:3626140/1‑38 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                      >  1:2882276/1‑38 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                >  1:4983566/1‑44 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                >  1:4627231/1‑44 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                >  1:4340307/1‑44 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                >  1:3777310/1‑44 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                               >  1:4608047/1‑45 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                               >  1:1949034/1‑45 (MQ=255)
                                                              ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACgg     >  1:1301709/1‑71 (MQ=255)
                                                                 tACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt          <  1:1418166/63‑1 (MQ=255)
                                                                  aCCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCt  <  1:2288744/70‑1 (MQ=255)
                                                                     cGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                       >  1:276554/1‑46 (MQ=255)
                                                                     |                                                                  
TTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCT  >  minE/854792‑854927

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: