Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,915,008 C→T 5.5% intergenic (+4/‑91) ysaB → / → glyQ hypothetical protein/glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,915,0080CT 5.5% 226.9 / 2.3 72intergenic (+4/‑91)ysaB/glyQhypothetical protein/glycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (2/2);  total (40/32)
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

CCGTAAAGAGCAGTTTGTCTGTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAAT  >  W3110S.gb/3914938‑3915073
                                                                      |                                                                 
ccGTAAAGAGCAGTTTGTCTGTTCTTTTGACGCGGATGGCCATTTTTTGAATCTTTCCATGAGTTAAGccc                                                                   >  1:3865007/1‑71 (MQ=255)
    aaaGAGCAGTTTGTCTGTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCt                                                               >  1:1246386/1‑71 (MQ=255)
             tttGTCTGTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCttt                                                             <  1:431435/64‑1 (MQ=255)
               tGTCTGTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGttt                                                    >  1:2383654/1‑71 (MQ=255)
                    gTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGccc                                                                   >  1:3827324/1‑51 (MQ=255)
                    gTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGccc                                                                   >  1:4152902/1‑51 (MQ=255)
                    gTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCg                                                  <  1:3741464/68‑1 (MQ=255)
                    gTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCg                                                  <  1:4047400/68‑1 (MQ=255)
                     ttCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCg                                                       <  1:4697811/62‑1 (MQ=255)
                     ttCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACt                                              >  1:2353826/1‑71 (MQ=255)
                      tCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             >  1:3362693/1‑71 (MQ=255)
                      tCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             >  1:1769748/1‑71 (MQ=255)
                       cTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCTTGCTTTTTCCCGTTTCGTACtc                                             >  1:1002231/1‑69 (MQ=255)
                       cTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGtt                                                     >  1:2523150/1‑62 (MQ=255)
                       cTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             <  1:1386110/70‑1 (MQ=255)
                        ttttGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCg                                                  <  1:3263456/64‑1 (MQ=255)
                        ttttGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGt                                                 <  1:2105937/65‑1 (MQ=255)
                         tttGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCg                                                       >  1:3549667/1‑58 (MQ=255)
                           tGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtaca                                         <  1:4550274/70‑3 (MQ=255)
                             aCGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATc                                       <  1:3023552/70‑1 (MQ=255)
                             aCGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCt                                      <  1:2209655/71‑1 (MQ=255)
                               gcgGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtc                                    <  1:4036000/71‑1 (MQ=255)
                                 ggATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCg                                                  <  1:1630365/55‑1 (MQ=255)
                                  gATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             <  1:239041/59‑1 (MQ=255)
                                    tGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGta                                           >  1:3884084/1‑59 (MQ=255)
                                    tGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGta                                           >  1:639022/1‑59 (MQ=255)
                                      gCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTa                                                >  1:1673332/1‑52 (MQ=255)
                                         aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                          >  1:2642391/1‑71 (MQ=255)
                                         aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                          >  1:4096296/1‑71 (MQ=255)
                                         aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                          >  1:3314985/1‑71 (MQ=255)
                                          ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGt                                                 >  1:3620137/1‑47 (MQ=255)
                                          ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccat                                 <  1:2372657/63‑1 (MQ=255)
                                          ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccat                                 <  1:1754150/63‑1 (MQ=255)
                                          ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccat                                 <  1:4630467/63‑1 (MQ=255)
                                           tttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGt                                                      >  1:1207244/1‑41 (MQ=255)
                                           tttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTtccatcca                              >  1:3578575/1‑65 (MQ=255)
                                            ttttGCATGTTTCCATGCGTTAAGCTCTGCTTTTTCCCGTTTCGTACTGTacat                                        >  1:1431297/1‑54 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCttttt                                                           >  1:1255030/1‑35 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTcc                                                         >  1:4988402/1‑37 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             >  1:2110112/1‑49 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             >  1:2984810/1‑49 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTacat                                        >  1:2875569/1‑54 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       >  1:1787487/1‑71 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATc                                       >  1:2667539/1‑55 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                       >  1:2996543/1‑71 (MQ=255)
                                            ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                       >  1:1864590/1‑71 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       <  1:1521399/70‑1 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       >  1:2311354/1‑70 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       >  1:3941380/1‑70 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       <  1:3438760/70‑1 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       >  1:5038911/1‑70 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                       <  1:2716401/70‑1 (MQ=255)
                                             tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                       <  1:1794922/70‑1 (MQ=255)
                                                gCATCTTTCCATGCGTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATAc                    >  1:718155/1‑70 (MQ=255)
                                                 cATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCg                                                  <  1:1044061/39‑1 (MQ=255)
                                                  aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             <  1:3356143/43‑1 (MQ=255)
                                                  aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             <  1:1953873/43‑1 (MQ=255)
                                                  aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                             <  1:500268/43‑1 (MQ=255)
                                                  aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGtata                     <  1:1029262/67‑1 (MQ=255)
                                                  aTCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGtata                     <  1:4919944/67‑1 (MQ=255)
                                                    cTTTCCATGCGTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCggg                         <  1:653318/61‑1 (MQ=255)
                                                    cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                             <  1:2695137/57‑1 (MQ=255)
                                                    cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                             <  1:2289176/57‑1 (MQ=255)
                                                    cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                             <  1:1378241/57‑1 (MQ=255)
                                                     tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGta                                           >  1:3867598/1‑42 (MQ=255)
                                                     tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACt                   >  1:2865101/1‑66 (MQ=255)
                                                      ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCtt                                     >  1:1395177/1‑47 (MQ=255)
                                                      ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCg                           >  1:2240750/1‑57 (MQ=255)
                                                       tCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTAGTGTATATCTTCCATCCAGCGGGTATACTGATccc             >  1:2125753/1‑70 (MQ=255)
                                                           tGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGt                                            <  1:1967270/35‑1 (MQ=255)
                                                             cGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtac                                          >  1:2614077/1‑34 (MQ=255)
                                                             cGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtac                                          >  1:3790136/1‑34 (MQ=255)
                                                              gTTAAGCCTTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa     <  1:1947958/71‑1 (MQ=255)
                                                                 aaGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                     >  1:3182458/1‑52 (MQ=255)
                                                                 aaGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAAt  >  1:1718781/1‑71 (MQ=255)
                                                                      cTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                       <  1:3475198/45‑1 (MQ=255)
                                                                      |                                                                 
CCGTAAAGAGCAGTTTGTCTGTTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAAT  >  W3110S.gb/3914938‑3915073

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: