Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,747,558 A→G 6.2% K13K (AAA→AAG yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,747,5580AG 6.2% 166.4 / 2.9 64K13K (AAA→AAGyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (2/2);  total (30/34)
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.45e-01

CCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGT  >  W3110S.gb/2747491‑2747626
                                                                   |                                                                    
ccTATTTTTCCGGTTCAGAAACTGGAATGATGACACTCAACCGCTATCGCCTGAGACACATGGCGAAGCa                                                                    <  1:4218376/70‑1 (MQ=255)
ccTATTTTTCCGGTTCAGAAACTGGAATGATGACACTCAACCGCTATCGCCTGAGACACATGGCGAAGCa                                                                    <  1:1625872/70‑1 (MQ=255)
      tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTaa                                                              >  1:413360/1‑70 (MQ=255)
      tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTaa                                                              >  1:2686262/1‑70 (MQ=255)
           gggTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTc                                                        >  1:3509138/1‑71 (MQ=255)
            ggTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCg                                                       >  1:71724/1‑71 (MQ=255)
              tCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGc                                                     >  1:1708477/1‑71 (MQ=255)
                cGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGcc                                                    >  1:4177169/1‑70 (MQ=255)
                cGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGcc                                                    >  1:291878/1‑70 (MQ=255)
                   aaCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAc                                                >  1:552765/1‑71 (MQ=255)
                   aaCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAc                                                >  1:304389/1‑71 (MQ=255)
                      cGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAcgcg                                             <  1:1239252/71‑1 (MQ=255)
                         aaTGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTc                                           >  1:378357/1‑70 (MQ=255)
                          atgatgACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCg                                          >  1:1477265/1‑70 (MQ=255)
                             atgaCCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:2521027/70‑1 (MQ=255)
                                   cTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:648688/64‑1 (MQ=255)
                                   cTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:99429/64‑1 (MQ=255)
                                   cTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgct                                <  1:1013648/71‑1 (MQ=255)
                                   cTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGGCGaaa                                       <  1:1794534/64‑1 (MQ=255)
                                    tCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:292201/63‑1 (MQ=255)
                                     cAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                              >  1:2876169/1‑71 (MQ=255)
                                     cAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                              >  1:482916/1‑71 (MQ=255)
                                      aaCCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgct                                >  1:2477406/1‑68 (MQ=255)
                                      aaCCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCg                             >  1:3132036/1‑71 (MQ=255)
                                       aCCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCa                                                   <  1:2350633/48‑1 (MQ=255)
                                        ccGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGt                            >  1:3637976/1‑70 (MQ=255)
                                        ccGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTa                           <  1:451744/71‑1 (MQ=255)
                                         cGCTATCGTCTGCGACATATGGCGAGACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTaa                          <  1:1517256/71‑1 (MQ=255)
                                         cGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCaaa                                                 >  1:209795/1‑48 (MQ=255)
                                         cGCTATCGTCTGCGACATATGGCGAAAAAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTaa                          <  1:4156207/71‑1 (MQ=255)
                                          gCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:1605973/57‑1 (MQ=255)
                                           cTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                              <  1:1213321/65‑1 (MQ=255)
                                           cTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGc                        >  1:4049286/1‑71 (MQ=255)
                                           cTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGc                        <  1:966426/71‑1 (MQ=255)
                                           cTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGc                        <  1:2019036/71‑1 (MQ=255)
                                              tCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCa                                                   <  1:1422268/41‑1 (MQ=255)
                                               cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:1993000/52‑1 (MQ=255)
                                               cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                    >  1:1647719/1‑71 (MQ=255)
                                               cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                    >  1:2570019/1‑71 (MQ=255)
                                               cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                    >  1:3918556/1‑71 (MQ=255)
                                                gccTGAGACACATGGCGAAGCAAGGTAATCGCTCGGCAAAACGCGTCGAAAAATTACTGCGCAAGCCAGAc                   >  1:3804287/3‑71 (MQ=255)
                                                gccTGAGACACATGGCGAAGCAAGGTAATCGCTCGGCAAAACGCGTCGAAAAATTACTGCGCAAGCCAGAc                   >  1:1102440/3‑71 (MQ=255)
                                                gTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCa                      >  1:2106653/1‑68 (MQ=255)
                                                  cTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCg                 >  1:872863/1‑71 (MQ=255)
                                                   tGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCg                             <  1:488798/58‑1 (MQ=255)
                                                     cGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                    >  1:208327/1‑65 (MQ=255)
                                                     cGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCg                 <  1:1314939/68‑1 (MQ=255)
                                                      gACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                         <  1:3073498/59‑1 (MQ=255)
                                                      gACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                         <  1:1652704/59‑1 (MQ=255)
                                                      gACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                         <  1:2486472/59‑1 (MQ=255)
                                                        cATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                              <  1:2351533/52‑1 (MQ=255)
                                                        cATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAt           <  1:4325109/71‑1 (MQ=255)
                                                         atatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAg                     <  1:4102493/60‑1 (MQ=255)
                                                         atatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAg                     <  1:2112346/60‑1 (MQ=255)
                                                         atatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAt           >  1:2229418/1‑70 (MQ=255)
                                                         atatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAt           >  1:3214558/1‑70 (MQ=255)
                                                          tatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                                       <  1:4368946/41‑1 (MQ=255)
                                                          tatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGGCGaaaaa                                     <  1:709233/43‑1 (MQ=255)
                                                            tGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAAt                                    <  1:2514863/42‑1 (MQ=255)
                                                            tGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAAt                                    <  1:1334456/42‑1 (MQ=255)
                                                            tGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAAt                                    <  1:2889314/42‑1 (MQ=255)
                                                            tGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAg        <  1:1198206/70‑1 (MQ=255)
                                                              gCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                    >  1:1839875/1‑56 (MQ=255)
                                                              gCGAAACAGGGAAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTg             >  1:28212/1‑63 (MQ=255)
                                                                 aaaCAGGGTAATCGCTCGGCCAAACGCGTCGAAAAAtt                                   <  1:286097/38‑1 (MQ=255)
                                                                 aaaCAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGgt  >  1:16161/1‑71 (MQ=255)
                                                                   aCAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgc                                 >  1:3874115/1‑38 (MQ=255)
                                                                   aCAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgc                                 >  1:4303159/1‑38 (MQ=255)
                                                                   aCAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgc                                 >  1:3095501/1‑38 (MQ=255)
                                                                   aCAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgc                                 >  1:242029/1‑38 (MQ=255)
                                                                   |                                                                    
CCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGT  >  W3110S.gb/2747491‑2747626

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: