Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 626,700 A→G 6.3% T205A (ACA→GCA)  ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE626,7000AG 6.3% 392.0 / 13.2 143T205A (ACA→GCA) ycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (62/72);  new base G (3/6);  total (65/78)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

CCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGTGGCTGG  >  minE/626630‑626769
                                                                      |                                                                     
ccAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGACGATGAGCGCGATCTCGTTTCGCaaa                                                                       >  1:2732809/1‑71 (MQ=255)
 cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                                                       <  1:2431103/70‑1 (MQ=255)
 cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                                                       >  1:3051420/1‑70 (MQ=255)
 cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                                                       <  1:3449168/70‑1 (MQ=255)
  aGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                                                       >  1:215403/1‑69 (MQ=255)
  aGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACa                                                                     <  1:1973778/71‑1 (MQ=255)
   gcgcGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                                                    <  1:3484304/71‑1 (MQ=255)
    cgcgTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                                                    >  1:3616795/1‑70 (MQ=255)
    cgcgTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                                                    <  1:179120/70‑1 (MQ=255)
     gcgTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                                                                    >  1:1905140/1‑69 (MQ=255)
          tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                                              >  1:3886900/1‑70 (MQ=255)
          tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                                              >  1:980906/1‑70 (MQ=255)
          tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                                              >  1:1954360/1‑70 (MQ=255)
          tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                                              >  1:1847141/1‑70 (MQ=255)
           aTGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                                             >  1:2542816/1‑70 (MQ=255)
           aTGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                                             >  1:1135134/1‑70 (MQ=255)
               tttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatca                                                        <  1:2810462/71‑1 (MQ=255)
                ttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                                                       >  1:4423463/1‑71 (MQ=255)
                ttAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                                                       >  1:3163914/1‑71 (MQ=255)
                 tAACTCACCGTGTCATGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAg                                                                    <  1:1840279/57‑1 (MQ=255)
                 tAACTCACCGTGTCATGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAatcat                                                       <  1:2784506/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:292902/61‑1 (MQ=255)
                    cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                                             <  1:3312985/61‑1 (MQ=255)
                    cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                                             <  1:1172046/61‑1 (MQ=255)
                    cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                                                       <  1:489116/67‑1 (MQ=255)
                     ccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                                                     >  1:612569/1‑68 (MQ=255)
                     ccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGc                                                  <  1:4080274/71‑1 (MQ=255)
                         cGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTgg                                               <  1:2159232/70‑1 (MQ=255)
                          gtgtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                                             <  1:2666423/55‑1 (MQ=255)
                            gtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                                             <  1:1573366/53‑1 (MQ=255)
                            gtCATGCCTGACTGGGATGATGAGAGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                                                     <  1:565100/61‑1 (MQ=255)
                             tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           <  1:3147470/70‑1 (MQ=255)
                             tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATa                                            <  1:1392419/69‑1 (MQ=255)
                             tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAAc                                          <  1:4147106/71‑1 (MQ=255)
                              cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           >  1:1932846/1‑69 (MQ=255)
                              cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           >  1:2252527/1‑69 (MQ=255)
                              cATGCCTGACTGCGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           <  1:3913948/69‑1 (MQ=255)
                               aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAg                                                                    <  1:845258/43‑1 (MQ=255)
                               aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg                                        >  1:3562418/1‑71 (MQ=38)
                               aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg                                        >  1:1512972/1‑71 (MQ=38)
                                tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca                                                     >  1:2118663/1‑57 (MQ=255)
                                tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           <  1:2296613/67‑1 (MQ=255)
                                 gCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                                                       >  1:334367/1‑38 (MQ=255)
                                 gCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                                                       >  1:2730561/1‑38 (MQ=255)
                                     gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                                                           >  1:1266206/1‑46 (MQ=255)
                                     gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                                                           >  1:623408/1‑46 (MQ=255)
                                     gACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa                                                           >  1:2841969/1‑46 (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)
                                          ggATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGcc                             <  1:3300234/71‑1 (MQ=255)
                                          ggATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGc                              <  1:230713/70‑1 (MQ=255)
                                             tgatgaGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATcatcat                                                    <  1:3928808/45‑1 (MQ=255)
                                             tgatgaGCGCGATCTCGTTTCGCAAACAGAAGCAGaa                                                            <  1:1808275/37‑1 (MQ=255)
                                             tgatgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc                                                      <  1:2616816/43‑1 (MQ=255)
                                              gatgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCtcg                          >  1:1230292/1‑68 (MQ=255)
                                              gatgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTgg                          >  1:2055527/1‑70 (MQ=255)
                                              gatgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTgg                          >  1:3707288/1‑70 (MQ=255)
                                               atgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAat                                                          >  1:1736457/1‑37 (MQ=255)
                                               atgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAat                                                          >  1:719798/1‑37 (MQ=255)
                                               atgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATcatcat                                                    <  1:750185/43‑1 (MQ=255)
                                               atgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATcatcat                                                    <  1:219156/43‑1 (MQ=255)
                                               atgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           <  1:689913/52‑1 (MQ=255)
                                               atgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTggg                         <  1:44254/70‑1 (MQ=255)
                                                tgaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGaa                       >  1:3508993/1‑71 (MQ=255)
                                                 gaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAAc                                          >  1:945866/1‑51 (MQ=255)
                                                 gaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACa                                         >  1:1353276/1‑52 (MQ=255)
                                                 gaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACa                                         >  1:2904671/1‑52 (MQ=255)
                                                 gaGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACa                                         >  1:3106391/1‑52 (MQ=255)
                                                  aGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTg                                                >  1:3735133/1‑44 (MQ=255)
                                                  aGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAata                     <  1:122413/71‑1 (MQ=255)
                                                   gcgcGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATa                                            >  1:1591787/1‑47 (MQ=255)
                                                        aTCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAata                     <  1:3826811/65‑1 (MQ=255)
                                                        aTCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATAtt                   <  1:4654179/67‑1 (MQ=255)
                                                          ctcGTTTCGCAAACAGAAGCAGAAATCATCATGCTgg                                               <  1:4008794/37‑1 (MQ=255)
                                                          ctcGTTTCGCAAACAGAAGCAGAAATCATCATGCTg                                                >  1:912539/1‑36 (MQ=255)
                                                          ctcGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           >  1:2400173/1‑41 (MQ=255)
                                                            cGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                                           <  1:4526135/39‑1 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGt                                 >  1:2710288/1‑48 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAata                     >  1:2969140/1‑60 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCcgc               >  1:412266/1‑66 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGaa            >  1:2264450/1‑69 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGaa            >  1:4256186/1‑69 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGaa            <  1:1875589/69‑1 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGaa            <  1:95280/69‑1 (MQ=255)
                                                             gTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGaa            >  1:1257196/1‑69 (MQ=255)
                                                              tttCGCAAACAGAAGCAGAAATCATCATGCTGGATAAc                                          <  1:1680711/38‑1 (MQ=255)
                                                              tttCGCAAACAGAAGCAGAAATCATCATGCTGGATAACa                                         >  1:3786164/1‑39 (MQ=255)
                                                                tCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGtg       >  1:4067363/1‑71 (MQ=255)
                                                                 cGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAg         >  1:4548433/1‑68 (MQ=255)
                                                                  gCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGtg       <  1:3420503/69‑1 (MQ=255)
                                                                  gCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGtg       <  1:3244441/69‑1 (MQ=255)
                                                                     aaCAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGt        >  1:3260793/1‑65 (MQ=255)
                                                                      aCAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGtggctgg  >  1:2096619/1‑70 (MQ=255)
                                                                      aCAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGtggctgg  >  1:2002684/1‑70 (MQ=255)
                                                                      |                                                                     
CCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCCTGGGAATATTCCGCGAACAGTGGCTGG  >  minE/626630‑626769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: