Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,818,550 (C)7→8 intergenic (‑50/+55) chbF ← / ← chbR phospho‑chitobiase, general 6‑phospho‑beta‑glucosidase activity/repressor of chb operon for N,N'‑diacetylchitobiose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,818,5431.C92.9% 44.7 / ‑2.3 14intergenic (‑43/+62)chbF/chbRphospho‑chitobiase, general 6‑phospho‑beta‑glucosidase activity/repressor of chb operon for N,N'‑diacetylchitobiose utilization
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (8/5);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

TCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTAAACGT  >  NC_000913/1818412‑1818673
                                                                                                                                    |                                                                                                                                   
tCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGTTATATCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTACAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                                                                             <  1:473386/139‑1 (MQ=255)
           cTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACataaat                                                                                                                   <  2:510865/139‑1 (MQ=255)
                        tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg                                                                                                      >  1:420169/1‑139 (MQ=255)
                           aGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                                                                   <  1:75038/139‑1 (MQ=255)
                                      ggTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGaa                                                                                        <  1:84174/139‑1 (MQ=255)
                                                                     acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtgt                                                                                           >  1:144644/1‑105 (MQ=255)
                                                                     acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtgt                                                                                           <  2:144644/105‑1 (MQ=255)
                                                                                                             aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGt                 >  2:203619/1‑139 (MQ=255)
                                                                                                                  tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAAtt            >  1:165092/1‑139 (MQ=255)
                                                                                                                  tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAAtt            >  2:482211/1‑139 (MQ=255)
                                                                                                                       tCTGGCATATCTGCCCCGCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTa       >  2:321143/1‑139 (MQ=255)
                                                                                                                         tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGctct                                 <  1:112706/111‑1 (MQ=255)
                                                                                                                         tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGctct                                 >  2:112706/1‑111 (MQ=255)
                                                                                                                         tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTaaa     >  1:64437/1‑139 (MQ=255)
                                                                                                                            cATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTAAACGt  >  1:312596/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                   
TCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTAAACGT  >  NC_000913/1818412‑1818673

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