Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 121,347 T→C 70.0% M69V (ATG→GTG)  aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb121,3470TC70.0% 41.7 / 20.1 30M69V (ATG→GTG) aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base C (14/7);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATC  >  W3110S.gb/121280‑121411
                                                                   |                                                                
gCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTc                                                               <  1:2665663/71‑1 (MQ=255)
     aCTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCAcc                                                            >  1:2604874/1‑69 (MQ=255)
      cTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCACCCAg                                                         >  1:2232532/1‑71 (MQ=255)
                aTTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCa                                                              >  1:121540/1‑56 (MQ=255)
                aTTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCa                                                              >  1:145723/1‑56 (MQ=255)
                 tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGa                                              >  1:1613623/1‑71 (MQ=255)
                 tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGa                                              >  1:369186/1‑71 (MQ=255)
                 tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCACACAGCTGACGCATGa                                              >  1:672758/1‑71 (MQ=255)
                   tataAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCAc                                                             >  1:1874059/1‑54 (MQ=255)
                   tataAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGAt                                             >  1:2568010/1‑70 (MQ=255)
                                 cTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATaaa                              <  1:245965/71‑1 (MQ=255)
                                 cTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATaaa                              <  1:785934/71‑1 (MQ=255)
                                 cTAAAGGAGCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATaaa                              <  1:136289/71‑1 (MQ=255)
                                     aGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAc                            >  1:671906/1‑69 (MQ=255)
                                       gAGCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                        <  1:809932/71‑1 (MQ=255)
                                       gAGCCTGCGACAGGTTCTTCGACCACCACTTCACACAGCTGACGCATGATCAGAAAGGCGATAAAACCAg                         >  1:1765003/1‑70 (MQ=255)
                                        aGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                        >  1:2493861/1‑70 (MQ=255)
                                         gCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCaa                      >  1:1483458/1‑71 (MQ=255)
                                         gCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCaa                      >  1:2627602/1‑71 (MQ=255)
                                         gCCTGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCa                       >  1:116923/1‑70 (MQ=255)
                                            tGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCaa                      >  1:1399553/1‑68 (MQ=255)
                                            tGCGACAGGTTCTTCGACCACCACTTCACGCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATgg                   >  1:529240/1‑71 (MQ=255)
                                            tGCGACAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATgg                   >  1:1466209/1‑71 (MQ=255)
                                                aCAGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                        <  1:742622/62‑1 (MQ=255)
                                                  aGGTTCTTCGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAAc             >  1:1500900/1‑71 (MQ=255)
                                                       cttcGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAg         <  1:2304068/70‑1 (MQ=255)
                                                          cGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATTGCGTAACCCAGGAt      <  1:2618157/70‑1 (MQ=255)
                                                          cGACCACCACTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt      <  1:1583457/70‑1 (MQ=255)
                                                            accaccATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCa                       <  1:603435/51‑1 (MQ=255)
                                                             ccaccaCTTCACCCAGCTGACGCACGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATc  <  1:1021716/71‑1 (MQ=255)
                                                              caccacTTCACCCAGCTGACGCATGATCAGAAAGGCGATaa                               <  1:2275620/41‑1 (MQ=255)
                                                              caccacTTCACCCAGCTGACGCATGATCAGAAAGGCGATaa                               <  1:1544795/41‑1 (MQ=255)
                                                              caccacTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCaa                      >  1:384898/1‑50 (MQ=255)
                                                                   |                                                                
GCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATC  >  W3110S.gb/121280‑121411

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: