Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 904,996 G→A 84.0% intergenic (‑107/‑19) ybjP ← / → ybjQ predicted lipoprotein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb904,9960GA84.0% 183.2 / 27.6 75intergenic (‑107/‑19)ybjP/ybjQpredicted lipoprotein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base A (38/25);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 3.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAA  >  W3110S.gb/904929‑905062
                                                                   |                                                                  
cGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACa                                                                  <  1:3408594/70‑1 (MQ=255)
cGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACa                                                                  <  1:3046422/70‑1 (MQ=255)
   taataaTGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAAc                                                               >  1:21975/1‑70 (MQ=255)
   taataaTGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAAc                                                               >  1:2033585/1‑70 (MQ=255)
   taataaTGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAAc                                                               >  1:1601452/1‑70 (MQ=255)
   taataaTGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAAc                                                               >  1:3213278/1‑70 (MQ=255)
   taataaTGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAAc                                                               >  1:3078754/1‑70 (MQ=255)
     ataatGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAg                                                            >  1:3707613/1‑71 (MQ=255)
        atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAg                                                            >  1:634204/1‑68 (MQ=255)
        atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAg                                                            >  1:167323/1‑68 (MQ=255)
        atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAg                                                            >  1:610006/1‑68 (MQ=255)
        atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGgag                                                         <  1:253370/71‑1 (MQ=255)
        atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGgag                                                         <  1:2276532/71‑1 (MQ=255)
        atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGgag                                                         <  1:2303159/71‑1 (MQ=255)
         tGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGgaga                                                        <  1:1773891/71‑1 (MQ=255)
         tGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGgaga                                                        <  1:1158912/71‑1 (MQ=255)
         tGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGgaga                                                        <  1:2455628/71‑1 (MQ=255)
              tCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAg                                                            >  1:1987669/1‑62 (MQ=255)
              tCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcgat                                                   <  1:3164362/71‑1 (MQ=255)
              tCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcgat                                                   <  1:2956896/71‑1 (MQ=255)
              tCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcgat                                                   <  1:2664372/71‑1 (MQ=255)
              tCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcgat                                                   <  1:2271367/71‑1 (MQ=255)
                   tCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATg                                               <  1:1185279/70‑1 (MQ=255)
                   tCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATg                                               <  1:2776816/70‑1 (MQ=255)
                    cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGc                                              >  1:2665117/1‑70 (MQ=255)
                    cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGc                                              >  1:993321/1‑70 (MQ=255)
                       taataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAt                                           >  1:1255299/1‑70 (MQ=255)
                         ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTaacaaaca                                                              <  1:2305181/49‑1 (MQ=255)
                         ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTaacaaaca                                                              <  1:624279/49‑1 (MQ=255)
                         ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTaacaaaca                                                              <  1:1186868/49‑1 (MQ=255)
                         ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAttt                                         <  1:2586889/70‑1 (MQ=255)
                          taGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTc                                       <  1:1415218/71‑1 (MQ=255)
                          taGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGgaga                                                        <  1:2471213/54‑1 (MQ=255)
                           aGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACaaa                                                                >  1:274281/1‑45 (MQ=255)
                            gCGTATCAGTCTGATAATGCTTTTGAGATCGAAGTCTTAACAAACAAGGAGATCGATCATGCAATTTTCTa                                     >  1:2971434/1‑71 (MQ=255)
                            gCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTa                                     >  1:2543020/1‑71 (MQ=255)
                            gCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTa                                     >  1:1469358/1‑71 (MQ=255)
                             cGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTAc                                    <  1:3535347/71‑1 (MQ=255)
                              gTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACa                                   <  1:50271/71‑1 (MQ=255)
                              gTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACa                                   <  1:1958481/71‑1 (MQ=255)
                              gTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACa                                   <  1:468283/71‑1 (MQ=255)
                                 tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                >  1:641209/1‑71 (MQ=255)
                                 tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                >  1:652430/1‑71 (MQ=255)
                                 tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                >  1:913556/1‑71 (MQ=255)
                                 tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                >  1:1453892/1‑71 (MQ=255)
                                 tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                >  1:1258996/1‑71 (MQ=255)
                                 tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcc                                >  1:1823436/1‑71 (MQ=255)
                                   aGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAAccc                               >  1:3589508/1‑70 (MQ=255)
                                     tCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcg                                                     <  1:858465/46‑1 (MQ=255)
                                      cTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcgat                                                   <  1:255061/47‑1 (MQ=255)
                                           aaTGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTc                                       <  1:938315/54‑1 (MQ=255)
                                            aTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       >  1:3334902/1‑53 (MQ=255)
                                            aTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       >  1:247326/1‑53 (MQ=255)
                                            aTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       >  1:3364857/1‑53 (MQ=255)
                                             tGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAc                                 >  1:1032277/1‑58 (MQ=255)
                                             tGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAc                                 >  1:773608/1‑58 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:2597017/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:411870/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:3608706/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:335397/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:2290647/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:2673445/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:3113556/1‑71 (MQ=255)
                                              gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                   >  1:169486/1‑71 (MQ=255)
                                                   tGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCa               <  1:1228832/70‑1 (MQ=255)
                                                   tGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCa               <  1:1502256/70‑1 (MQ=255)
                                                     agaTCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGAc            <  1:590707/71‑1 (MQ=255)
                                                     agaTCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGAc            <  1:1905667/71‑1 (MQ=255)
                                                       aTCGAAGGCTTAACAAACAAGGAGATCGATCATTCAATTTTCTACAACCCCAActc                         <  1:3581233/56‑1 (MQ=255)
                                                       aTCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc                         <  1:2046096/56‑1 (MQ=255)
                                                          gAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACaa                                  >  1:1104542/1‑44 (MQ=255)
                                                          gAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACaa                                  >  1:2999245/1‑44 (MQ=255)
                                                             ggCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTg    <  1:2523441/71‑1 (MQ=255)
                                                             ggCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTg    <  1:2519301/71‑1 (MQ=255)
                                                               cTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  <  1:284528/71‑1 (MQ=255)
                                                               cTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  <  1:59205/71‑1 (MQ=255)
                                                                ttAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  >  1:325744/1‑70 (MQ=255)
                                                                ttAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  >  1:2416632/1‑70 (MQ=255)
                                                                   |                                                                  
CGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAA  >  W3110S.gb/904929‑905062

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: