Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,728,118 (C)8→9 coding (1254/1455 nt) xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,1101.C85.2% 79.5 / 7.3 27coding (1262/1455 nt)xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (10/13);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATG  >  NC_000913/3727974‑3728231
                                                                                                                                              |                                                                                                                          
cTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑cc                                                                                                                         <  2:1243731/139‑1 (MQ=255)
   aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccc                                                                                                                      <  1:740733/139‑1 (MQ=255)
              aaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCGTACGGt                                                                                                            <  2:456961/139‑1 (MQ=255)
                     cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑cccccccccg                                                                                                                 <  2:658647/126‑3 (MQ=255)
                     cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑cccccccccg                                                                                                                 >  1:658647/1‑124 (MQ=255)
                           tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAgctgct                                                                                                <  2:423895/139‑1 (MQ=255)
                                    cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccccccccg                                                                                                                >  2:1266536/1‑109 (MQ=255)
                                    cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccccccccg                                                                                                                <  1:1266536/112‑4 (MQ=255)
                                               aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccc                                                                                                                      <  2:27578/95‑1 (MQ=255)
                                               aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccc                                                                                                                      >  1:27578/1‑95 (MQ=255)
                                                   aaGTTCCGGAAACAATTAATTGGGCGTTTTTCCTGGTTTCCCCCCGTCCGCCGCCGGCCTTGCTGGCCCCGTGGTGGGCCCACCACC‑‑C‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                        <  2:557075/138‑1 (MQ=255)
                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAgctgct                                                                                                <  1:384994/114‑1 (MQ=255)
                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAgctgct                                                                                                >  2:384994/1‑114 (MQ=255)
                                                      ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                    >  1:215220/1‑139 (MQ=255)
                                                           gCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                              <  2:805880/109‑1 (MQ=255)
                                                           gCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                              >  1:805880/1‑109 (MQ=255)
                                                             aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAc                                                             <  2:215220/139‑1 (MQ=255)
                                                             aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                  <  1:1403529/117‑1 (MQ=255)
                                                             aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                  >  2:1403529/1‑117 (MQ=255)
                                                              acaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCCTACCGAAATCGAGCTGCTGACCGCTGCTATCCGCCACCCTCTGCCCCCAGTAc                                                             >  1:855265/1‑139 (MQ=255)
                                                               caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACt                                                            <  1:214807/139‑1 (MQ=255)
                                                               caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccccccccgta                                                                                                              >  2:1206443/1‑82 (MQ=255)
                                                               caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccccccccgta                                                                                                              <  1:1206443/87‑6 (MQ=255)
                                                                       tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccccccccgtac                                                                                                             >  1:1196363/1‑74 (MQ=255)
                                                                       tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑ccccccccccgtac                                                                                                             <  2:1196363/80‑7 (MQ=255)
                                                                          gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                   <  2:1462369/121‑1 (MQ=255)
                                                                          gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                   >  1:1462369/1‑121 (MQ=255)
                                                                                          cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACACCC‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACg                                 <  2:1328040/139‑1 (MQ=255)
                                                                                           gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat                                                                                    <  2:1014889/87‑1 (MQ=255)
                                                                                           gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat                                                                                    >  1:1014889/1‑87 (MQ=255)
                                                                                           gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                    <  1:575028/103‑1 (MQ=255)
                                                                                           gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                    >  2:575028/1‑103 (MQ=255)
                                                                                           gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTa                               >  1:1498751/1‑139 (MQ=255)
                                                                                                    tGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                       <  2:133883/139‑1 (MQ=255)
                                                                                                                       ccAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa    <  1:690621/139‑1 (MQ=255)
                                                                                                                        cAGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt   <  2:855265/139‑1 (MQ=255)
                                                                                                                         aGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAg                                                                                                     <  1:506294/40‑1 (MQ=255)
                                                                                                                         aGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCGTACGGTAATCGAg                                                                                                     >  2:506294/1‑40 (MQ=255)
                                                                                                                         aGTGCTGG‑‑CCCCACA‑‑TCCCCCCCCCCCCACGGCAATCGAAATGTTGAACGCTTATAACCCCCCGCCTCCGACGCCAGTACTCAAAACCCCGCCCGCCCCCAAACCACGTAACACATTGCCGCTTAATAACGcctgcacg  >  2:1295201/1‑132 (MQ=255)
                                                                                                                                              |                                                                                                                          
CTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGG‑‑CCCCACA‑‑T‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATG  >  NC_000913/3727974‑3728231

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