Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 626,673 G→A 5.8% D196N (GAT→AAT)  ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE626,6730GA 5.8% 322.2 / 6.4 104D196N (GAT→AAT) ycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (47/50);  new base A (3/3);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 113 total aligned reads displayed.

TGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCC  >  minE/626605‑626742
                                                                    |                                                                     
tgatggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                                     <  1:3819018/71‑1 (MQ=255)
tgatggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                                     <  1:2225984/71‑1 (MQ=255)
tgatggtgatTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                                     <  1:226071/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)
        atTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgc                                                              <  1:2758017/70‑1 (MQ=255)
         ttGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgcgc                                                            <  1:491042/71‑1 (MQ=255)
          tGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgcgc                                                            >  1:4401394/1‑70 (MQ=255)
           gAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGAtgatga                                                                <  1:3761772/65‑1 (MQ=255)
             aCGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                                        <  1:1813816/71‑1 (MQ=255)
             aCGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                                        <  1:2693692/71‑1 (MQ=255)
                  gcAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGGGCGATCTCGtt                                                   <  1:2970791/71‑1 (MQ=255)
                  gcAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGtt                                                   <  1:1237825/71‑1 (MQ=255)
                    aaCTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTc                                                 >  1:3434084/1‑71 (MQ=255)
                      cTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAg                                                               >  1:4095940/1‑55 (MQ=255)
                        tCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCa                                              >  1:1201661/1‑70 (MQ=255)
                         ccAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGACGATGAGCGCGATCTCGTTTCGCaaa                                            >  1:2732809/1‑71 (MQ=255)
                          cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                            >  1:3051420/1‑70 (MQ=255)
                          cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                            <  1:3449168/70‑1 (MQ=255)
                          cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                            <  1:2431103/70‑1 (MQ=255)
                          cAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                            <  1:4242362/70‑1 (MQ=255)
                           aGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                            >  1:215403/1‑69 (MQ=255)
                           aGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACa                                          <  1:1973778/71‑1 (MQ=255)
                            gcgcGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCa                                              <  1:4363935/66‑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)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgc                                                              >  1:3970652/1‑43 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgc                                                              >  1:631150/1‑43 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAg                                                               >  1:2258857/1‑42 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACa                                          >  1:4580616/1‑63 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                   >  1:3886900/1‑70 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                   >  1:675383/1‑70 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                   >  1:980906/1‑70 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                   >  1:1847141/1‑70 (MQ=255)
                                   tATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                                   >  1:1954360/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)
                                          tAACTCACCGTGTCATGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAg                                         <  1:1840279/57‑1 (MQ=255)
                                          tAACTCACCGTGTCATGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAatcat                            <  1:2784506/70‑1 (MQ=255)
                                          tAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                            >  1:2351352/1‑70 (MQ=255)
                                          tAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                            <  1:1778724/70‑1 (MQ=255)
                                          tAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                            <  1:3209769/70‑1 (MQ=255)
                                           aaCCCACCGTGTCATGCCTGATTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc                           <  1:4342943/70‑1 (MQ=255)
                                           aaCCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc                           <  1:3883832/70‑1 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCa                                    >  1:3601949/1‑59 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCa                                    >  1:3596656/1‑59 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                  <  1:3312985/61‑1 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                  <  1:1172046/61‑1 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                  <  1:292902/61‑1 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                            <  1:2895403/67‑1 (MQ=255)
                                             cccACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat                            <  1:489116/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)
                                                  cGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTgg                    <  1:2159232/70‑1 (MQ=255)
                                                   gtgtCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                                  <  1:2666423/55‑1 (MQ=255)
                                                    tgtCATGCCTGACTGGGATGATGAGCGCGATCTCGtt                                                   <  1:1275705/37‑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)
                                                      tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAAc               <  1:125329/71‑1 (MQ=255)
                                                      tCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAAc               <  1:4147106/71‑1 (MQ=255)
                                                       cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                >  1:2252527/1‑69 (MQ=255)
                                                       cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                >  1:1932846/1‑69 (MQ=255)
                                                       cATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATa                 <  1:2185536/68‑1 (MQ=255)
                                                       cATGCCTGACTGCGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATaa                <  1:3913948/69‑1 (MQ=255)
                                                        aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAg                                         <  1:845258/43‑1 (MQ=255)
                                                        aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg             >  1:1512972/1‑71 (MQ=38)
                                                        aTGCCTGAGTGGAATGATGAGCGCGATCTCGTTTCGCAAGCAGAAGCAGAAATCATCATGCTGGATAACAg             >  1:4034851/1‑71 (MQ=38)
                                                        aTGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaa                                             >  1:2312185/1‑39 (MQ=255)
                                                        aTGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaa                                             >  1:4576149/1‑39 (MQ=255)
                                                         tGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                                            >  1:970844/1‑39 (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)
                                                              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)
                                                                   ggATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAg    <  1:1830011/69‑1 (MQ=255)
                                                                   ggATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGcc  <  1:3300234/71‑1 (MQ=255)
                                                                   ggATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGc   <  1:230713/70‑1 (MQ=255)
                                                                    |                                                                     
TGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGATAACAGTGCGCGTAGCC  >  minE/626605‑626742

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: