Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 938,015 A→G 14.1% Q75Q (CAA→CAG lolA → chaperone for lipoproteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb938,0150AG14.1% 173.2 / 21.3 71Q75Q (CAA→CAGlolAchaperone for lipoproteins
Reads supporting (aligned to +/- strand):  ref base A (35/26);  new base G (6/4);  total (41/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 = 8.47e-01

GAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGC  >  W3110S.gb/937956‑938083
                                                           |                                                                    
gAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAGCCTGATGaa                                                             <  1:177515/69‑1 (MQ=255)
   ggTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGc                                                          >  1:3635131/1‑69 (MQ=255)
    gTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCa                                                         >  1:2894159/1‑69 (MQ=255)
      cAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCt                                                                   <  1:2566652/57‑1 (MQ=255)
      cAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCAtt                                                       >  1:2536798/1‑69 (MQ=255)
          gCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCt                                                                   >  1:3356729/1‑53 (MQ=255)
            gATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGtt                                                 <  1:2177785/69‑1 (MQ=255)
            gATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGtt                                                 <  1:2789028/69‑1 (MQ=255)
               cTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              >  1:2906394/1‑69 (MQ=255)
                    ggTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              <  1:3618835/64‑1 (MQ=255)
                    ggTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              <  1:1187017/64‑1 (MQ=255)
                     gTGAAACGTCCAAACTTATTCAACTTGCATATGACACAACCtgatga                                                              <  1:2851702/47‑1 (MQ=255)
                     gTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCtgatga                                                              >  1:1084166/1‑47 (MQ=255)
                     gTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCtgatga                                                              <  1:2751883/47‑1 (MQ=255)
                     gTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCtgatga                                                              <  1:3721466/47‑1 (MQ=255)
                     gTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCtgatga                                                              >  1:623947/1‑47 (MQ=255)
                     gTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCtgatga                                                              >  1:832424/1‑47 (MQ=255)
                      tGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              <  1:82189/62‑1 (MQ=255)
                      tGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              <  1:1237424/62‑1 (MQ=255)
                      tGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGt                                        <  1:2211602/68‑1 (MQ=255)
                       gAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTaa                                      >  1:561354/1‑69 (MQ=255)
                        aaaCGTCCAAACTTATTCAACTGGCATATGACACAACCt                                                                   <  1:3561435/39‑1 (MQ=255)
                         aaCGTCCAAACTTATTCAACTGGCATATGACACAACCtgatg                                                               >  1:163680/1‑42 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCt                                                     <  1:842228/51‑1 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCt                                                     <  1:1893743/51‑1 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCt                                                     <  1:2615390/51‑1 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTaaaa                                    <  1:3046947/68‑1 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTaaaa                                    <  1:3535698/68‑1 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAac                                   >  1:3662174/1‑69 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAac                                   >  1:1478397/1‑69 (MQ=255)
                          aCGTCCAAACTTATTCAACTGGCATATGAAACAACCTGATGAAAGCATTCt                                                     <  1:1704748/51‑1 (MQ=255)
                           cGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAaca                                  <  1:1671051/69‑1 (MQ=255)
                            gTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAacac                                 >  1:3099859/1‑69 (MQ=255)
                            gTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAacac                                 >  1:2826764/1‑69 (MQ=255)
                             tCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACt                                >  1:2452314/1‑69 (MQ=255)
                             tCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACt                                >  1:1767185/1‑69 (MQ=255)
                             tCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACt                                >  1:1800737/1‑69 (MQ=255)
                               cAAACTTATTCAACTGGCATATGACACAGCCTGATGaa                                                             <  1:3189228/38‑1 (MQ=255)
                               cAAACTTATTCAACTGGCATATGACACAGCCTGATGaa                                                             <  1:2937441/38‑1 (MQ=255)
                               cAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTaa                                      >  1:1342136/1‑61 (MQ=255)
                               cAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTaa                                      >  1:2954193/1‑61 (MQ=255)
                                  aCTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              <  1:1936340/50‑1 (MQ=255)
                                  aCTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCt                                              <  1:2007996/50‑1 (MQ=255)
                                   cTTATTGAACTGGCATATGACACAACCTGATGAAAg                                                           >  1:476590/1‑36 (MQ=255)
                                   cTTATTCAACTGGCATATGACACAACCTTATGAAAg                                                           >  1:1233060/1‑36 (MQ=255)
                                   cTTATTCAACTGGCATATGACACAACCTGATGAAAg                                                           >  1:1381639/1‑36 (MQ=255)
                                   cTTATTCAACTGGCATATGACACAACCTGATGAAAg                                                           >  1:3538862/1‑36 (MQ=255)
                                   cTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTATAACACTGTGGtt                          <  1:2124961/69‑1 (MQ=255)
                                   cTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGtt                          <  1:3135860/69‑1 (MQ=255)
                                   cTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGtt                          <  1:3719779/69‑1 (MQ=255)
                                     tATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTc                         >  1:29471/1‑68 (MQ=255)
                                     tATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTc                         >  1:3517046/1‑68 (MQ=255)
                                         cAACTGGCATATGACACAGCCTGATGAAAGCATTCTTGTTTCTGACGGTAAAaca                                  >  1:3257937/1‑55 (MQ=255)
                                         cAACTGGCATATGACACAGCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAaca                                  >  1:1871493/1‑55 (MQ=255)
                                         cAACTGGCATATGACACAGCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAaca                                  >  1:2580754/1‑55 (MQ=255)
                                         cAACTGGCATATGACACAGCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAaca                                  >  1:1201332/1‑55 (MQ=255)
                                         cAACTGGCATATGACACAGCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAaca                                  >  1:1469400/1‑55 (MQ=255)
                                         cAACTGGCATATGACACAGCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATaa                    <  1:3693039/69‑1 (MQ=255)
                                         cAACTGGCATATGACACAAGCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCta                       <  1:671194/66‑1 (MQ=255)
                                         cAACAGGCATATGACACAGCCTGATGAAAGCATTCt                                                     >  1:1894531/1‑36 (MQ=37)
                                                  tatGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCcgttcg            >  1:3392650/1‑68 (MQ=255)
                                                   atGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACtg                               <  1:3642443/48‑1 (MQ=255)
                                                   atGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCcgttcg            >  1:3267114/1‑67 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:2623635/1‑69 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:2518316/1‑69 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:3656027/1‑69 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:1997642/1‑69 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:3613372/1‑69 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:2741737/1‑69 (MQ=255)
                                                    tGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:2082342/1‑69 (MQ=255)
                                                     gACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTg         >  1:342328/1‑68 (MQ=255)
                                                          aaCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCa     >  1:3413076/1‑67 (MQ=255)
                                                          aaCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCa     >  1:1247949/1‑67 (MQ=255)
                                                          aaCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCa     >  1:1045748/1‑67 (MQ=255)
                                                           aCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGc  <  1:2451921/69‑1 (MQ=255)
                                                           aCCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGc  <  1:1687330/69‑1 (MQ=255)
                                                           |                                                                    
GAAGGTCAGGGCGATCTGTGGGTGAAACGTCCAAACTTATTCAACTGGCATATGACACAACCTGATGAAAGCATTCTGGTTTCTGACGGTAAAACACTGTGGTTCTATAACCCGTTCGTTGAGCAAGC  >  W3110S.gb/937956‑938083

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: