Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,839,384 T→A intergenic (‑120/‑140) ascG ← / → ascF asc operon transcriptional repressor, prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,839,3840TA98.0% 154.3 / NA 51intergenic (‑120/‑140)ascG/ascFasc operon transcriptional repressor, prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (29/21);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATAA  >  NC_000913/2839253‑2839517
                                                                                                                                   |                                                                                                                                     
cATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                                <  2:387308/139‑1 (MQ=255)
   cgtcgtCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGCAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCtt                                                                                                                             >  2:91094/1‑139 (MQ=255)
    gtcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTg                                                                                                                            >  2:962660/1‑139 (MQ=255)
         catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGCAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTAt                                                                                                                       <  1:91094/139‑1 (MQ=255)
                                      aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                                    >  1:105249/1‑97 (MQ=255)
                                      aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                                    <  2:105249/97‑1 (MQ=255)
                                       cGCAGTGTACTACACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                                    >  2:75208/1‑96 (MQ=255)
                                       cGCAGTGTACTACACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                                    <  1:75208/96‑1 (MQ=255)
                                           gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                                     <  1:962660/139‑1 (MQ=255)
                                           gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAAACGTCGGGAAAAATATTCAGGTGACCGTTTTCACAAATATAAAAAATGAACAAATCGCTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                                     <  1:1177489/139‑1 (MQ=255)
                                               aCTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                           <  1:471299/97‑1 (MQ=255)
                                               aCTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                           >  2:471299/1‑97 (MQ=255)
                                                 tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                                >  2:1167492/1‑90 (MQ=255)
                                                 tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                                <  1:1167492/90‑1 (MQ=255)
                                                 tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTAt                                                                                                                       >  2:783266/1‑99 (MQ=255)
                                                 tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTAt                                                                                                                       <  1:783266/99‑1 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACGATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  1:1055297/1‑139 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                                  >  2:383575/1‑81 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                                  <  1:383575/81‑1 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGt                                                                                                                  <  2:462089/97‑1 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGt                                                                                                                  >  1:462089/1‑97 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAc                                                                                        <  1:293396/123‑1 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  1:361902/1‑139 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  2:137650/1‑139 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTAATCATGATTTTGTGAAAc                                                                                        >  2:293396/1‑123 (MQ=255)
                                                        tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGATACACAAATATAAAAAATGAACAAATCAATCTCTTGCTTATTTAGTGACAACTATTCATGAATTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  2:518334/1‑139 (MQ=255)
                                                        tttGCCGATTGTCCGTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATTATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  1:322760/1‑139 (MQ=255)
                                                                     cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGt                                                                                                                  >  1:311241/1‑84 (MQ=255)
                                                                     cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGt                                                                                                                  <  2:311241/84‑1 (MQ=255)
                                                                             aTCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc                                                   >  2:434280/1‑139 (MQ=255)
                                                                             aTCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc                                                   >  1:30661/1‑139 (MQ=255)
                                                                             aTCGGCGGGAAAAATATTCAGGTGACCGGATTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc                                                   >  2:1205161/1‑139 (MQ=255)
                                                                              tCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCg                                                  <  2:30661/139‑1 (MQ=255)
                                                                                    ggAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGt                                            <  2:159883/139‑1 (MQ=255)
                                                                                            atTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGAt                                    <  1:434280/139‑1 (MQ=255)
                                                                                                aGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATaaca                                >  1:1042399/1‑139 (MQ=255)
                                                                                                  gTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATaacaac                              >  2:528861/1‑139 (MQ=255)
                                                                                                   tGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATaacaac                              >  2:170973/1‑138 (MQ=255)
                                                                                                   tGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATaacaac                              <  1:170973/138‑1 (MQ=255)
                                                                                                        ggTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATc                        <  1:137650/139‑1 (MQ=255)
                                                                                                          tttCACAAATATAAAAAATGAACAATTCACTCGCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTa                      <  1:518334/139‑1 (MQ=255)
                                                                                                            tCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTAcc                    >  1:1080552/1‑139 (MQ=255)
                                                                                                            tCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTAcc                    >  2:894525/1‑139 (MQ=255)
                                                                                                                 aaTATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGACACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACt               <  1:1205161/139‑1 (MQ=255)
                                                                                                                    ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGttt                                                                     <  2:159604/82‑1 (MQ=255)
                                                                                                                    ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGttt                                                                     >  1:159604/1‑82 (MQ=255)
                                                                                                                    ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGAt                                    <  2:488560/115‑1 (MQ=255)
                                                                                                                    ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGAt                                    >  1:488560/1‑115 (MQ=255)
                                                                                                                             gAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATa   >  2:1037157/1‑139 (MQ=255)
                                                                                                                              aaCAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATaa  >  2:104262/1‑139 (MQ=255)
                                                                                                                              aaCAAATCACTCTCTTGCTTATTTAGTGACAAATATTCATGATTTTGTGAAACCGGATTCTTAATTCGGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAAGAACTAAATCTACACTACTAATAAAAGGATaa  >  1:57009/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
CATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATAA  >  NC_000913/2839253‑2839517

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