Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,106 G→T 70.3% E31D (GAG→GAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1060GT70.3% 107.9 / 110.3 118E31D (GAG→GATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/16);  new base T (36/47);  total (55/63)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCA  >  minE/218036‑218175
                                                                      |                                                                     
 aagaagAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag                                                                       >  1:1398727/1‑70 (MQ=255)
  agaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga                                                                      <  1:953943/70‑1 (MQ=255)
  agaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga                                                                      <  1:35170/70‑1 (MQ=255)
   gaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAg                                                                    <  1:3688200/71‑1 (MQ=255)
     agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg                                                                   >  1:3544019/1‑70 (MQ=255)
      gagGAAATGCAGAAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTa                                                                 >  1:468475/1‑71 (MQ=255)
       agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAt                                                                <  1:148946/71‑1 (MQ=255)
       agGAAATGCAGCAGAATTATTATCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAt                                                                >  1:2042261/1‑71 (MQ=255)
       agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt                                                                >  1:2218775/1‑71 (MQ=255)
        ggAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTAtt                                                               >  1:590098/1‑71 (MQ=255)
          aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATaa                                                                     >  1:3004111/1‑63 (MQ=255)
          aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGt                                                                  <  1:2092416/66‑1 (MQ=255)
           aatgcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                                              >  1:1713405/4‑69 (MQ=255)
              gcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGc                                                         <  1:357132/71‑1 (MQ=255)
                  cagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAggg                                                     >  1:3282911/1‑71 (MQ=255)
                   agAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAtt                                                               <  1:2851185/60‑1 (MQ=255)
                    gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGt                                                    >  1:3750623/1‑70 (MQ=255)
                     aaTTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGtt                                                   >  1:106149/1‑70 (MQ=255)
                     aaTTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTa                                                  <  1:2946905/71‑1 (MQ=255)
                       ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga                                                                      <  1:2267441/49‑1 (MQ=255)
                       ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAg                                                <  1:2251150/71‑1 (MQ=255)
                         attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGAc                                              >  1:2504320/1‑71 (MQ=255)
                         attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGAc                                              >  1:3921119/1‑71 (MQ=255)
                          ttttCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGAc                                              <  1:1024937/70‑1 (MQ=255)
                          ttttCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGAc                                              <  1:3688807/70‑1 (MQ=255)
                            ttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                                            <  1:1830179/70‑1 (MQ=255)
                            ttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGt                                            <  1:3534068/70‑1 (MQ=255)
                             tCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGt                                            >  1:814337/1‑69 (MQ=255)
                             tCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTa                                           <  1:709753/70‑1 (MQ=255)
                               aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCa                                                             <  1:346821/50‑1 (MQ=255)
                                gggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCa                                                        <  1:1375320/54‑1 (MQ=255)
                                 ggATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAg                                                >  1:1663374/1‑61 (MQ=255)
                                 ggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGc                                                         <  1:1538305/52‑1 (MQ=255)
                                  gATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGc                                     <  1:3494260/71‑1 (MQ=255)
                                   aTACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCg                                    <  1:2167799/71‑1 (MQ=255)
                                     aCGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGa                                   >  1:62942/1‑70 (MQ=255)
                                      cGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa                                  <  1:2105306/70‑1 (MQ=255)
                                      cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaaa                                 <  1:1962751/71‑1 (MQ=255)
                                       ggCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaaaa                                <  1:818590/71‑1 (MQ=255)
                                        gCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag                               <  1:2000860/71‑1 (MQ=255)
                                        gCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag                               <  1:3701228/71‑1 (MQ=255)
                                         cacaTTGCTATTATTATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag                               >  1:91430/1‑70 (MQ=255)
                                         cacaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagc                              <  1:2067929/71‑1 (MQ=255)
                                           caTCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTaa                                          <  1:1048425/57‑1 (MQ=255)
                                           caTCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTaa                                          <  1:1898900/57‑1 (MQ=255)
                                            aTTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACg                                             <  1:3178609/53‑1 (MQ=255)
                                             ttGCTATTAATATTTTGACGAATGATAAGGTATTCAAgg                                                          <  1:1036586/39‑1 (MQ=255)
                                             ttGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGcc                          <  1:1063705/71‑1 (MQ=255)
                                             ttGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGcc                          <  1:1390508/71‑1 (MQ=255)
                                             ttGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACg                                             >  1:3125694/1‑52 (MQ=255)
                                             tcgCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGa                                        >  1:3359388/3‑57 (MQ=255)
                                              tGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCa                         >  1:3720397/1‑71 (MQ=255)
                                              tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGAt                                       >  1:845740/1‑57 (MQ=255)
                                               gCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATg                                      >  1:3278508/1‑57 (MQ=255)
                                               gCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCa                         <  1:552225/70‑1 (MQ=255)
                                                cTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGAt                                       <  1:2649186/55‑1 (MQ=255)
                                                cTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATg                       >  1:1260558/1‑71 (MQ=255)
                                                cTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATg                       >  1:3834812/1‑71 (MQ=255)
                                                 tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATg                       <  1:3537645/70‑1 (MQ=255)
                                                  aTTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGc                                     >  1:932122/1‑55 (MQ=255)
                                                  aTTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGa                     >  1:130959/1‑71 (MQ=255)
                                                  aTTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTATGATGCGAAAAGCAGCCATGGa                     >  1:2790156/1‑71 (MQ=255)
                                                   ttAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGa                     <  1:966734/70‑1 (MQ=255)
                                                   ttAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTaa                                          <  1:270869/49‑1 (MQ=255)
                                                    tAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                                            >  1:755133/1‑46 (MQ=255)
                                                    tAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                                            >  1:1340150/1‑46 (MQ=255)
                                                     aaTATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGa                     >  1:3335187/1‑68 (MQ=255)
                                                     aaTATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagca                             >  1:3294999/1‑60 (MQ=255)
                                                     aaTATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGa                     >  1:1241869/1‑68 (MQ=255)
                                                     aaTATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGa                     >  1:1615305/1‑68 (MQ=255)
                                                     aaTATTTTAAAGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagca                             >  1:3251252/1‑60 (MQ=255)
                                                      atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGa                                        <  1:3094054/48‑1 (MQ=255)
                                                      atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                  <  1:1122330/70‑1 (MQ=255)
                                                      atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                  >  1:3891349/1‑70 (MQ=255)
                                                      atatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag                               >  1:2011799/1‑57 (MQ=255)
                                                      atatTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCa                         <  1:21122/63‑1 (MQ=255)
                                                       tatTTTGACGAATGATAAGGTATTCAAGGCAGGGtt                                                   <  1:1273849/36‑1 (MQ=255)
                                                       tatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagcag                            <  1:2992649/59‑1 (MQ=255)
                                                       tatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag                               <  1:2666338/56‑1 (MQ=255)
                                                       tatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACa                   <  1:3529269/68‑1 (MQ=255)
                                                       tatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                 <  1:3083070/70‑1 (MQ=255)
                                                       tatTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGaa                <  1:1368838/71‑1 (MQ=255)
                                                        atTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATg                                      >  1:2987773/1‑48 (MQ=255)
                                                         ttttGACGAATGATAAGGTATTCAAGGCTGGGTTAAGACGTAAGATGCGaa                                  <  1:2653340/51‑1 (MQ=255)
                                                         ttttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa                                  <  1:3496166/51‑1 (MQ=255)
                                                         ttttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGcc                          <  1:1151720/59‑1 (MQ=255)
                                                         ttttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc              >  1:1118617/1‑71 (MQ=255)
                                                          tttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                                            <  1:1737905/40‑1 (MQ=255)
                                                           ttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACt             <  1:2634741/70‑1 (MQ=255)
                                                           ttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTa            >  1:2226934/1‑71 (MQ=255)
                                                           ttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTa            >  1:3356191/1‑71 (MQ=255)
                                                            taaCGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaaa                                 <  1:623721/47‑1 (MQ=255)
                                                             gACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAc           <  1:1429119/70‑1 (MQ=255)
                                                              aCGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCt         >  1:642518/1‑71 (MQ=255)
                                                              aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCt         <  1:486881/71‑1 (MQ=255)
                                                                gAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATg                       <  1:2315423/55‑1 (MQ=255)
                                                                gAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTg        >  1:1157819/1‑70 (MQ=255)
                                                                   tGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATgg                      <  1:2343478/53‑1 (MQ=255)
                                                                   tGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTat           >  1:3368471/1‑63 (MQ=255)
                                                                     agaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCa  <  1:2567343/71‑1 (MQ=255)
                                                                      |                                                                     
TAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCA  >  minE/218036‑218175

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: