Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,747,549 T→C 7.1% H10H (CAT→CAC yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,747,5490TC 7.1% 147.1 / 3.1 56H10H (CAT→CACyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base C (2/2);  total (27/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 = 7.60e-01

GTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAT  >  W3110S.gb/2747481‑2747617
                                                                    |                                                                    
gtCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACatat                                                                    >  1:2513306/1‑71 (MQ=255)
gtCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACata                                                                     <  1:3093835/70‑1 (MQ=255)
 tCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATg                                                                   <  1:3262703/71‑1 (MQ=255)
 tCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATg                                                                   <  1:4040936/71‑1 (MQ=255)
     ttCAGCATATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCg                                                                >  1:1970528/1‑70 (MQ=255)
      tCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACata                                                                     >  1:255211/1‑64 (MQ=255)
      tCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGa                                                               <  1:510676/70‑1 (MQ=255)
          ccTATTTTTCCGGTTCAGAAACTGGAATGATGACACTCAACCGCTATCGCCTGAGACACATGGCGAAGCa                                                           <  1:4218376/70‑1 (MQ=255)
          ccTATTTTTCCGGTTCAGAAACTGGAATGATGACACTCAACCGCTATCGCCTGAGACACATGGCGAAGCa                                                           <  1:1625872/70‑1 (MQ=255)
                tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTaa                                                     >  1:2686262/1‑70 (MQ=255)
                tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTaa                                                     >  1:413360/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:482916/1‑71 (MQ=255)
                                               cAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                     >  1:2876169/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:2019036/71‑1 (MQ=255)
                                                     cTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGc               <  1:966426/71‑1 (MQ=255)
                                                     cTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGc               >  1:4049286/1‑71 (MQ=255)
                                                        tCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCa                                          <  1:1422268/41‑1 (MQ=255)
                                                         cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaa                              <  1:1993000/52‑1 (MQ=255)
                                                         cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa           >  1:3918556/1‑71 (MQ=255)
                                                         cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa           >  1:1647719/1‑71 (MQ=255)
                                                         cGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa           >  1:2570019/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:2486472/59‑1 (MQ=255)
                                                                gACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                <  1:3073498/59‑1 (MQ=255)
                                                                gACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                <  1:1652704/59‑1 (MQ=255)
                                                                  cATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                     <  1:2351533/52‑1 (MQ=255)
                                                                  cATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAt  <  1:4325109/71‑1 (MQ=255)
                                                                   atatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAg            <  1:2112346/60‑1 (MQ=255)
                                                                   atatGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAg            <  1:4102493/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)
                                                                    |                                                                    
GTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAT  >  W3110S.gb/2747481‑2747617

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: