Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,1102.C29.2% 32.9 / 25.9 24coding (1262/1455 nt)xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base . (9/8);  new base C (4/3);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAC  >  NC_000913/3727980‑3728234
                                                                                                                                    |                                                                                                                             
aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑cccccccc                                                                                                                       <  1:203737/139‑1 (MQ=255)
     cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccccccgta                                                                                                                  <  2:281078/139‑6 (MQ=255)
               cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAg                                                                                                         <  2:1126685/139‑1 (MQ=255)
                    cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAgctgct                                                                                                    <  1:778418/139‑1 (MQ=255)
                    cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑c                                                                                                                              >  1:1122622/1‑112 (MQ=255)
                    cTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑c                                                                                                                              <  2:1122622/112‑1 (MQ=255)
                           aTGCTACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGGAATCTAGCTGCTGACCGc                                                                                              <  1:866852/139‑1 (MQ=255)
                                 cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATc                                                                                       >  1:346900/1‑139 (MQ=255)
                                 atGTTCTAACGGTAGTTGCGGCAAAAATTCAATGAGCGATTTCTCTGGGTTCGCCGCGTTCCGCGTCCGCCTTGCTGCGCCCGGAGGTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat                                                                                        <  2:955278/138‑1 (MQ=255)
                                         aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGc                                                                                >  1:78278/1‑139 (MQ=255)
                                            gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccc                                                                                                                          >  1:198710/1‑92 (MQ=255)
                                            gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccc                                                                                                                          <  2:198710/92‑1 (MQ=255)
                                              aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCCAGCGGCTGACCGCTGATATCCGCCAGCATCTg                                                                           >  1:497255/1‑139 (MQ=255)
                                              aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                          >  2:260508/1‑139 (MQ=255)
                                                  gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                       <  2:78278/139‑1 (MQ=255)
                                                     gCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGt                                                                   >  1:949089/1‑139 (MQ=255)
                                                                   aGCGATTTCTCTGGATTCGCCGCGAACTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGTTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcg‑cc                                                    >  2:70506/1‑139 (MQ=255)
                                                                    gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGata                                                                                         >  2:647420/1‑102 (MQ=255)
                                                                    gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGata                                                                                         <  1:647420/102‑1 (MQ=255)
                                                                            tctGGATTTGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGACCCACATC‑CCCCCCCCGTACGGTAATCGCGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCa                                           <  1:675601/139‑1 (MQ=255)
                                                                              tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAAt                                         <  1:241873/139‑1 (MQ=255)
                                                                              tGGATTCGCCGCGAACTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGTTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAAt                                         <  1:70506/139‑1 (MQ=255)
                                                                                     gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTa                                  >  1:634857/1‑139 (MQ=255)
                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCCCGCCCGGCAATCGCGCCGCTCACCGCTGATTTTCCCCACCACCTCTCTCCCGTCCCCCCTACCCGCCCCCGCCCCCatcccccca                                   >  2:753238/1‑131 (MQ=255)
                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCCTACGGTAACCGCACTGCCGGCCCCTGCGATCCGCCAGCATCCTGCGCCAGTACTCACCCCCCG‑cccccccccccaccaccgca                                  >  2:462050/1‑129 (MQ=255)
                                                                                                 gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                     <  2:130253/93‑1 (MQ=255)
                                                                                                 gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                     >  1:130253/1‑93 (MQ=255)
                                                                                                         tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc               <  1:462050/139‑1 (MQ=255)
                                                                                                          gctggGCCCAGTGCTGGCCCCACAT‑‑cccccccccc                                                                                                                     <  1:112574/35‑3 (MQ=255)
                                                                                                          gctgcGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg             >  2:929325/1‑139 (MQ=255)
                                                                                                          gctgcGCCCAGTGCTGGCCCCACAT‑‑cccccccccc                                                                                                                     >  2:112574/1‑33 (MQ=255)
                                                                                                             gcgcCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg          <  2:949089/139‑1 (MQ=255)
                                                                                                                     tGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc  >  1:157352/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                             
AGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCG‑CCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAC  >  NC_000913/3727980‑3728234

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