Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 865,658 A→G 9.1% F377L (TTT→CTT)  moeA ← molybdopterin biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb865,6580AG 9.1% 144.1 / 6.0 55F377L (TTT→CTT) moeAmolybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base G (2/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGC  >  W3110S.gb/865590‑865723
                                                                    |                                                                 
cACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAgg                                                                 <  1:507189/71‑1 (MQ=255)
    cATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGc                                                              >  1:2222525/1‑70 (MQ=255)
                tCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTaaa                                                                   >  1:321346/1‑53 (MQ=255)
                tCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTaaa                                                                   >  1:1453106/1‑53 (MQ=255)
                tCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTaaa                                                                   >  1:1594332/1‑53 (MQ=255)
                tCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTaaa                                                                   >  1:2291841/1‑53 (MQ=255)
                     aTTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCt                                            <  1:2073135/71‑1 (MQ=255)
                     aTTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCt                                            <  1:3855574/71‑1 (MQ=255)
                      ttGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCt                                            <  1:547734/70‑1 (MQ=255)
                      ttGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCt                                            >  1:462370/1‑70 (MQ=255)
                      ttGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTg                                           <  1:1014588/71‑1 (MQ=255)
                      ttGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTg                                           <  1:251226/71‑1 (MQ=255)
                      ttGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTg                                           <  1:1421628/71‑1 (MQ=255)
                      ttGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTg                                           <  1:181190/71‑1 (MQ=255)
                        gCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:2482083/1‑71 (MQ=255)
                        gCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:1713856/1‑71 (MQ=255)
                        gCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:1505966/1‑71 (MQ=255)
                        gCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:54945/1‑71 (MQ=255)
                         ccGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         <  1:2459220/70‑1 (MQ=255)
                         ccGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         <  1:2207928/70‑1 (MQ=255)
                               tCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTaaa                                                          <  1:1040908/47‑1 (MQ=255)
                                cgcgTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGa                                          >  1:3703449/1‑62 (MQ=255)
                                cgcgTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGa                                          >  1:3926632/1‑62 (MQ=255)
                                cgcgTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCggtg                                 >  1:2410678/1‑71 (MQ=255)
                                   gTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtc                              >  1:3845793/1‑71 (MQ=255)
                                   gTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtc                              >  1:121275/1‑71 (MQ=255)
                                      ccAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGa                                          >  1:321738/1‑56 (MQ=255)
                                      ccAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGa                                          >  1:851251/1‑56 (MQ=255)
                                      ccAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCgg                                   >  1:2915526/1‑63 (MQ=255)
                                      ccAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCgg                                   >  1:1619427/1‑63 (MQ=255)
                                          cACGATAAAGCAGTTTCCGAGGCTAAGGGAGCTAAATATATGTGAACCCTGa                                          >  1:3442671/1‑52 (MQ=255)
                                          cACGATAAAGCAGTTGCCGAGGCTAAGGGAGCTAAATATATGTGAACCCTGa                                          >  1:3727495/1‑52 (MQ=255)
                                           aCGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTc                      >  1:3072214/1‑71 (MQ=255)
                                           aCGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTc                      >  1:3808173/1‑71 (MQ=255)
                                             gATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                        <  1:1044738/35‑1 (MQ=255)
                                             gATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                        <  1:2593012/35‑1 (MQ=255)
                                             gATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                        <  1:1081504/35‑1 (MQ=255)
                                             gATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                        <  1:1175081/35‑1 (MQ=255)
                                             gATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                        <  1:3886710/35‑1 (MQ=255)
                                               tAAAGCAGTTGCCGAGGCTAAAGGAGCTAAatata                                                      >  1:1244498/1‑35 (MQ=255)
                                               tAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGa                                                 >  1:3331415/1‑40 (MQ=255)
                                               tAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGt                  >  1:3608480/1‑71 (MQ=255)
                                                aaaGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:3536956/1‑47 (MQ=255)
                                                aaaGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:2933528/1‑47 (MQ=255)
                                                aaaGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGAt                                         >  1:4072177/1‑47 (MQ=255)
                                                    cAGTTGCCGAGGCTAAAGTAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGc              <  1:20944/70‑1 (MQ=255)
                                                    cAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGc              <  1:3404300/70‑1 (MQ=255)
                                                    cAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGc              <  1:2937705/70‑1 (MQ=255)
                                                         gCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCgg        >  1:1217962/1‑71 (MQ=255)
                                                         gCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCg         >  1:370559/1‑70 (MQ=255)
                                                          ccGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCgg        <  1:1847112/70‑1 (MQ=255)
                                                          ccGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgc       <  1:2762562/71‑1 (MQ=255)
                                                          ccGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgc       <  1:1638360/71‑1 (MQ=255)
                                                          ccGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgc       <  1:2247896/71‑1 (MQ=255)
                                                          ccGAGGCTAAAGGAGCTAAATATATGTGAACCATGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCgg        <  1:2548694/70‑1 (MQ=255)
                                                               gCTAAGGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgtt    <  1:1254375/69‑1 (MQ=255)
                                                               gCTAAGGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgtt    <  1:3878806/69‑1 (MQ=255)
                                                               gCTAAGGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgtt    <  1:3238468/69‑1 (MQ=255)
                                                               gCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgttgc  <  1:2666168/71‑1 (MQ=255)
                                                                   aaGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCAc                         >  1:2362986/1‑44 (MQ=255)
                                                                    |                                                                 
CACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGC  >  W3110S.gb/865590‑865723

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: