Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,336 A→C 20.0% intergenic (+34/+26) ptsG → / ← fhuE fused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,3360AC20.0% 52.8 / 12.4 25intergenic (+34/+26)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (10/10);  new base C (2/3);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

GCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTA  >  NC_000913/1159199‑1159439
                                                                                                                                         |                                                                                                       
gcgcAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                               <  3:136581/148‑1 (MQ=255)
                               tgtTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGc                                                                                                     >  1:563464/1‑111 (MQ=255)
                               tgtTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGc                                                                                                     <  2:563464/111‑1 (MQ=255)
                                         aTTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                       <  1:896368/99‑1 (MQ=255)
                                         aTTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                       >  2:896368/1‑99 (MQ=255)
                                                       aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                       <  1:38997/85‑1 (MQ=255)
                                                       aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                       >  2:38997/1‑85 (MQ=255)
                                                                                 atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGa                                                                                                         >  1:1091424/1‑57 (MQ=255)
                                                                                 atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGa                                                                                                         <  2:1091424/57‑1 (MQ=255)
                                                                                          tCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                      <  3:793724/51‑1 (MQ=255)
                                                                                          tCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                      <  3:70122/51‑1 (MQ=255)
                                                                                          tCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                      >  4:70122/1‑51 (MQ=255)
                                                                                          tCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                      >  4:793724/1‑51 (MQ=255)
                                                                                                   aCTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                    >  1:1093692/1‑108 (MQ=255)
                                                                                                   aCTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                    <  2:1093692/108‑1 (MQ=255)
                                                                                                                               aGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACg                            >  2:57280/1‑88 (MQ=255)
                                                                                                                               aGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACg                            <  1:57280/88‑1 (MQ=255)
                                                                                                                                 gCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTa  <  2:853225/112‑1 (MQ=255)
                                                                                                                                 gCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTa  >  1:853225/1‑112 (MQ=255)
                                                                                                                                     ccAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgc                                           <  2:830987/67‑1 (MQ=255)
                                                                                                                                     ccAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgc                                           >  1:830987/1‑67 (MQ=255)
                                                                                                                                       aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGt                                 <  3:445242/75‑1 (MQ=255)
                                                                                                                                       aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGt                                 >  4:445242/1‑75 (MQ=255)
                                                                                                                                       aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAg                               <  4:876049/77‑1 (MQ=255)
                                                                                                                                       aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAg                               >  3:876049/1‑77 (MQ=255)
                                                                                                                                       aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                           <  2:826519/81‑1 (MQ=255)
                                                                                                                                       aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                           >  1:826519/1‑81 (MQ=255)
                                                                                                                                         aTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGcc                                                              <  4:503711/44‑1 (MQ=255)
                                                                                                                                         aTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGcc                                                              >  3:503711/1‑44 (MQ=255)
                                                                                                                                         |                                                                                                       
GCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTA  >  NC_000913/1159199‑1159439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: