Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,013,011 C→T 7.0% D235N (GAC→AAC)  yqeC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,013,0110CT 7.0% 170.4 / 2.2 57D235N (GAC→AAC) yqeCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/23);  new base T (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGTT  >  W3110S.gb/3012943‑3013077
                                                                    |                                                                  
cTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcgt                                                                  <  1:1001695/71‑1 (MQ=255)
cTAATTCACAAATCTGCGCGCGATTGCAGGATCTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcgt                                                                  <  1:3412563/71‑1 (MQ=255)
   aTTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgtt                                                               >  1:157094/1‑71 (MQ=255)
     tCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgttgt                                                             >  1:1962116/1‑71 (MQ=255)
      cacaAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgttgt                                                             >  1:173014/1‑70 (MQ=255)
      cacaAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTgttgtt                                                            >  1:622218/1‑71 (MQ=255)
      cacaAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTgttgtt                                                            >  1:626662/1‑71 (MQ=255)
       acaAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTgttgtt                                                            >  1:933866/1‑70 (MQ=255)
          aaTCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAg                                                        >  1:3308629/1‑71 (MQ=255)
                cgcgcgATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg                                                  <  1:2515794/71‑1 (MQ=255)
                cgcgcgATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg                                                  <  1:3303490/71‑1 (MQ=255)
                      aTTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGa                                                 <  1:3160276/66‑1 (MQ=255)
                          cAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                                            >  1:3481894/1‑67 (MQ=255)
                          cAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                        >  1:1083695/1‑71 (MQ=255)
                            ggATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                      <  1:2683370/71‑1 (MQ=255)
                             gATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTTGTGTTGTTGCAGCGGTTGAAg                                               <  1:2035279/61‑1 (MQ=255)
                             gATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTTGTGTTGTTGCAGCGGTTGAAg                                               <  1:2890127/61‑1 (MQ=255)
                             gATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                           <  1:2015724/65‑1 (MQ=255)
                             gATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                           <  1:398349/65‑1 (MQ=255)
                              aTGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTTGTGTTGTTGCAGCGGTTGAAGGAGCTCGCttt                                    >  1:1773953/1‑71 (MQ=255)
                               tGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                         <  1:967751/49‑1 (MQ=255)
                                gTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgtt                                                               <  1:1968025/42‑1 (MQ=255)
                                gTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgtt                                                               <  1:1542127/42‑1 (MQ=255)
                                   cTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAgt                                              >  1:1181360/1‑55 (MQ=255)
                                   cTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGc                                <  1:2430491/70‑1 (MQ=255)
                                    ttGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCaa                              >  1:2905827/1‑71 (MQ=255)
                                      gTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                            >  1:3282229/1‑71 (MQ=255)
                                       tataTCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCgg                                                     >  1:2336320/1‑45 (MQ=255)
                                         tatCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                         <  1:2417409/71‑1 (MQ=255)
                                         tatCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                         <  1:3216547/71‑1 (MQ=255)
                                           tCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCaa                              >  1:675890/1‑64 (MQ=255)
                                            cGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtc                      >  1:1193257/1‑71 (MQ=255)
                                            cGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtc                      >  1:3822239/1‑71 (MQ=255)
                                             gCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct                     >  1:3578436/1‑71 (MQ=255)
                                             gCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct                     >  1:986039/1‑71 (MQ=255)
                                             gCCAAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct                     >  1:1452860/1‑71 (MQ=255)
                                                  gccagATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg                                                  <  1:2247536/37‑1 (MQ=255)
                                                  gccagATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg                                                  <  1:917438/37‑1 (MQ=255)
                                                    cagATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACAtt               >  1:2121656/1‑70 (MQ=255)
                                                     agATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                               >  1:3830912/1‑37 (MQ=255)
                                                      gATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                               <  1:2691921/36‑1 (MQ=255)
                                                      gATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                               <  1:210890/36‑1 (MQ=255)
                                                      gATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg                                               <  1:3750959/36‑1 (MQ=255)
                                                       actGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCttt                                    <  1:554698/44‑1 (MQ=255)
                                                         tGCCTCTACGTTGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGCGCAATCGCATTCTCACATTGAGaa          >  1:1034552/1‑70 (MQ=255)
                                                         tGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGaa          >  1:3842105/1‑70 (MQ=255)
                                                         tGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGaa          >  1:2604041/1‑70 (MQ=255)
                                                         tGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGaa          >  1:2326337/1‑70 (MQ=255)
                                                              ctACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACgg    >  1:3683563/1‑71 (MQ=255)
                                                              ctACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACgg    >  1:679132/1‑71 (MQ=255)
                                                              ctACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACgg    >  1:1459799/1‑71 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                      <  1:1969379/36‑1 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                      <  1:1666703/36‑1 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                      <  1:1511530/36‑1 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGc                                <  1:3757625/42‑1 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGc                                <  1:122859/42‑1 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                            >  1:126756/1‑46 (MQ=255)
                                                               tACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct                     >  1:3764852/1‑53 (MQ=255)
                                                                aCGTCGTGTTGTTGCAGCGGTTGAAGGAGCTAGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGt   >  1:1005392/1‑70 (MQ=255)
                                                                 cgtcgtGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGtt  <  1:2264735/70‑1 (MQ=255)
                                                                  gtcgtGTTGTTTCAGCGGTTGAAGGAGCTCGCTTTGTGc                                <  1:3653276/39‑1 (MQ=255)
                                                                    |                                                                  
CTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGTT  >  W3110S.gb/3012943‑3013077

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: