Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,416,062 A→G 33.3% 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,0620AG33.3% 72.2 / 57.7 63N27D (AAT→GAT) yjfKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/18);  new base G (12/9);  total (36/27)
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

TTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGG  >  W3110S.gb/4415993‑4416131
                                                                     |                                                                     
tttttCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCaa                                                                      >  1:3139041/1‑71 (MQ=255)
  tttCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAAt                                                                     <  1:2213280/70‑1 (MQ=255)
       aGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGc                                                               >  1:3658598/1‑71 (MQ=255)
             tGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACg                                                         >  1:1106385/1‑71 (MQ=255)
                ttGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCaa                                                                      <  1:2876175/55‑1 (MQ=255)
                ttGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATa                                                                    >  1:4611170/1‑57 (MQ=255)
                ttGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTc                                                      >  1:532448/1‑71 (MQ=255)
                 tGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTc                                                      <  1:3532129/70‑1 (MQ=255)
                 tGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTc                                                      <  1:385210/70‑1 (MQ=255)
                 tGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCT‑‑GACTTCATCTCAATAGTGGCTTTACGCTc                                                      <  1:2072243/68‑1 (MQ=255)
                  ggCAAGGATAATAAGCCAGCTAACGCTCGTGGTCCGCTGGGACTTCATCTCAATa                                                                    >  1:1755948/1‑55 (MQ=255)
                        gATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGt                                               >  1:1308165/1‑70 (MQ=255)
                            ataaGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcg                                          >  1:4337776/1‑71 (MQ=255)
                               aGCCACCTATCGCTCGTGGTCCGCTGCGACTTCATCTCGAt                                                                     >  1:3008215/1‑41 (MQ=255)
                                 ccAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGa                                                    <  1:4490239/56‑1 (MQ=255)
                                   aGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgtt                                   <  1:3603696/71‑1 (MQ=255)
                                   aGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgtt                                   >  1:465649/1‑71 (MQ=255)
                                   aGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgtt                                   <  1:787052/71‑1 (MQ=255)
                                      tATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTga                                >  1:4566419/1‑70 (MQ=255)
                                      tATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGc                                >  1:2907536/1‑71 (MQ=255)
                                        tCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTg                              >  1:1860799/1‑71 (MQ=255)
                                        tCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTg                              <  1:3399339/71‑1 (MQ=255)
                                        tCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCtt                                                             <  1:364490/40‑1 (MQ=255)
                                        tCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCtt                                                             <  1:276066/40‑1 (MQ=255)
                                        tCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgtt                                   <  1:4444409/66‑1 (MQ=255)
                                        tCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGc                                <  1:522217/69‑1 (MQ=255)
                                         cGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTg                                 <  1:3014097/67‑1 (MQ=255)
                                         cGCTCGTGGTCCGCTGGGACTTAATCTCAATAGTGGCt                                                              >  1:305616/1‑38 (MQ=255)
                                          gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTa                                                           >  1:2573220/1‑40 (MQ=255)
                                          gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGac                                                   >  1:4719883/1‑48 (MQ=255)
                                          gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGac                                                   >  1:1367960/1‑48 (MQ=255)
                                          gCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTgg                             >  1:816763/1‑70 (MQ=255)
                                           cTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCt                               >  1:154607/1‑67 (MQ=255)
                                            tCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTgg                             >  1:519214/1‑68 (MQ=255)
                                            tCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTgg                             >  1:513465/1‑68 (MQ=255)
                                               tGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                       <  1:4300747/71‑1 (MQ=255)
                                               tGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                       <  1:4453933/71‑1 (MQ=255)
                                               tGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                       <  1:1108198/71‑1 (MQ=255)
                                                ggTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                       >  1:998395/1‑70 (MQ=255)
                                                   ccGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGc                                >  1:4063250/1‑58 (MQ=255)
                                                    cGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTg                  >  1:3137239/1‑71 (MQ=255)
                                                    cGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcg                                     <  1:2369201/52‑1 (MQ=255)
                                                     gCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGTTGAATTGc                 <  1:4242138/71‑1 (MQ=255)
                                                      cTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt                <  1:1333747/71‑1 (MQ=255)
                                                        gggACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt                >  1:2592679/1‑69 (MQ=255)
                                                        gggACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGa              >  1:4656750/1‑71 (MQ=255)
                                                           aCTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                     >  1:3116860/1‑61 (MQ=255)
                                                           aCTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                     >  1:659282/1‑61 (MQ=255)
                                                           aCTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                     >  1:1565081/1‑61 (MQ=255)
                                                           aCTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTg           <  1:4727082/71‑1 (MQ=255)
                                                            cTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTgc          <  1:1519976/71‑1 (MQ=255)
                                                             ttCATGTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTg                  <  1:1897592/62‑1 (MQ=255)
                                                             ttCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                       <  1:2032833/57‑1 (MQ=255)
                                                             ttCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTg                  <  1:4799426/62‑1 (MQ=255)
                                                             ttCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTg                  <  1:2068531/62‑1 (MQ=255)
                                                             ttCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgtt                                   >  1:4009317/1‑45 (MQ=255)
                                                             ttCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgtt                                   >  1:3561575/1‑45 (MQ=255)
                                                             ttCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAATATg                       <  1:2810246/57‑1 (MQ=255)
                                                              tCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTgcgc        >  1:2747213/1‑71 (MQ=255)
                                                                aTCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTg      >  1:1336530/1‑71 (MQ=255)
                                                                 tctcGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTgc          <  1:3101992/66‑1 (MQ=255)
                                                                   tcAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgt                                    >  1:4025090/1‑38 (MQ=255)
                                                                   tcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                     >  1:1555366/1‑53 (MQ=255)
                                                                   tcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCg   >  1:3543271/1‑71 (MQ=255)
                                                                    cAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCgg  >  1:2729569/1‑71 (MQ=255)
                                                                     |                                                                     
TTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGG  >  W3110S.gb/4415993‑4416131

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: