Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 965,646 C→T 8.9% G83D (GGC→GAC)  ydeM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE965,6460CT 8.9% 315.0 / 14.5 101G83D (GGC→GAC) ydeMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/53);  new base T (3/6);  total (42/59)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTGGGTTCACC  >  minE/965579‑965713
                                                                   |                                                                   
aaaGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTg                                                                  <  1:3624547/71‑1 (MQ=255)
  aGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCa                                                                >  1:3797910/1‑71 (MQ=255)
  aGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCa                                                                >  1:416174/1‑71 (MQ=255)
   ggCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               <  1:4535664/71‑1 (MQ=255)
    gCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATa                                                              <  1:313503/71‑1 (MQ=255)
     cacaCCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAg                                                             <  1:3236791/71‑1 (MQ=255)
      acacCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGc                                                                 <  1:2191605/66‑1 (MQ=255)
      acacCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGc                                                            <  1:3138502/71‑1 (MQ=255)
        acCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCa                                                                >  1:3142641/1‑65 (MQ=255)
        acCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               >  1:4284484/1‑66 (MQ=255)
        acCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               >  1:2224367/1‑66 (MQ=255)
        acCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgt                                                          <  1:3654183/71‑1 (MQ=255)
        acCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgt                                                          <  1:3407338/71‑1 (MQ=255)
         ccATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgtt                                                         >  1:1158337/1‑71 (MQ=255)
          cattcattATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgttg                                                        <  1:2282/71‑1 (MQ=255)
            ttcattATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAg                                                             <  1:3437677/64‑1 (MQ=255)
            ttcattATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAg                                                             <  1:2505698/64‑1 (MQ=255)
              cattATTCAATAAAATGCCATTAGTTTGTAATGCATTAAAAAAACGTTATTGTCCTGc                                                                 >  1:2511746/1‑58 (MQ=38)
                      aaTAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGaa                                             >  1:4339535/1‑70 (MQ=255)
                       aTAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGTATGCATAGCGTTGTTGATAGTGAAt                                            <  1:3310844/70‑1 (MQ=255)
                       aTAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgttg                                                        >  1:278473/1‑58 (MQ=255)
                       aTAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgttg                                                        >  1:4703402/1‑58 (MQ=255)
                        tAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTg                                                                  <  1:3453377/47‑1 (MQ=255)
                        tAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTg                                                                  <  1:4093703/47‑1 (MQ=255)
                        tAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTg                                                                  <  1:218818/47‑1 (MQ=255)
                        tAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATaa                                          >  1:2650890/1‑71 (MQ=255)
                         aaaaTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGACAGTGAATAAc                                         >  1:2320346/1‑71 (MQ=255)
                          aaaTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               >  1:3173023/1‑48 (MQ=255)
                          aaaTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               >  1:4130173/1‑48 (MQ=255)
                          aaaTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               >  1:3934992/1‑48 (MQ=255)
                          aaaTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCAt                                                               >  1:2240723/1‑48 (MQ=255)
                           aaTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACtt                                       <  1:4030711/71‑1 (MQ=255)
                            aTGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACttt                                      <  1:805488/71‑1 (MQ=255)
                             tGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTa                                     <  1:4216762/71‑1 (MQ=255)
                             tGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTa                                     <  1:4043756/71‑1 (MQ=255)
                                cATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACgg                                  <  1:554133/71‑1 (MQ=255)
                                  ttCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAg                                                             <  1:1553363/42‑1 (MQ=255)
                                  ttCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGc                                                            <  1:304187/43‑1 (MQ=255)
                                   tCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATa                                                  <  1:3742102/52‑1 (MQ=255)
                                   tCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATa                                                  <  1:3100820/52‑1 (MQ=255)
                                    cGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCg                                                           >  1:3032696/1‑42 (MQ=255)
                                    cGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCg                                                           >  1:4279856/1‑42 (MQ=255)
                                    cGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGttgttg                                                     >  1:2898406/1‑48 (MQ=255)
                                    cGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaa                               >  1:1756840/1‑70 (MQ=255)
                                    cGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaa                               >  1:3559710/1‑70 (MQ=255)
                                      tttGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaaa                            <  1:2630350/71‑1 (MQ=255)
                                      tttGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaaa                            <  1:346300/71‑1 (MQ=255)
                                      tttGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaaa                            <  1:3818530/71‑1 (MQ=255)
                                        tGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAAt                           <  1:3122803/70‑1 (MQ=255)
                                         gTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATc                          <  1:4757292/70‑1 (MQ=255)
                                          tAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATcc                         <  1:2754532/70‑1 (MQ=255)
                                          tAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATcc                         <  1:3498319/70‑1 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGTCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaa                             <  1:1102307/64‑1 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGTCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaa                             <  1:4053097/64‑1 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGTCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaa                             <  1:2013593/64‑1 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGTCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaa                             <  1:1878530/64‑1 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGTCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaa                             <  1:1660388/64‑1 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAgg                      >  1:4387277/1‑71 (MQ=255)
                                            aTGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAgg                      >  1:1629943/1‑71 (MQ=255)
                                             tGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgc                     <  1:3468674/71‑1 (MQ=255)
                                             tGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgc                     <  1:2359000/71‑1 (MQ=255)
                                             tGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgc                     <  1:1485514/71‑1 (MQ=255)
                                             tGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgc                     <  1:374060/71‑1 (MQ=255)
                                               cATTAAAAATACGTTTTTGGCCTGCATAGCGttgttg                                                     <  1:3160746/37‑1 (MQ=255)
                                               cATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgc                     >  1:583509/1‑69 (MQ=255)
                                                 ttAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTg                                               <  1:3823301/41‑1 (MQ=255)
                                                  tAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:838688/1‑70 (MQ=255)
                                                   aaaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaaaaa                            >  1:2860158/1‑58 (MQ=255)
                                                   aaaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:3191647/1‑69 (MQ=255)
                                                   aaaaaTACGTTTTTGGCCTGCACAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:2322061/1‑69 (MQ=255)
                                                    aaaaTACGTTTTTGGTCTGCATAGCGTTGTTGATa                                                  >  1:922253/1‑35 (MQ=255)
                                                    aaaaTACGTTTTTGGTCTGCATAGCGTTGTTGATa                                                  >  1:3236658/1‑35 (MQ=255)
                                                    aaaaTACGTTTTTGGTCTGCATAGCGTTGTTGATa                                                  >  1:3231282/1‑35 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagcca               <  1:630486/70‑1 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagcca               <  1:2537879/70‑1 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:108122/1‑68 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:3627041/1‑68 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:3621812/1‑68 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:1313779/1‑68 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:151067/1‑68 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGgccagc                 >  1:4628995/1‑68 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGccagccag              <  1:3740763/71‑1 (MQ=255)
                                                    aaaaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGccagccag              <  1:3068130/71‑1 (MQ=255)
                                                      aaTACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCagag            >  1:3216239/1‑71 (MQ=255)
                                                       aTACGTTTTTGGCCTGCATAACGTTGTTGATAGTGAATAACTTTACGGaaaaaa                            <  1:768553/54‑1 (MQ=255)
                                                        tACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTg          <  1:2740340/71‑1 (MQ=255)
                                                         aCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTgg         >  1:327260/1‑71 (MQ=255)
                                                          cGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTggg        <  1:4679781/71‑1 (MQ=255)
                                                           gTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTGGGt       >  1:3825159/1‑71 (MQ=255)
                                                           gTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTGGGt       <  1:4187411/71‑1 (MQ=255)
                                                            tttttGGCCTGCATAGCGTTGTTGATAGTGAATAACttt                                      <  1:2681045/39‑1 (MQ=255)
                                                            tttttGGCCTGCATAGCGTTGTTGATAGTGAATAACttt                                      <  1:1232985/39‑1 (MQ=255)
                                                            tttttGGCCTGCATAGCGTTGTTGATAGTGAATAACttt                                      <  1:1253689/39‑1 (MQ=255)
                                                            tttttGGCCTGCATAGCGTTGTTGATAGTGAATAACttt                                      <  1:688739/39‑1 (MQ=255)
                                                            tttttGGCCTGCATAGCGTTGTTGATAGTGAATAACttt                                      <  1:3177231/39‑1 (MQ=255)
                                                            tttttGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTGGGtt      >  1:4695323/1‑71 (MQ=255)
                                                            tttttGGCCTGCATAGAGTTGTTGATAGTGAATAACttt                                      <  1:1641644/39‑1 (MQ=255)
                                                              tttGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAg                       >  1:1692874/1‑52 (MQ=255)
                                                                 ggCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTGGGTTCAcc  <  1:2425779/70‑1 (MQ=255)
                                                                  gCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaaa                               <  1:4726686/40‑1 (MQ=255)
                                                                   |                                                                   
AAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAAAAAATCCAGGCCAGCCAGAGTGGGTTCACC  >  minE/965579‑965713

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: