Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,915,028 C→T 6.7% intergenic (+24/‑71) 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,0280CT 6.7% 287.7 / 8.1 90intergenic (+24/‑71)ysaB/glyQhypothetical protein/glycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base T (3/3);  total (47/43)
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.75e-01

TTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAAT  >  W3110S.gb/3914959‑3915098
                                                                     |                                                                      
ttCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACt                                                                       <  1:2222393/71‑1 (MQ=255)
ttCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACt                                                                       <  1:312897/71‑1 (MQ=255)
  cTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                                                      <  1:2222160/70‑1 (MQ=255)
  cTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                                                      >  1:1254243/1‑70 (MQ=255)
     ttGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTATTGtat                                                                   <  1:4179687/70‑1 (MQ=255)
     ttGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtac                                                                   <  1:1344240/70‑2 (MQ=255)
      tGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGtata                                                                  <  1:471750/70‑1 (MQ=255)
       gACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTatat                                                                 <  1:3470372/70‑1 (MQ=255)
       gACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTatat                                                                 <  1:153361/70‑1 (MQ=255)
           cgGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                                                      <  1:673524/61‑1 (MQ=255)
           cgGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtcc                                                            <  1:2596600/71‑1 (MQ=255)
            ggATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtcca                                                           >  1:1537624/1‑71 (MQ=255)
               tGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtccatc                                                         >  1:2898750/1‑70 (MQ=255)
                  ccATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtc                                                             <  1:2931001/63‑1 (MQ=255)
                  ccATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     <  1:581006/71‑1 (MQ=255)
                  ccATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     <  1:742026/71‑1 (MQ=255)
                  ccATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     <  1:952634/71‑1 (MQ=255)
                  ccATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGc                                                     <  1:3542176/71‑1 (MQ=255)
                   cATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCg                                                    <  1:1810019/71‑1 (MQ=255)
                    aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtc                                                             <  1:2007493/61‑1 (MQ=255)
                    aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                   >  1:59472/1‑71 (MQ=255)
                    aTTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCgg                                                   >  1:337061/1‑71 (MQ=255)
                     ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAg                                                      <  1:3148622/67‑1 (MQ=255)
                     ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCggg                                                  >  1:3616327/1‑71 (MQ=255)
                     ttttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATc                                                                <  1:2870367/57‑1 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTATTGTATATCTTCCATCCAGCggg                                                  >  1:4220556/1‑69 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTATTGTATATCTTCCATCCAGCGGGt                                                 <  1:1167584/70‑1 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGta                                                                    >  1:661060/1‑51 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGt                                                                     >  1:2744915/1‑50 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:184412/1‑71 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:1082023/1‑71 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:2901767/1‑71 (MQ=255)
                       ttttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGta                                                >  1:1617887/1‑71 (MQ=255)
                        tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:1581747/1‑70 (MQ=255)
                        tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                <  1:1866086/70‑1 (MQ=255)
                        tttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                <  1:833750/70‑1 (MQ=255)
                         ttGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:2538318/1‑71 (MQ=255)
                           gCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATAc                                             <  1:1121336/70‑1 (MQ=255)
                            cATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTg                                           <  1:228956/71‑1 (MQ=255)
                            cATCTATCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTg                                           <  1:2672976/71‑1 (MQ=255)
                               cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTtc                                                             <  1:2948998/50‑1 (MQ=255)
                               cTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATc                                        >  1:1033574/1‑71 (MQ=255)
                                tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:2306290/1‑64 (MQ=255)
                                tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATcc                                       >  1:4173227/1‑71 (MQ=255)
                                tttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATcc                                       >  1:55925/1‑71 (MQ=255)
                                 ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATccc                                      >  1:2511640/1‑71 (MQ=255)
                                 ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATccc                                      >  1:2916517/1‑71 (MQ=255)
                                 ttCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATccc                                      >  1:1772628/1‑71 (MQ=255)
                                  tCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:3497502/1‑60 (MQ=255)
                                  tCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGAt                                         <  1:1787207/67‑1 (MQ=255)
                                    cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACt                                                                       >  1:2418859/1‑35 (MQ=255)
                                    cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTg                                           >  1:3117821/1‑63 (MQ=255)
                                    cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttc                                   <  1:3972280/71‑1 (MQ=255)
                                    cATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttc                                   >  1:679971/1‑71 (MQ=255)
                                     aTGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                                                      <  1:2297071/35‑1 (MQ=255)
                                     aTGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTg                                                                      <  1:589530/35‑1 (MQ=255)
                                     aTGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     >  1:196744/1‑52 (MQ=255)
                                     aTGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     >  1:1583211/1‑52 (MQ=255)
                                     aTGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttcc                                  <  1:3442423/71‑1 (MQ=255)
                                       gCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGta                                                                    >  1:3028948/1‑35 (MQ=255)
                                        cGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCccttcctt                                <  1:914224/70‑1 (MQ=255)
                                        cGTTAAGCCCTGCTTTTTCCCGCTTCGTACTGTATATCTtcg                                                            >  1:3515465/1‑41 (MQ=39)
                                         gTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCtt                                                              >  1:708499/1‑39 (MQ=255)
                                         gTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTACATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                              <  1:2232836/71‑1 (MQ=255)
                                             aGCCCTGCTTTTTCCCGTTTCGTATTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATc                          <  1:1787174/71‑1 (MQ=255)
                                             aGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATc                          <  1:2797543/71‑1 (MQ=255)
                                             aGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATc                          <  1:1377686/71‑1 (MQ=255)
                                              gCCCTGCTTTTTCCCGTTTCGTATTGTATATCTtccat                                                          >  1:2991712/1‑38 (MQ=255)
                                              gCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGta                                                >  1:2135617/1‑48 (MQ=255)
                                               cccTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     <  1:1337230/42‑1 (MQ=255)
                                                 cTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCcaca                      >  1:1019345/1‑71 (MQ=255)
                                                  tGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGc                                                     <  1:3330589/39‑1 (MQ=255)
                                                  tGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAAt                           <  1:4048059/65‑1 (MQ=255)
                                                   gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:293858/1‑45 (MQ=255)
                                                   gCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:2659240/1‑45 (MQ=255)
                                                    cTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCggg                                                  >  1:1972097/1‑40 (MQ=255)
                                                    cTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACg                    >  1:4056110/1‑70 (MQ=255)
                                                     tttttCCCGTTTCGTATTGTATATCTTCCATCCAGCGGGTATACTGAt                                         >  1:745547/1‑48 (MQ=255)
                                                     tttttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTaaa                            >  1:3248453/1‑61 (MQ=255)
                                                     tttttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAAt                           <  1:1361834/62‑1 (MQ=255)
                                                       tttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:2331660/1‑41 (MQ=255)
                                                        ttCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATcc                                       <  1:4252796/47‑1 (MQ=255)
                                                         tCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTa                              >  1:218499/1‑55 (MQ=255)
                                                          cccGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCa              <  1:1306328/70‑1 (MQ=255)
                                                             gTTTCGTACTGTATATCTTCCATCCAGCGGGtata                                              >  1:4067981/1‑35 (MQ=255)
                                                             gTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCAc          >  1:2972373/1‑71 (MQ=255)
                                                              tttCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAg             <  1:1092727/67‑1 (MQ=255)
                                                               ttCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCcct                                     >  1:594780/1‑42 (MQ=255)
                                                                tCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCcct                                     >  1:1097892/1‑41 (MQ=255)
                                                                tCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGa        >  1:4116856/1‑70 (MQ=255)
                                                                     cTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAaataat  <  1:3986839/71‑1 (MQ=255)
                                                                     |                                                                      
TTCTTTTGACGCGGATGGCCATTTTTTGCATCTTTCCATGCGTTAAGCCCTGCTTTTTCCCGTTTCGTACTGTATATCTTCCATCCAGCGGGTATACTGATCCCTTCCTTTAAATCCACACGTATCCAGCACGAAATAAT  >  W3110S.gb/3914959‑3915098

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: