Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 4,241,676 A→G 8.9% S165S (TCT→TCC)  xylE ← D‑xylose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,241,6760AG 8.9% 104.7 / 2.3 44S165S (TCT→TCC) xylED‑xylose transporter
Reads supporting (aligned to +/- strand):  ref base A (15/24);  new base G (2/2);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

TTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATAAACAGGCACAGT  >  NC_000913_3_pae_tpiA/4241528‑4241822
                                                                                                                                                    |                                                                                                                                                  
ttAAGAACAGCAGTGAAGGGATGAGCTGCGGGGCAAACATAGAACGCCAGCCGTCAGTATTCAGCCAGCTGGAATAACCGGACCGGGCAATAAAATAGTTTACGAAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTaaaa                                                                                                                                                    <  2:150195/149‑1 (MQ=255)
           agTGCAGGTATACATCCCGAGGCAAACATATAACCCCAGCCGTCAGTAGTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGttt                                                                                                                                         <  1:591732/149‑1 (MQ=255)
                aGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCTCAAACTGGTTAAAAGAGACCAGTTTCCcgc                                                                                                                                    >  1:601053/1‑149 (MQ=255)
                  ggATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGa                                                                                                                                  >  1:108085/1‑149 (MQ=255)
                   gATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGaa                                                                                                                                 >  2:637404/1‑149 (MQ=255)
                   gATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGaa                                                                                                                                 >  2:637224/1‑149 (MQ=255)
                   gATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGaa                                                                                                                                 <  1:163508/149‑1 (MQ=255)
                   gATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGaa                                                                                                                                 <  1:299708/149‑1 (MQ=255)
                       cATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATg                                                                                                                             <  1:115086/149‑1 (MQ=255)
                        atcccGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGa                                                                                                                            >  2:684313/4‑149 (MQ=255)
                         ttCCGAGGTAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACGGGCACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAg                                                                                                                           <  2:12665/149‑1 (MQ=255)
                         ttCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCACTAAAATAGTTTACGCAGTAAACTAAAAATTGCCCAAAAATAATCACAACATGGCGAAAAGAGACCAGTTTCCCGCGAATATGac                                                                                                                           >  2:20257/1‑148 (MQ=255)
                         ttCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGac                                                                                                                           <  1:734658/149‑2 (MQ=255)
                                 cAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAGGAGACCAGTTTCCCGCGaa                                                                                                                                 >  2:573208/1‑135 (MQ=255)
                                 cAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAGGAGACCAGTTTCCCGCGaa                                                                                                                                 >  2:455409/1‑135 (MQ=255)
                                 cAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAGGAGACCAGTTTCCCGCGaa                                                                                                                                 <  1:455409/135‑1 (MQ=255)
                                 cAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAGGAGACCAGTTTCCCGCGaa                                                                                                                                 <  1:573208/135‑1 (MQ=255)
                                  aaaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGGTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCa                                                                                                                  >  1:285137/1‑149 (MQ=255)
                                   aaCATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCCGTTTCCCGCCAATATGCGCTGGCGCCAg                                                                                                                 >  2:297817/1‑149 (MQ=255)
                                   aaCACATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGGTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAg                                                                                                                 <  2:555150/149‑1 (MQ=255)
                                             gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAATTGAGGCCGGTTTCCCGCGAGTATGATCTGGAGTCAGTTTCggtata                                                                                                       >  1:59343/1‑144 (MQ=255)
                                             gCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAata                                                                                                       <  2:447641/149‑1 (MQ=255)
                                               cAGCCGTCAGGATTAAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACCGGTTAATAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAAtata                                                                                                     <  1:105112/149‑1 (MQ=255)
                                                      cAGTATTCAGCCAGCTGGCATCAACGGAACGGGCAATAAAATAGTTGACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGc                                                                                              <  1:6972/149‑1 (MQ=255)
                                                           ttCAGCCAGCTGGCATCACCGGACCGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGCCAAATTTCCCACGAATAGAAGCTCACACTAGCTCCGCAAAATACATTgaccgaaa                                                                                         >  2:299811/1‑142 (MQ=255)
                                                               gccagcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATCCGAATTATCTCCAGGTCCGCAATATACATTGGCGAGAGCAtt                                                                                     <  1:356904/149‑1 (MQ=255)
                                                               gccagcTGGCAACACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCAtt                                                                                     <  1:375706/149‑1 (MQ=255)
                                                                  agcTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAg                                                                                  <  2:309082/149‑1 (MQ=255)
                                                                             ccGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAAc                                                                       <  2:130321/149‑1 (MQ=255)
                                                                             ccGGAACGGGCAATAAAATAGGTTACGACGTCAAGCAGAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAAc                                                                       >  1:729732/1‑149 (MQ=255)
                                                                                      gCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCAGCAATATACATTGGCGAGAACATTGCAGCTAAACCAACCCCAATAcc                                                              >  2:12363/1‑149 (MQ=255)
                                                                                          tAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCa                                                          >  2:222019/1‑149 (MQ=255)
                                                                                           aaaaTAGTTTACGCAGTAAACTAAAAGTTGGCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACTTTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCaa                                                         <  2:221729/149‑1 (MQ=255)
                                                                                                  tttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGGGACCAGTTTCCCTCGAATCGGAGCCGGAGCCAGTTCCGCAATATACATTGGCGAAAGCATTGAAGTCAACACCACGCCCATACCGCAaaaaataca                                                  >  2:273655/1‑146 (MQ=255)
                                                                                                   ttACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGa                                                 <  1:660268/149‑1 (MQ=255)
                                                                                                         aGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGCGAATTTAGTGTAAACCAACGCCAATACCGCCAATAATGCGataaata                                           >  1:18250/1‑149 (MQ=255)
                                                                                                                   aaGTTGCCCGAAAAGAACCGCAAACTGGTGAAAAGAGACCAGTTGCCCGCGAATATGAGCTGGAGCCCGTTCCGCAATATACATTGGCGCGAGCATGGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCg                                 <  2:109274/149‑1 (MQ=255)
                                                                                                                      ttGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGa                              <  1:589580/149‑1 (MQ=255)
                                                                                                                                ataatCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGc                    <  1:306329/149‑1 (MQ=255)
                                                                                                                                ataagCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATCGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGc                    <  1:154206/149‑1 (MQ=255)
                                                                                                                                  aatCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCa                  <  1:122683/149‑1 (MQ=255)
                                                                                                                                       cAAACTGGTTAGAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATaa             >  1:449000/1‑149 (MQ=255)
                                                                                                                                           cTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATAAACAg         <  1:575323/149‑1 (MQ=255)
                                                                                                                                            tGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAAGATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATAAACAgg        <  2:394069/149‑1 (MQ=255)
                                                                                                                                              gttagcagagACCAGGTGCCCTCGAAGATGAGCGGGAGCCAGTTGCGCACTATACATTGGCGAGAGCATTGAGGCTAAACAAACGCCAATACCGCCAATAATGCGATAAATAACAGATTCCGGGACATAACCTGCCAGATAAACAGGca      <  2:234026/143‑1 (MQ=255)
                                                                                                                                                tAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATAAACAGGcaca    <  1:16659/149‑1 (MQ=255)
                                                                                                                                                  aaaGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATAAACAGGCACAgt  <  1:640152/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTAAGAACAGCAGTGCAGGGATACATTCCGAGGCAAACATATAACGCCAGCCGTCAGTATTCAGCCAGCTGGCATCACCGGAACGGGCAATAAAATAGTTTACGCAGTAAACTAAAAGTTGCCCGAAAATAATCGCAAACTGGTTAAAAGAGACCAGTTTCCCGCGAATATGAGCTGGAGCCAGTTCCGCAATATACATTGGCGAGAGCATTGAGGCTAAACCAACGCCAATACCGCCAATAATGCGATAAATAACAAATTCCGGGACATAACCTGCCAGATAAACAGGCACAGT  >  NC_000913_3_pae_tpiA/4241528‑4241822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑