Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,818,550 (C)7→8 intergenic (‑50/+55) chbF ← / ← chbR monoacetylchitobiose‑6‑phosphate hydrolase/DNA‑binding transcriptional dual regulator ChbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,818,5431.C93.9% 105.6 / 0.4 33intergenic (‑43/+62)chbF/chbRmonoacetylchitobiose‑6‑phosphate hydrolase/DNA‑binding transcriptional dual regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (17/14);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTAAACG  >  NC_000913/1818416‑1818672
                                                                                                                                |                                                                                                                                   
gATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccccggac                                                                                                                       <  1:111228/140‑6 (MQ=255)
           ataaATCCTTCCAGTCACTCCGGGGCATATCCGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                                                             <  1:277066/140‑1 (MQ=255)
           ataaATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                                                             <  2:372946/140‑1 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAaataat                                                                                                               >  2:10297/1‑129 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAaataat                                                                                                               <  1:10297/129‑1 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                        >  1:447463/1‑136 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                        <  2:447463/136‑1 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAgg                                                                                                    <  1:127970/140‑1 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAgg                                                                                                    >  1:461031/1‑140 (MQ=255)
                    tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCAGATCTGCCACCCGCCAGACATAGATAATCCAGCAACAg                                                                                                     >  1:173522/1‑140 (MQ=255)
                     ccAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                                   <  2:361557/140‑1 (MQ=255)
                     ccAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                                   <  1:612064/140‑1 (MQ=255)
                       aGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACa                                                                                                 <  2:349954/140‑1 (MQ=255)
                           aCTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                                                             <  2:461031/140‑1 (MQ=255)
                              ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtgt                                                                                          <  2:53755/140‑1 (MQ=255)
                                      taGCTGCTCCCGCCACCGATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                                                                  <  2:283180/140‑1 (MQ=255)
                                       aGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                                                                 >  1:390612/1‑140 (MQ=255)
                                         ctgctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGACATAAAAAAATCAACAACAAGGAAGATATGTGAATTTTCCggg                                                                              >  2:667110/1‑139 (MQ=255)
                                            ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                                                                            >  2:590268/1‑140 (MQ=255)
                                              cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATaa                                                                          >  2:27696/1‑140 (MQ=255)
                                                              acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaa                                                          >  2:392586/1‑140 (MQ=255)
                                                                   ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAg                                                     >  2:401007/1‑140 (MQ=255)
                                                                                 tCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGa                                       <  1:88327/140‑1 (MQ=255)
                                                                                             ccTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGc                           >  2:373940/1‑140 (MQ=255)
                                                                                                 cAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTaa                       >  1:568079/1‑140 (MQ=255)
                                                                                                  aGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTaaa                      <  2:97775/140‑1 (MQ=255)
                                                                                                  aGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTaaa                      >  2:335971/1‑140 (MQ=255)
                                                                                                   gTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAg                     <  1:392586/140‑1 (MQ=255)
                                                                                                    tACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAgg                    >  1:431307/1‑140 (MQ=255)
                                                                                                         aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGtt               >  1:384275/1‑140 (MQ=255)
                                                                                                              tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAAttt          >  1:498206/1‑140 (MQ=255)
                                                                                                              tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATATCTTACGTTAGGTCCTAGGCGTAAAGGACGTTAAttt          >  1:239749/1‑140 (MQ=255)
                                                                                                                   tCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTaa     >  1:657008/1‑140 (MQ=255)
                                                                                                                      ggCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTAAACg  >  2:595405/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                   
GATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTAAACG  >  NC_000913/1818416‑1818672

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