Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,294,438 T→C 10.6% L395S (TTG→TCG)  yraM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,294,4380TC10.6% 175.6 / 12.7 66L395S (TTG→TCG) yraMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base C (4/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATCGCG  >  W3110S.gb/3294369‑3294503
                                                                     |                                                                 
tAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTg                                                                  >  1:1572546/1‑71 (MQ=255)
tAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTg                                                                  >  1:333228/1‑71 (MQ=255)
     cAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa                                                             <  1:943094/71‑1 (MQ=255)
      aGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAag                                                            <  1:77216/71‑1 (MQ=255)
       ggTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa                                                             >  1:2329198/1‑69 (MQ=255)
        gTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGtt                                                                   <  1:1810743/62‑1 (MQ=255)
         ttCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc                                                         <  1:1702360/71‑1 (MQ=255)
         ttCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGc                                                         <  1:1417300/71‑1 (MQ=255)
          tCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGc                                                                 >  1:4340069/1‑62 (MQ=255)
          tCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGc                                                                 >  1:4334840/1‑62 (MQ=255)
           caggagGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAaga                                                           <  1:2298645/67‑1 (MQ=255)
            agcagGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa                                                             <  1:3118853/64‑1 (MQ=255)
                 gATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTcc                                                 <  1:1616955/71‑1 (MQ=255)
                     gcgcGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc                                                         <  1:2163795/59‑1 (MQ=255)
                         gAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGATCg                                         <  1:2032248/71‑1 (MQ=255)
                          aGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGt                                        <  1:2429191/71‑1 (MQ=255)
                           gTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTa                                       <  1:166508/71‑1 (MQ=255)
                            tATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGtt                                                                   <  1:2470767/42‑1 (MQ=255)
                              ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCt                                                                >  1:2868854/1‑43 (MQ=255)
                              ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACt                                     <  1:1976112/70‑1 (MQ=255)
                              ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg                                    <  1:3341617/71‑1 (MQ=255)
                               tgtgGTCGGTCTGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAagag                                                          >  1:119606/1‑48 (MQ=255)
                                gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTgg                                   >  1:3293590/1‑70 (MQ=255)
                                gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTgg                                   <  1:2640447/70‑1 (MQ=255)
                                 tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc                                  <  1:2451477/70‑1 (MQ=255)
                                 tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc                                  <  1:4522395/70‑1 (MQ=255)
                                 tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc                                  <  1:2598185/70‑1 (MQ=255)
                                    tcggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc                                  <  1:3463361/67‑1 (MQ=255)
                                     cggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACt                               >  1:1726086/1‑69 (MQ=255)
                                      ggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGa                             <  1:4176299/70‑1 (MQ=255)
                                      ggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGa                             >  1:27152/1‑70 (MQ=255)
                                       gtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCt                                              <  1:146362/52‑1 (MQ=255)
                                        tcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAc                           >  1:1006567/1‑70 (MQ=255)
                                        tcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAc                           >  1:4304310/1‑70 (MQ=255)
                                          cGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg                                    >  1:924032/1‑59 (MQ=255)
                                          cGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAg                        >  1:2906663/1‑71 (MQ=255)
                                           gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGa                             <  1:1730537/65‑1 (MQ=255)
                                           gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGc                       <  1:2169685/71‑1 (MQ=255)
                                            ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc                      >  1:1955400/1‑71 (MQ=255)
                                            ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc                      <  1:288238/71‑1 (MQ=255)
                                            ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc                      >  1:2023554/1‑71 (MQ=255)
                                            ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc                      <  1:1332797/71‑1 (MQ=255)
                                            ttGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAc                           >  1:971584/1‑66 (MQ=255)
                                            ttGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc                      >  1:976957/1‑71 (MQ=255)
                                              gCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCg                     >  1:3555336/1‑70 (MQ=255)
                                                tGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTgg                                   >  1:3244128/1‑54 (MQ=255)
                                                tGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCgg                    <  1:2460516/69‑1 (MQ=255)
                                                 gAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag                  >  1:1687017/1‑70 (MQ=255)
                                                 gAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag                  >  1:1912328/1‑70 (MQ=255)
                                                  aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACt                                     >  1:1540741/1‑50 (MQ=255)
                                                  aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag                  <  1:1713540/69‑1 (MQ=255)
                                                  aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGaa                >  1:4314521/1‑71 (MQ=255)
                                                   aaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCa                         <  1:1394849/61‑1 (MQ=255)
                                                   aaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAat               <  1:2269904/71‑1 (MQ=255)
                                                    aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGt                                        <  1:3308640/45‑1 (MQ=255)
                                                    aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGaa                            >  1:836327/1‑57 (MQ=255)
                                                    aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAata              >  1:3387809/1‑71 (MQ=255)
                                                      aTAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATc            >  1:4608042/1‑71 (MQ=255)
                                                       tAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACg                                         <  1:2069577/41‑1 (MQ=255)
                                                        aaCGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgaga                 >  1:522272/1‑64 (MQ=255)
                                                         aCGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGaa                                           >  1:4137259/1‑37 (MQ=255)
                                                         aCGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGaa         <  1:4540947/71‑1 (MQ=255)
                                                           gTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGaaaa       <  1:1872369/71‑1 (MQ=255)
                                                           gTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGaaaa       <  1:4421539/71‑1 (MQ=255)
                                                              gAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag                  <  1:4318029/57‑1 (MQ=255)
                                                              gAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATc     <  1:4328613/70‑1 (MQ=255)
                                                                agagTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgaga                 >  1:3257322/1‑56 (MQ=255)
                                                                agagTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATcgcg  >  1:2137646/1‑71 (MQ=255)
                                                                     |                                                                 
TAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATCGCG  >  W3110S.gb/3294369‑3294503

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: