Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,165,573 (C)6→5 intergenic (+52/‑95) yegP → / → yegQ DUF1508 domain‑containing protein YegP/putative peptidase YegQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,5680C.97.1% 134.2 / ‑2.4 34intergenic (+47/‑100)yegP/yegQDUF1508 domain‑containing protein YegP/putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (16/17);  total (16/18)
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

GATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGG  >  NC_000913/2165434‑2165699
                                                                                                                                      |                                                                                                                                   
gATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGccccca                                                                                                                                <  1:768289/140‑2 (MQ=255)
gATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGccccca                                                                                                                                <  1:577856/140‑2 (MQ=255)
 aTGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGcccccac                                                                                                                               >  1:514979/1‑138 (MQ=255)
      cGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAgcgc                                                                                                                         >  2:840330/1‑140 (MQ=255)
           ccGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACaa                                                                                                                    >  1:816704/1‑140 (MQ=255)
           ccGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACaa                                                                                                                    >  1:405495/1‑140 (MQ=255)
               gcAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAATCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATgcc                                                                                                                <  1:534028/140‑1 (MQ=255)
                cAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATgtc                                                                                                                <  1:414148/140‑3 (MQ=255)
                            aCCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTg                                                                                                   >  1:117524/1‑140 (MQ=255)
                             ccGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTgg                                                                                                  >  2:562876/1‑140 (MQ=255)
                             ccGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTgg                                                                                                  >  2:593741/1‑140 (MQ=255)
                                 aaTTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                                              <  1:636753/140‑1 (MQ=255)
                                  aTTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCa                                                                                             <  2:739402/140‑1 (MQ=255)
                                    tccctcctTTAAAGCCATGGCCAAAAGCCAGCGGGTGAAAGACAATACGAAATTACGTGGCCGGGCGTGGGTTTTTGGCGTCCGGCGCTTTACCAGACCCCCCACAGCGCTACAAGGCCCGCCCTTAAAGTGGGGGCACt                                                                                           <  2:514979/138‑1 (MQ=255)
                                          cGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTa                                                                                     <  2:537022/140‑1 (MQ=255)
                                            ttAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAAc                                                                                   <  2:82789/140‑1 (MQ=255)
                                             tAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAAcc                                                                                  <  2:117524/140‑1 (MQ=255)
                                                      tGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatca                                                                         >  2:526790/1‑140 (MQ=255)
                                                        gCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                                       <  1:476551/140‑1 (MQ=255)
                                                            aaGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGaa                                                                   >  1:898891/1‑140 (MQ=255)
                                                             aGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                                  <  2:279590/140‑1 (MQ=255)
                                                                    cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCt                                                           <  1:562876/140‑1 (MQ=255)
                                                                        gAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAcc                                                       <  1:221773/140‑1 (MQ=255)
                                                                            gACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtc                                                   <  1:840330/140‑1 (MQ=255)
                                                                                 tACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAg                                              <  2:816704/140‑1 (MQ=255)
                                                                                 tACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAg                                              >  1:240430/1‑140 (MQ=255)
                                                                                       aTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACg                                          >  2:791285/1‑138 (MQ=255)
                                                                                       aTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACg                                          <  1:791285/138‑1 (MQ=255)
                                                                                        ttACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                                       <  1:804152/140‑1 (MQ=255)
                                                                                           cGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGaa                                    >  2:915308/1‑140 (MQ=255)
                                                                                            gATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAt                                   <  1:526790/140‑1 (MQ=255)
                                                                                                   ggCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGttt                            >  1:186576/1‑140 (MQ=255)
                                                                                                          gttgttGGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCgg                     <  2:898891/140‑1 (MQ=255)
                                                                                                            tgttgGCGTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                   <  2:405495/140‑1 (MQ=255)
                                                                                                                   gTCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt            >  2:152226/1‑140 (MQ=255)
                                                                                                                    tCCGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTc           >  1:184226/1‑140 (MQ=255)
                                                                                                                      cGGCGCTTTACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTccc         >  2:360691/1‑140 (MQ=255)
                                                                                                                             ttACCAGAGCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCggg  >  1:720365/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGG  >  NC_000913/2165434‑2165699

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