Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 904,996 G→A 70.5% 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,9960GA70.5% 112.2 / 63.5 78intergenic (‑107/‑19)ybjP/ybjQpredicted lipoprotein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/8);  new base A (34/21);  total (49/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 = 1.00e+00

AATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAA  >  W3110S.gb/904936‑905062
                                                            |                                                                  
aatGATATCCTTTCAATAATAGCGTAACAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAg                                                            <  1:2157629/69‑1 (MQ=255)
 atGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAg                                                            <  1:1990329/68‑1 (MQ=255)
    atatCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACa                                                              >  1:901387/1‑63 (MQ=255)
        ccTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcg                                                     >  1:214990/1‑68 (MQ=255)
         cTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcgat                                                   <  1:3462207/69‑1 (MQ=255)
             cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgat                                                       >  1:1604947/1‑61 (MQ=255)
             cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCAt                                                >  1:2045065/1‑68 (MQ=255)
             cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCAt                                                >  1:804761/1‑68 (MQ=255)
             cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCAt                                                >  1:1629914/1‑68 (MQ=255)
             cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCAt                                                >  1:635852/1‑68 (MQ=255)
             cAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCAt                                                >  1:3461442/1‑68 (MQ=255)
               ataataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGc                                              <  1:2329131/68‑1 (MQ=255)
                  ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCaa                                            >  1:3125940/1‑67 (MQ=255)
                  ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtt                                          >  1:1969006/1‑69 (MQ=255)
                  ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtt                                          <  1:1675328/69‑1 (MQ=255)
                  ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtt                                          >  1:2361934/1‑69 (MQ=255)
                  ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtt                                          >  1:245433/1‑69 (MQ=255)
                  ataGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtt                                          >  1:308538/1‑69 (MQ=255)
                    aGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAttt                                         >  1:2479430/1‑68 (MQ=255)
                    aGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAAttt                                         >  1:1147284/1‑68 (MQ=255)
                    aGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGAATTTTc                                       >  1:2751164/1‑69 (MQ=255)
                     gCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGTAGATCGATCATGCAATTTTc                                       <  1:171851/69‑1 (MQ=255)
                     gCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       <  1:2310531/69‑1 (MQ=255)
                        tATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTAc                                    <  1:1336588/69‑1 (MQ=255)
                        tATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTAc                                    <  1:2653292/69‑1 (MQ=255)
                        tATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTAc                                    <  1:2890468/69‑1 (MQ=255)
                          tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACaa                                  >  1:2293683/1‑69 (MQ=255)
                          tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcg                                                     >  1:2136364/1‑50 (MQ=255)
                          tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAgatcg                                                     >  1:2948867/1‑50 (MQ=255)
                          tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACaa                                  >  1:383932/1‑69 (MQ=255)
                          tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACaa                                  >  1:1721237/1‑69 (MQ=255)
                          tCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACaa                                  >  1:3646977/1‑69 (MQ=255)
                                tGATAATGCTTTTGAGATCGAAGGCTTAGCAAACa                                                              <  1:2995085/35‑1 (MQ=255)
                                tGATAATGCTTTTGAGATCGAAGGCTTAGCAAACa                                                              >  1:1531162/1‑35 (MQ=255)
                                tGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcccc                              >  1:1748325/1‑67 (MQ=255)
                                tGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAAcccc                              >  1:462862/1‑67 (MQ=255)
                                tGATAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCa                             >  1:2706883/1‑68 (MQ=255)
                                  aTAATGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAct                          >  1:2755968/1‑69 (MQ=255)
                                   tAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcg                                                     >  1:2816230/1‑41 (MQ=255)
                                   tAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcg                                                     >  1:304881/1‑41 (MQ=255)
                                   tAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcg                                                     >  1:3103762/1‑41 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGgaga                                                        >  1:2296271/1‑37 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc                         >  1:3694000/1‑68 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc                         >  1:2446042/1‑68 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc                         >  1:1526896/1‑68 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc                         >  1:2306846/1‑68 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc                         >  1:3707816/1‑68 (MQ=255)
                                    aaTGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAct                          >  1:571296/1‑67 (MQ=255)
                                      tGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActct                        >  1:3404137/1‑67 (MQ=255)
                                      tGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       >  1:927236/1‑52 (MQ=255)
                                      tGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       >  1:1643640/1‑52 (MQ=255)
                                      tGCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTc                                       >  1:1703544/1‑52 (MQ=255)
                                       gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTgg                      <  1:3278631/68‑1 (MQ=255)
                                       gCTTTTGAGATCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTgg                      <  1:2824138/68‑1 (MQ=255)
                                        cTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGaa                    >  1:973862/1‑69 (MQ=255)
                                        cTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGaa                    >  1:2395543/1‑69 (MQ=255)
                                        cTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGaa                    >  1:2176214/1‑69 (MQ=255)
                                        cTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGaa                    >  1:232982/1‑69 (MQ=255)
                                          tttGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAgg                  <  1:1953500/69‑1 (MQ=255)
                                          tttGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAgg                  <  1:1957776/69‑1 (MQ=255)
                                             gagaTCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCa               >  1:3285275/1‑69 (MQ=255)
                                              agaTCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCa               <  1:1388615/68‑1 (MQ=255)
                                                aTCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGa             <  1:2661291/68‑1 (MQ=255)
                                                 tCGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGAcc           <  1:2230029/69‑1 (MQ=255)
                                                  cGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtttt                                        >  1:2354270/1‑39 (MQ=255)
                                                  cGAAGGCTTAACAAACAAGGAGATCGATCATGCAAtttt                                        >  1:1816259/1‑39 (MQ=255)
                                                  cGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGc                 <  1:3123272/62‑1 (MQ=255)
                                                  cGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGc                 <  1:1192926/62‑1 (MQ=255)
                                                  cGAAGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGAcc           >  1:1304074/1‑68 (MQ=255)
                                                     aGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc        <  1:1769968/68‑1 (MQ=255)
                                                     aGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc        <  1:1656258/68‑1 (MQ=255)
                                                     aGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc        <  1:777421/68‑1 (MQ=255)
                                                     aGGCTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc        <  1:2967040/68‑1 (MQ=255)
                                                        cTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTg    <  1:229717/69‑1 (MQ=255)
                                                        cTTAACAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTg    <  1:665456/69‑1 (MQ=255)
                                                          tAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  <  1:2334125/69‑1 (MQ=255)
                                                          tAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  <  1:1627072/69‑1 (MQ=255)
                                                          tAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa  <  1:1483195/69‑1 (MQ=255)
                                                            |                                                                  
AATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAA  >  W3110S.gb/904936‑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: