Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,313,027 T→G intergenic (‑36/‑323) alsR ← / → rpiB DNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,0270TG94.3% 119.7 / ‑6.2 35intergenic (‑36/‑323)alsR/rpiBDNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (15/18);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTG  >  NC_000913/4312789‑4313244
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                         
gATACCTTAACTATCATCGCGTCAGATACCGCCAGAGCTTCTGCGACAGCTTTAATTGCGGGTGCACAACTCATGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTCGCATTCCTTCCTGCTTCATTCGCAGGTAAGGCTCTAACCCTATACCGTGCTGAAGCCCTGAATCAAACCCTGACTGGCACGTCACGTTCTCTTGTCCTTGAATTTAAGTCTCTACTTGATCACCACAGCa                                                                                                                                                                                                               <  2:86081/251‑1 (MQ=255)
              catcGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGa                                                                                                                                                                                                  <  1:7968/250‑1 (MQ=255)
                              gCCAGAGCTTCTGCGCCATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGCTGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTcaca                                                                                                                                                                                  >  1:33053/1‑250 (MQ=255)
                               ccAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTCCCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACaa                                                                                                                                                                                 >  2:96408/1‑250 (MQ=255)
                                  gatcTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGAGCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTca                                                                                                                                                                                    <  2:313900/241‑1 (MQ=255)
                                  gagCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTca                                                                                                                                                                                    >  1:313900/1‑244 (MQ=255)
                                    gCTTCTGCGACATCTTTAATTGCGGGTTCCCAACTCAGGTTACCGGGTTTGAGTAACCACGCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGTCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACGAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAaaaat                                                                                                                                                                           <  2:443814/251‑1 (MQ=255)
                                               aTCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTtata                                                                                                                                                                 >  1:318822/1‑250 (MQ=255)
                                                     aaTTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCTGCTGTCATTTTCTGTCATTCCTTCCTGCTTCAATCGCAGGTAAGGCGCTAACCCTAGACCGTGCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAAGTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCGCAAAAAAATAAAACTGAGAAAGttt                                                                                                                                                           <  2:317746/250‑1 (MQ=255)
                                                      aTTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttat                                                                                                                                                         <  1:239271/251‑1 (MQ=255)
                                                        tGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttatt                                                                                                                                                        <  1:155595/250‑1 (MQ=255)
                                                         gCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttatttt                                                                                                                                                      <  1:265457/251‑1 (MQ=255)
                                                          cGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttatttt                                                                                                                                                      >  1:387530/1‑250 (MQ=255)
                                                                        aGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGc                                                                                                                                       >  1:266391/1‑251 (MQ=255)
                                                                        aGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTAATTTTATGTCCTTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACACTCTACCGTTCGGAAGCGATGACTCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGc                                                                                                                                       >  2:196645/1‑251 (MQ=255)
                                                                                             cACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTGATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCACCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACACCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATACAATTTATAAAGTTTATTTTATTGTGTACATTGGCGAGTGGCTTAACATTTTAttt                                                                                                                  >  1:271352/1‑251 (MQ=255)
                                                                                                 ccACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCGTTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACtt                                                                                                                                                                                      >  1:84152/1‑179 (MQ=255)
                                                                                                 ccACGATGCGGCTTTCATTTTCTGTCATGCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACGAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACtt                                                                                                                                                                                      <  2:84152/179‑1 (MQ=255)
                                                                                                    cGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTaa                                                                                                           <  1:54277/251‑1 (MQ=255)
                                                                                                     gATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATAACCACAGCAAGAACAGAAAAGATGCTGGCGTCTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAAt                                                                                                          <  1:354703/251‑1 (MQ=255)
                                                                                                                 aTTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAAc                                                                                              <  1:100950/251‑1 (MQ=255)
                                                                                                                          cATTCCTTCCTGCTTCATACGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCGGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATATATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTa                                                                                     <  2:266391/251‑1 (MQ=255)
                                                                                                                           atgccttcctGCGTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAAGCAAACTCTGACTGGCTCATCACGGTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTGCACAAAAAAATAAAATTTATAAAGTTTATTGTATTGTGAACAATGGCGAGTGGCTTAAAATTTGATTTCATTTAATTAGTTTAAAAACATGATGTTAt                                                                                    <  2:194564/248‑1 (MQ=255)
                                                                                                                             tccttcctGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTAttt                                                                                  >  1:398060/1‑251 (MQ=255)
                                                                                                                                            tCGCAGGTAAGGCGCTAACCCTATACCGTTCCGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttt                                                                    <  2:111269/250‑1 (MQ=255)
                                                                                                                                             cGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGCATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGCACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt                                                                  >  2:502022/1‑251 (MQ=255)
                                                                                                                                                aggaaaGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATCTCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt                                                                  >  2:138957/5‑248 (MQ=255)
                                                                                                                                                 ggTAAGGCGCTAACCCTTTACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTATCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt                                                                  >  1:436762/1‑247 (MQ=255)
                                                                                                                                                    aaGGCGCTAACCCTATTCCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAgatga                                                            <  2:110788/250‑1 (MQ=255)
                                                                                                                                                    aaGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACAGGATGTTATTTGTTTTTCTGGTTTTTTAgatgat                                                           <  2:318822/251‑1 (MQ=255)
                                                                                                                                                                       gTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTgg                                        >  2:146522/1‑250 (MQ=255)
                                                                                                                                                                          cGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAgt                                      <  1:359777/250‑1 (MQ=255)
                                                                                                                                                                                            ctctGACTGGCTCATAACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAACTTTATCAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTgcgc                   >  1:342373/1‑251 (MQ=255)
                                                                                                                                                                                                            aCGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAActct   <  1:106883/251‑1 (MQ=255)
                                                                                                                                                                                                             cGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGCTCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTg  >  2:217351/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                         
GATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTG  >  NC_000913/4312789‑4313244

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑