Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,017,194 G→A 12.0% A310T (GCA→ACA)  manA → mannose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,017,1940GA12.0% 155.1 / 10.6 58A310T (GCA→ACA) manAmannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (29/22);  new base A (4/3);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAA  >  minE/1017125‑1017262
                                                                     |                                                                    
cTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGc                                                                     <  1:3434943/71‑1 (MQ=255)
   gTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAg                                                                   <  1:2542374/70‑1 (MQ=255)
        cAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCGGCCGGTGAAACAAGGTGCAGAACTg                                                              >  1:2231681/1‑70 (MQ=255)
        cAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTg                                                              >  1:558381/1‑70 (MQ=255)
        cAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTACAGAAc                                                                >  1:3477218/1‑68 (MQ=255)
          aTGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGAc                                                           >  1:926206/1‑71 (MQ=255)
          aTGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGAc                                                           >  1:234715/1‑71 (MQ=255)
           tgtgAAATTCGAAGCCAAACAGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTACAGAACTGGACt                                                          >  1:2096693/1‑71 (MQ=255)
             tgAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTc                                                        <  1:1171663/71‑1 (MQ=255)
                 aTTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCa                                                                    >  1:1310915/1‑55 (MQ=255)
                 aTTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGa                                                    >  1:3359387/1‑71 (MQ=255)
                  ttCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGAt                                                   >  1:3351997/1‑71 (MQ=255)
                  ttCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGAt                                                   >  1:1971291/1‑71 (MQ=255)
                   tCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGAtt                                                  >  1:3361571/1‑71 (MQ=255)
                     gAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGa                                                    >  1:888656/1‑67 (MQ=255)
                     gAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTc                                                 >  1:3078940/1‑70 (MQ=255)
                         ccAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTcc                                                <  1:3731272/67‑1 (MQ=255)
                         ccAAACCGGCAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTg                                            >  1:3864601/1‑70 (MQ=255)
                           aaaCCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTg                                            <  1:3882903/69‑1 (MQ=255)
                           aaaCCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTACAGAACTGGACTTCCCGATTCCAGTGga                                          <  1:2410041/71‑1 (MQ=255)
                            aaCCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgat                                         >  1:2506117/1‑71 (MQ=255)
                             aCCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgat                                         >  1:1830558/1‑70 (MQ=255)
                             aCCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgat                                         >  1:1953684/1‑70 (MQ=255)
                              ccGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTACAGAACTGGACTTCCCGAtt                                                  >  1:1370330/1‑60 (MQ=255)
                                  cTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtttt                                   <  1:1845297/71‑1 (MQ=255)
                                  cTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtttt                                   <  1:2958817/71‑1 (MQ=255)
                                   tAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAg                                              >  1:3613155/1‑59 (MQ=255)
                                     aCCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGcc                                >  1:824929/1‑71 (MQ=255)
                                        aGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtc                             <  1:1316856/71‑1 (MQ=255)
                                         gttgttGACCCAGCCGGTGAAACAAGGTACAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtct                            >  1:2291322/1‑71 (MQ=255)
                                          ttgttgACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCt                               >  1:1934565/1‑67 (MQ=255)
                                          ttgttgACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCt                               >  1:3651466/1‑67 (MQ=255)
                                          ttgttgACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtctc                           <  1:80998/71‑1 (MQ=255)
                                           tgttgACCCAGCCGGTGAAACAAGGTGCAGAACTg                                                              >  1:1310912/1‑35 (MQ=255)
                                           tgttgACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtt                                     >  1:1214638/1‑60 (MQ=255)
                                            gttgACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCg                          <  1:1807653/70‑1 (MQ=255)
                                               gACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCa                                               <  1:1572234/46‑1 (MQ=255)
                                               gACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTg                       >  1:1027823/1‑70 (MQ=255)
                                                aCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtct                            <  1:2470160/64‑1 (MQ=255)
                                                aCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGc                      <  1:2140267/70‑1 (MQ=255)
                                                    aGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAttt                                    >  1:674183/1‑52 (MQ=255)
                                                    aGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAttt                                    >  1:2938176/1‑52 (MQ=255)
                                                    aGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGAc                 >  1:2253189/1‑71 (MQ=255)
                                                       cGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCtt              <  1:2768126/71‑1 (MQ=255)
                                                        ggTGAAACAAGGTACAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTa             <  1:3161604/71‑1 (MQ=255)
                                                         gTGAAACAAGGTACAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGcc                                <  1:3872962/51‑1 (MQ=255)
                                                          tGAAACAAGGTGCTGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCa                     <  1:3294512/61‑1 (MQ=255)
                                                          tGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGc                         <  1:3063838/57‑1 (MQ=255)
                                                          tGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGt           <  1:1038486/71‑1 (MQ=255)
                                                            aaaCAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCtt                              <  1:3241889/50‑1 (MQ=255)
                                                            aaaCAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGACTTCTCGc                         >  1:196665/1‑55 (MQ=255)
                                                             aaCAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGa         >  1:981975/1‑70 (MQ=255)
                                                              aCAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtt                                     <  1:224164/41‑1 (MQ=255)
                                                               cAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtctc                           <  1:510817/50‑1 (MQ=255)
                                                                aaGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTa             >  1:50366/1‑63 (MQ=255)
                                                                  ggTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAg    >  1:3728274/1‑70 (MQ=255)
                                                                  ggTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGa   <  1:1976539/71‑1 (MQ=255)
                                                                   gTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGaa  <  1:2170019/71‑1 (MQ=255)
                                                                     |                                                                    
CTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAA  >  minE/1017125‑1017262

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: