Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,800,301 A→C 25.0% F119C (TTC→TGC)  pheS ← phenylalanine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,800,3010AC25.0% 59.0 / 37.6 52F119C (TTC→TGC) pheSphenylalanine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (16/23);  new base C (5/8);  total (21/31)
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

GCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCA  >  W3110S.gb/1800232‑1800370
                                                                     |                                                                     
gCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCgaaga                                                                       >  1:3638270/1‑70 (MQ=255)
gCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCgaaga                                                                       >  1:4040290/1‑70 (MQ=255)
  aTCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAAc                                                                    >  1:3585341/1‑71 (MQ=255)
     gAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACtt                                                                  <  1:3310161/70‑1 (MQ=255)
     gAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACttt                                                                 >  1:3332966/1‑71 (MQ=255)
      aaGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTc                                                                <  1:3620032/71‑1 (MQ=255)
           aTGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATAc                                                           >  1:2569971/1‑71 (MQ=255)
            tGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATAc                                                           <  1:1056600/70‑1 (MQ=255)
             gATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACgg                                                         >  1:1103935/1‑71 (MQ=255)
               tAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAAc                                                                    >  1:1281766/1‑58 (MQ=255)
               tAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTc                                                       <  1:3375176/71‑1 (MQ=255)
               tAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTc                                                       <  1:1635662/71‑1 (MQ=255)
                aaTCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTc                                                       >  1:2142944/1‑70 (MQ=255)
                 aTCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGa                                                     <  1:2120692/71‑1 (MQ=255)
                  tCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGAt                                                    <  1:921644/71‑1 (MQ=255)
                    gTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATgg                                                  >  1:3986706/1‑71 (MQ=255)
                    gTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATgg                                                  >  1:2086271/1‑71 (MQ=255)
                    gTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATgg                                                  <  1:2962864/71‑1 (MQ=255)
                    gTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATgg                                                  <  1:475184/71‑1 (MQ=255)
                    gTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATgg                                                  <  1:1349159/71‑1 (MQ=255)
                     tCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGa                                                              <  1:1747275/58‑1 (MQ=255)
                     tCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGa                                                              <  1:3811799/58‑1 (MQ=255)
                       ttCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATGGTAc                                               >  1:620719/1‑71 (MQ=255)
                       ttCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTAc                                               >  1:3102632/1‑71 (MQ=255)
                       ttCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTAc                                               >  1:3931701/1‑71 (MQ=255)
                       ttCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTAc                                               >  1:2900563/1‑71 (MQ=255)
                       ttCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTAc                                               >  1:2952409/1‑71 (MQ=255)
                         cGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGa                                                     >  1:3763693/1‑63 (MQ=255)
                           aTTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATGGTACGGGt                                           <  1:1114910/71‑1 (MQ=255)
                              tCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAAc                                        >  1:2124448/1‑71 (MQ=255)
                                cGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:399861/71‑1 (MQ=255)
                                cGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:700946/71‑1 (MQ=255)
                                cGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:31177/71‑1 (MQ=255)
                                cGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:2060839/71‑1 (MQ=255)
                                          gCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGa                                                     >  1:925902/1‑46 (MQ=255)
                                          gCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGa                                                     >  1:2228157/1‑46 (MQ=255)
                                                gTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:641870/55‑1 (MQ=255)
                                                gTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:1629971/55‑1 (MQ=255)
                                                gTAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:2880966/55‑1 (MQ=255)
                                                 tAAAGCCAAGCTCACCGAAGCAACTTTCGATACGGTCGAt                                                    <  1:3645802/40‑1 (MQ=255)
                                                  aaaGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGttt                     >  1:862346/1‑70 (MQ=255)
                                                  aaaGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGttt                     >  1:2025473/1‑70 (MQ=255)
                                                   aaGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                   <  1:1175585/71‑1 (MQ=255)
                                                   aaGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                   <  1:3882382/71‑1 (MQ=255)
                                                   aaGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                   <  1:2583908/71‑1 (MQ=255)
                                                     gCCAAGCTCACCGAAGAAACTTTCGATACGGTCGa                                                     >  1:1766837/1‑35 (MQ=255)
                                                     gCCAAGCTCACCGAAGAAACTTTCGATACGGTCGa                                                     >  1:1386877/1‑35 (MQ=255)
                                                     gCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                   <  1:2961502/69‑1 (MQ=255)
                                                       cAAGCTCACCGAAGAAACTTTCGATACGGTCTATGGTACGGGTAACCg                                      <  1:1606716/48‑1 (MQ=255)
                                                       cAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:3190011/48‑1 (MQ=255)
                                                       cAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:1370780/48‑1 (MQ=255)
                                                       cAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:3498610/48‑1 (MQ=255)
                                                         aGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                      <  1:3500097/46‑1 (MQ=255)
                                                                cGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                   <  1:850972/58‑1 (MQ=255)
                                                                    gaaACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCa  <  1:1627885/71‑1 (MQ=255)
                                                                     aaaCTTTCGATTCGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCa  <  1:3539785/70‑1 (MQ=255)
                                                                     aaaCTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCa  <  1:3049451/70‑1 (MQ=255)
                                                                     |                                                                     
GCATCGAAGTTATGATAATCGTCTTCGATTTCCGGCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCA  >  W3110S.gb/1800232‑1800370

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: