Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 626,669 C→G 5.5% D194E (GAC→GAG ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE626,6690CG 5.5% 328.1 / 5.4 112D194E (GAC→GAGycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/59);  new base G (3/3);  total (49/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 116 total aligned reads displayed.

AGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGT  >  minE/626602‑626738
                                                                   |                                                                     
   tgatggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                                 <  1:226071/71‑1 (MQ=255)
   tgatggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                                 <  1:3819018/71‑1 (MQ=255)
   tgatggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                                 <  1:2225984/71‑1 (MQ=255)
      tggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgatgat                                                              <  1:369710/71‑1 (MQ=255)
      tggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgatga                                                               >  1:1903008/1‑70 (MQ=255)
        gtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgatgat                                                              <  1:432842/69‑1 (MQ=255)
        gtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGTatgatg                                                             <  1:10109/70‑1 (MQ=255)
           atTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTgg                                                                    <  1:2496167/60‑1 (MQ=255)
           atTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTgg                                                                    <  1:4441021/60‑1 (MQ=255)
           atTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgc                                                          <  1:2758017/70‑1 (MQ=255)
            ttGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgcgc                                                        <  1:491042/71‑1 (MQ=255)
              gAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGAtgatga                                                            <  1:3761772/65‑1 (MQ=255)
              gAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGa                                                      >  1:4543797/1‑71 (MQ=255)
                aCGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                                    <  1:1813816/71‑1 (MQ=255)
                aCGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                                    <  1:2693692/71‑1 (MQ=255)
                aCGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                                    <  1:4227533/71‑1 (MQ=255)
                     gcAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGGGCGATCTCGtt                                               <  1:2970791/71‑1 (MQ=255)
                     gcAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGtt                                               <  1:1237825/71‑1 (MQ=255)
                        aCTTCCAGCGCTTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgcg                                                         <  1:188129/58‑1 (MQ=255)
                         cTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTgg                                                                    <  1:2216650/46‑1 (MQ=255)
                         cTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAg                                                           >  1:4095940/1‑55 (MQ=255)
                          ttCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTg                                                                     <  1:3487940/44‑1 (MQ=255)
                           tCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCa                                          >  1:1201661/1‑70 (MQ=255)
                            ccAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGACGATGAGCGCGATCTCGTTTCGCaaa                                        >  1:2732809/1‑71 (MQ=255)
                             cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        <  1:3449168/70‑1 (MQ=255)
                             cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        <  1:4242362/70‑1 (MQ=255)
                             cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        <  1:2431103/70‑1 (MQ=255)
                             cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        >  1:3051420/1‑70 (MQ=255)
                              aGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        >  1:215403/1‑69 (MQ=255)
                               gcgcGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCa                                          <  1:4363935/66‑1 (MQ=255)
                               gcgcGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                     <  1:3484304/71‑1 (MQ=255)
                                cgcgTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                     <  1:179120/70‑1 (MQ=255)
                                cgcgTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                     >  1:3616795/1‑70 (MQ=255)
                                 gcgTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                     >  1:1905140/1‑69 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgc                                                          >  1:3970652/1‑43 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAg                                                           >  1:2258857/1‑42 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACa                                      >  1:4580616/1‑63 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                               >  1:675383/1‑70 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                               >  1:980906/1‑70 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                               >  1:3886900/1‑70 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                               >  1:1954360/1‑70 (MQ=255)
                                      tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                               >  1:1847141/1‑70 (MQ=255)
                                       aTGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              >  1:1135134/1‑70 (MQ=255)
                                       aTGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              >  1:2542816/1‑70 (MQ=255)
                                          aTTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatc                          <  1:79980/71‑1 (MQ=255)
                                           tttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatca                         <  1:2810462/71‑1 (MQ=255)
                                           tttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatca                         >  1:3410706/1‑71 (MQ=255)
                                            ttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        >  1:4423463/1‑71 (MQ=255)
                                            ttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        >  1:3163914/1‑71 (MQ=255)
                                            ttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatc                          >  1:21407/1‑69 (MQ=255)
                                             tAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        <  1:3209769/70‑1 (MQ=255)
                                             tAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        <  1:1778724/70‑1 (MQ=255)
                                             tAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        >  1:2351352/1‑70 (MQ=255)
                                              aaCCCACCGTGTCATGCCTGATTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc                       <  1:4342943/70‑1 (MQ=255)
                                              aaCCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc                       <  1:3883832/70‑1 (MQ=255)
                                                cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCa                                >  1:3596656/1‑59 (MQ=255)
                                                cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              <  1:1172046/61‑1 (MQ=255)
                                                cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              <  1:3312985/61‑1 (MQ=255)
                                                cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              <  1:292902/61‑1 (MQ=255)
                                                cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        <  1:489116/67‑1 (MQ=255)
                                                cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        <  1:2895403/67‑1 (MQ=255)
                                                 ccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        <  1:2715940/66‑1 (MQ=255)
                                                 ccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      >  1:612569/1‑68 (MQ=255)
                                                 ccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGc                   <  1:4080274/71‑1 (MQ=255)
                                                   aCCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAat                           >  1:2269405/1‑61 (MQ=255)
                                                     cGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTgg                <  1:2159232/70‑1 (MQ=255)
                                                      gtgtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              <  1:2666423/55‑1 (MQ=255)
                                                        gtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                            >  1:3927254/1‑55 (MQ=255)
                                                        gtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                              <  1:1573366/53‑1 (MQ=255)
                                                        gtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      <  1:3478437/61‑1 (MQ=255)
                                                        gtCATGCCTGACTGGGATGATGAGAGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      <  1:565100/61‑1 (MQ=255)
                                                         tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                        <  1:2743315/58‑1 (MQ=255)
                                                         tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa            <  1:3147470/70‑1 (MQ=255)
                                                         tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATa             <  1:1392419/69‑1 (MQ=255)
                                                         tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAAc           <  1:125329/71‑1 (MQ=255)
                                                          cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa            >  1:2252527/1‑69 (MQ=255)
                                                          cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa            >  1:1932846/1‑69 (MQ=255)
                                                          cATGCCTGACTGCGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa            <  1:3913948/69‑1 (MQ=255)
                                                           aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAg                                     <  1:845258/43‑1 (MQ=255)
                                                           aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg         >  1:4034851/1‑71 (MQ=38)
                                                           aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg         >  1:1512972/1‑71 (MQ=38)
                                                           aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg         >  1:3562418/1‑71 (MQ=38)
                                                           aTGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaa                                         >  1:4576149/1‑39 (MQ=255)
                                                            tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        >  1:970844/1‑39 (MQ=255)
                                                            tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaa                             >  1:3279965/1‑50 (MQ=255)
                                                            tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      >  1:2118663/1‑57 (MQ=255)
                                                            tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa            <  1:2296613/67‑1 (MQ=255)
                                                             gCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        >  1:2730561/1‑38 (MQ=255)
                                                             gCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                        >  1:334367/1‑38 (MQ=255)
                                                              ccTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAat                           <  1:717951/50‑1 (MQ=255)
                                                              ccTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTg       <  1:4281426/70‑1 (MQ=255)
                                                                tGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      >  1:1989635/1‑53 (MQ=255)
                                                                 gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                            >  1:1266206/1‑46 (MQ=255)
                                                                 gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                            >  1:623408/1‑46 (MQ=255)
                                                                 gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                            >  1:2841969/1‑46 (MQ=255)
                                                                 gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGcgcg   >  1:2323947/1‑71 (MQ=255)
                                                                   cTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGa                                    <  1:4301114/36‑1 (MQ=255)
                                                                   cTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      >  1:4581606/1‑50 (MQ=255)
                                                                   cTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                      >  1:546094/1‑50 (MQ=255)
                                                                   cTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGt  <  1:2044423/70‑1 (MQ=255)
                                                                   |                                                                     
AGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGT  >  minE/626602‑626738

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: