Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,416,062 A→G 100% N27D (AAT→GAT)  yjfK → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,416,0620AG91.5% 142.6 / 3.9 47N27D (AAT→GAT) yjfKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (17/26);  total (18/29)
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.67e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTG  >  W3110S.gb/4416002‑4416127
                                                            |                                                                 
cGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCtt                                                         <  1:1036333/71‑1 (MQ=255)
 gTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCtt                                                         <  1:329698/70‑1 (MQ=255)
               gATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGAt                                                                 <  1:1259083/48‑1 (MQ=255)
                    taaGCCGGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGc                                       >  1:1909189/1‑69 (MQ=255)
                    taaGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgt                                     >  1:442824/1‑71 (MQ=255)
                    taaGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgt                                     >  1:1604491/1‑71 (MQ=255)
                    taaGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgt                                     >  1:677609/1‑71 (MQ=255)
                    taaGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcg                                      >  1:654539/1‑70 (MQ=255)
                       gCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGAAACGCTAAcgtgtc                                  <  1:2310285/71‑1 (MQ=255)
                             tATCGCTCGTTGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgtttc                                  <  1:1451146/65‑1 (MQ=255)
                             tATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgtttc                                  <  1:84709/65‑1 (MQ=255)
                             tATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgtttc                                  <  1:1980553/65‑1 (MQ=255)
                               tCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCttt                                                        <  1:683852/41‑1 (MQ=255)
                                cGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgtttc                                  <  1:1832085/62‑1 (MQ=255)
                                 gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTAc                                                      >  1:1278518/1‑41 (MQ=255)
                                 gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTAc                                                      >  1:512068/1‑41 (MQ=255)
                                 gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTg                             <  1:1484774/66‑1 (MQ=255)
                                 gCTCGTGGGCCGCTGGGACTTCATCTCGATAGTGGCTTTAc                                                      >  1:2436449/1‑41 (MQ=255)
                                   tCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGc                            >  1:44723/1‑65 (MQ=255)
                                   tCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCt                           >  1:767448/1‑66 (MQ=255)
                                   tCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCt                           >  1:2074738/1‑66 (MQ=255)
                                     gTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGaa                       <  1:2379987/68‑1 (MQ=255)
                                      tGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                   <  1:2118348/71‑1 (MQ=255)
                                      tGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                   <  1:588294/71‑1 (MQ=255)
                                      tGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                   <  1:1220864/71‑1 (MQ=255)
                                      tGGTCCGCCGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                   <  1:778450/71‑1 (MQ=255)
                                       ggTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGa                     >  1:546063/1‑68 (MQ=255)
                                       ggTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                   >  1:369114/1‑70 (MQ=255)
                                        gTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                 <  1:1714017/71‑1 (MQ=255)
                                        gTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                 <  1:136636/71‑1 (MQ=255)
                                        gTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                 <  1:44386/71‑1 (MQ=255)
                                           cGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGt                                           >  1:2222859/1‑42 (MQ=255)
                                           cGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGt                                           >  1:2640052/1‑42 (MQ=255)
                                           cGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGt                                           >  1:1752487/1‑42 (MQ=255)
                                           cGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGt                                           >  1:1333854/1‑42 (MQ=255)
                                             cTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt            <  1:273254/71‑1 (MQ=255)
                                             cTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt            <  1:1997678/71‑1 (MQ=255)
                                             cTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt            <  1:687688/71‑1 (MQ=255)
                                              tGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTa                                         <  1:156571/41‑1 (MQ=255)
                                                ggACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGGAGATGAATTGCTg           <  1:604993/69‑1 (MQ=255)
                                                ggACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTg           <  1:1001888/69‑1 (MQ=255)
                                                 gACTTCATCTCGATAGTGTCTTTACGCTCGACACGTTAGcgtttc                                  <  1:90934/45‑1 (MQ=255)
                                                    ttCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgtttc                                  <  1:1777325/42‑1 (MQ=255)
                                                    ttCATCTCGATAGTGGCTTTACGCTCGACACGTTAGcgtttc                                  <  1:1277842/42‑1 (MQ=255)
                                                    ttCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                   <  1:2661918/57‑1 (MQ=255)
                                                       aTCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCt   <  1:2189508/70‑1 (MQ=255)
                                                       aTCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTg  >  1:2127156/1‑71 (MQ=255)
                                                            |                                                                 
CGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTG  >  W3110S.gb/4416002‑4416127

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: