Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 771,768 C→T 15.5% intergenic (+735/‑112) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb771,7680CT15.5% 147.3 / 17.4 58intergenic (+735/‑112)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (4/5);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

TATATACGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGTGTGACCCTTTCTAAC  >  W3110S.gb/771699‑771835
                                                                     |                                                                   
tatataCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCa                                                                    <  1:4169359/71‑1 (MQ=255)
  tataCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACATTTCTTTACGCCGTTTTTGTGTGCAt                                                                   >  1:2996796/1‑70 (MQ=255)
   ataCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCAtt                                                                  >  1:2540883/1‑70 (MQ=255)
    taCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTc                                                                 <  1:587331/70‑1 (MQ=255)
    taCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTc                                                                 <  1:1532349/70‑1 (MQ=255)
     aCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCa                                                                    >  1:460000/1‑66 (MQ=255)
     aCGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTca                                                                >  1:305889/1‑70 (MQ=255)
          ggATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTcaca                                                              <  1:920116/67‑1 (MQ=255)
           gATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCATTACGCCGTTTTTGTGTGCATTAACAt                                                             <  1:2254795/67‑1 (MQ=255)
              cAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCa                                                                    >  1:3957877/1‑57 (MQ=255)
              cAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTatga                                                      >  1:41017/1‑71 (MQ=255)
                aTTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGTATTCACATTGTatgatg                                                    <  1:3703340/71‑1 (MQ=255)
                 tttGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCAGTTTTTGTGTGCATTCACATGGTatgatg                                                    <  1:4516287/70‑1 (MQ=255)
                   tGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTatgatg                                                    >  1:1469959/1‑68 (MQ=255)
                       tttttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCa                                                                    <  1:762251/48‑1 (MQ=255)
                         tttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCAt                                                                   >  1:2878452/1‑47 (MQ=255)
                         tttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTcaca                                                              >  1:2692625/1‑52 (MQ=255)
                         tttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTatgat                                                     >  1:1647461/1‑61 (MQ=255)
                         tttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTatga                                                      >  1:2896155/1‑60 (MQ=255)
                         tttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGtt                                            >  1:1860158/1‑70 (MQ=255)
                         tttATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGtt                                            >  1:1291201/1‑70 (MQ=255)
                           tattaACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGTATTCACATGGTATGATGAAAGTGTTca                                          >  1:2385800/1‑70 (MQ=255)
                             ttaACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGTATTCACATGGTATGATGAAAGTGTTcaaac                                       <  1:3138458/71‑1 (MQ=255)
                                aCATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTAtgatga                                                   <  1:4356339/56‑1 (MQ=255)
                                 cATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAAt                                   <  1:882404/71‑1 (MQ=255)
                                 cATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTCCATTCACATGGTATGATGAAAg                                                >  1:2017331/1‑58 (MQ=255)
                                  aTGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTatgatg                                                    <  1:2438291/53‑1 (MQ=255)
                                   tGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGTATTCACATGGTATGATGaa                                                  <  1:3965964/54‑1 (MQ=255)
                                   tGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCa                                                                    <  1:76991/36‑1 (MQ=255)
                                   tGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGaaa                                                 >  1:2805230/1‑55 (MQ=255)
                                     tttGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTc                                <  1:114953/70‑1 (MQ=255)
                                     tttGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCt                               <  1:2048198/71‑1 (MQ=255)
                                        gCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTAtt                            >  1:3050757/1‑71 (MQ=255)
                                        gCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTAtt                            >  1:3045176/1‑71 (MQ=255)
                                        gCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTAt                             <  1:2639940/70‑1 (MQ=255)
                                            ccTTTCTTTACGCCGTTTTTGTGTGTATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTgggg                        <  1:1442619/71‑1 (MQ=255)
                                             ctttctttACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAg                                                >  1:2286696/1‑46 (MQ=255)
                                             ctttctttACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTcaaacaa                                     <  1:344425/57‑1 (MQ=255)
                                              tttctttACGCCGTTTTTGTGTGTATTCACATGGTatga                                                      >  1:731875/1‑39 (MQ=255)
                                              tttctttACGCCGTTTTTGTGTGTATTCACATGGTatga                                                      >  1:679196/1‑39 (MQ=255)
                                              tttctttACGCCGTTTTTGTGTGTATTCACATGGTatga                                                      >  1:401647/1‑39 (MQ=255)
                                              tttctttACGCCGTTTTTGTGTGTATTCACATGGTatg                                                       <  1:1209678/38‑1 (MQ=255)
                                              tttctttACGCCGTTTTTGTGTGCATTCACATGGTAtgatga                                                   <  1:2252208/42‑1 (MQ=255)
                                               ttctttACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAAt                                   >  1:1176435/1‑57 (MQ=255)
                                                tctttACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTAtt                            >  1:1812408/1‑63 (MQ=255)
                                                 ctttACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTgg                          >  1:2743081/1‑64 (MQ=255)
                                                  tttACGACGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTAtt                            <  1:544219/61‑1 (MQ=255)
                                                      cGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCgtgt               <  1:2628694/70‑1 (MQ=255)
                                                      cGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCgtgt               <  1:3532469/70‑1 (MQ=255)
                                                      cGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGtgtg              <  1:1240377/71‑1 (MQ=255)
                                                             tttGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTAtt                            >  1:4449847/1‑50 (MQ=255)
                                                             tttGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTgggg                        >  1:1864102/1‑54 (MQ=255)
                                                             tttGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCgt                 >  1:610138/1‑61 (MQ=255)
                                                                gtgtgCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGTGTGACCCTTTCt     <  1:3141140/70‑1 (MQ=255)
                                                                 tgtgCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGTGTGACCCTTTCTa    <  1:4314068/70‑1 (MQ=255)
                                                                   tgCATTCACATGGTATGATTAAAGTGTTCAAACAAATTTCTATTGGGGCATGCgtgt               <  1:1812208/57‑1 (MQ=255)
                                                                   tgCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGTGTGACCCTTTCt     <  1:1638397/67‑1 (MQ=255)
                                                                   tgCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGTGTGACCCTTTCTAAc  <  1:424325/70‑1 (MQ=255)
                                                                     |                                                                   
TATATACGTGGGATCAATTTGAGTTTTTATTAACATGTTTGCAACCTTTCTTTACGCCGTTTTTGTGTGCATTCACATGGTATGATGAAAGTGTTCAAACAAATTTCTATTGGGGCATGCGTGTGACCCTTTCTAAC  >  W3110S.gb/771699‑771835

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: