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.C91.2% 104.6 / 3.4 34coding (1262/1455 nt)xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (14/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCA  >  NC_000913/3727974‑3728242
                                                                                                                                          |                                                                                                                                      
cTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑cc                                                                                                                                     <  2:1244345/139‑1 (MQ=255)
  gATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑cccc                                                                                                                                   <  1:94464/139‑1 (MQ=255)
          gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑cccccccccgta                                                                                                                          <  2:855343/139‑5 (MQ=255)
                       cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAgct                                                                                                              <  2:933191/139‑1 (MQ=255)
                        aTCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAgctg                                                                                                             <  1:1307284/139‑1 (MQ=255)
                          cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAgctgct                                                                                                           <  2:55665/139‑1 (MQ=255)
                                    ggACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                 <  1:1089983/138‑1 (MQ=255)
                                    cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGat                                                                                                 >  1:1005836/1‑139 (MQ=255)
                                      aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat                                                                                               >  2:1336439/1‑139 (MQ=255)
                                      aCTGGTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCACCTTACGGTAATCGAGCTGCTGACCGCTGagat                                                                                               >  2:1043678/1‑139 (MQ=255)
                                       cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCGCAGTGCTGGCCCCACAT‑‑cc                                                                                                                                     <  2:981820/100‑1 (MQ=255)
                                       cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATc                                                                                              >  2:131159/1‑139 (MQ=255)
                                       cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑cc                                                                                                                                     >  1:981820/1‑100 (MQ=255)
                                          ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                           <  1:131159/139‑1 (MQ=255)
                                              aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAg                                                                                        <  1:475926/139‑1 (MQ=255)
                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACC‑‑CCCCCCCCGTAAGGTAATCGAGCTGCTGCCCGCTGCTATCCGCCAGCATCTGa                                                                                 >  1:75074/1‑139 (MQ=255)
                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                 >  2:1458486/1‑139 (MQ=255)
                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccc                                                                                                                               <  1:721534/92‑1 (MQ=255)
                                                    aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccc                                                                                                                               >  2:721534/1‑92 (MQ=255)
                                                      ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCTAgctg                                                                                                             >  2:886694/1‑109 (MQ=255)
                                                      ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAgctg                                                                                                             <  1:886694/109‑1 (MQ=255)
                                                      ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                               >  2:446696/1‑139 (MQ=255)
                                                            caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                         <  1:446696/139‑1 (MQ=255)
                                                                          gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                           <  1:1458486/139‑1 (MQ=255)
                                                                          gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                           >  2:437578/1‑139 (MQ=255)
                                                                          gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGAGATCCGCCAGCATCTGCCGCCAGTACTCACTACGcgccc                                                           >  2:1376723/1‑139 (MQ=255)
                                                                                            ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCAATCAACGTaa                                         >  2:383931/1‑139 (MQ=255)
                                                                                                   cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                  <  2:1080763/139‑1 (MQ=255)
                                                                                                           gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                          <  2:75074/139‑1 (MQ=255)
                                                                                                            ccTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACtt                                    <  2:1476140/128‑1 (MQ=255)
                                                                                                            ccTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACtt                                    >  1:1476140/1‑128 (MQ=255)
                                                                                                             cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                        <  1:1336439/139‑1 (MQ=255)
                                                                                                                  tgcGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg                  <  1:311268/139‑1 (MQ=255)
                                                                                                                      cccAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAg                                                                                                                <  2:567662/42‑1 (MQ=255)
                                                                                                                      cccAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAg                                                                                                                >  1:567662/1‑42 (MQ=255)
                                                                                                                         aGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCCGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACCTAACACTTTCCGGCTTAATACCGCAAGCATg             >  2:1343440/1‑139 (MQ=255)
                                                                                                                               ggCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACa      <  1:1376723/139‑1 (MQ=255)
                                                                                                                                gCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                            >  2:569002/1‑68 (MQ=255)
                                                                                                                                gCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                            <  1:569002/68‑1 (MQ=255)
                                                                                                                                 ccccACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATc    >  2:810811/1‑139 (MQ=255)
                                                                                                                                   ccACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa  <  2:274300/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
CTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCA  >  NC_000913/3727974‑3728242

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