Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,165,573 (C)6→7 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,5671.C87.5% 120.9 / 11.0 40intergenic (+46/‑101)yegP/yegQDUF1508 domain‑containing protein YegP/putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (15/20);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

AGCCAGATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAA  >  NC_000913/2165429‑2165682
                                                                                                                                           |                                                                                                                    
agccagATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGc                                                                                                                     >  1:125723/1‑140 (MQ=255)
agccagATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGc                                                                                                                     >  1:212033/1‑140 (MQ=255)
     gATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGcccccc                                                                                                               <  2:200509/140‑1 (MQ=255)
       tGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGccccccca                                                                                                             >  2:640282/1‑138 (MQ=255)
       tGTACGCCACCGCGCAAACTCGTGAAACCGGTATTCCCTCCGTTAAAGTCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGccccccca                                                                                                             <  2:283880/140‑3 (MQ=255)
         tACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGcccccccaca                                                                                                           <  1:528401/140‑5 (MQ=255)
           cGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTGGTTGGCGCCCGTCGTTTTACCGGAGCCCCCCCACAgc                                                                                                          >  1:283880/1‑140 (MQ=255)
            gCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAgcg                                                                                                         >  2:502212/1‑140 (MQ=255)
                ccGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTAc                                                                                                     >  2:99783/1‑140 (MQ=255)
                   cgcAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAAt                                                                                                  >  1:999160/1‑140 (MQ=255)
                      aaTCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATgcc                                                                                               <  1:39971/140‑1 (MQ=255)
                      aaTCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATgcc                                                                                               <  1:122234/140‑1 (MQ=255)
                      aaTCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATgcc                                                                                               <  2:812909/140‑1 (MQ=255)
                        tctcGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATgcccg                                                                                             >  1:9458/1‑140 (MQ=255)
                                      aaTTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCGCCCCCACGGCGGTAGAATGCCAGCCCTTAAAGTgggg                                                                               >  1:275013/1‑140 (MQ=255)
                                        ttGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                             <  1:525643/140‑1 (MQ=255)
                                         tGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCa                                                                            <  1:176019/140‑1 (MQ=255)
                                          gCCTCCGTTAAAGCCAAGGGCATAAGCCAGAAGGTAAAAGACAATACGAAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCAc                                                                           <  2:395689/140‑1 (MQ=255)
                                             tCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                             <  1:175401/135‑1 (MQ=255)
                                             tCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGc                                                                             >  2:175401/1‑135 (MQ=255)
                                             tCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTcc                                                                        >  1:653198/1‑140 (MQ=255)
                                                 ttAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTa                                                                    <  1:99783/140‑1 (MQ=255)
                                                 ttAAAGCCAATGGCACAAGCCACACGGTGAACGACCATAAGTAATTACAATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTa                                                                    <  2:275013/140‑1 (MQ=255)
                                                    aaGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTTCAATGCCGGCCCTTGAAGTGGGGGCACTCCCCTAAcc                                                                 >  1:359403/1‑140 (MQ=255)
                                                         aaTGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtc                                                            >  1:772465/1‑140 (MQ=255)
                                                         aaTGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCACGCGCTTAAAGTGGGGGAACTCCCCTAACCGCTtc                                                            >  2:349085/1‑140 (MQ=255)
                                                           tGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatc                                                         >  1:676506/1‑140 (MQ=255)
                                                           tGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatc                                                         >  1:881736/1‑140 (MQ=255)
                                                               ccaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                      <  2:608575/139‑1 (MQ=255)
                                                                   gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGaa                                                  >  1:956355/1‑140 (MQ=255)
                                                                    ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                 <  2:9458/140‑1 (MQ=255)
                                                                    ccAGACGGTGAAAGACAATACGTAATTACGAGGCCTGGCATGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                 <  2:373546/140‑1 (MQ=255)
                                                                     cAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCg                                               <  2:881736/140‑1 (MQ=255)
                                                                     cAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCg                                               <  2:676506/140‑1 (MQ=255)
                                                                              aaaGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                       <  2:772465/140‑1 (MQ=255)
                                                                                          tAATTACGATGCCGGACGTGGGTTGATGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGaa                           <  1:349085/140‑1 (MQ=255)
                                                                                           aaTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGGAc                          >  2:610812/1‑140 (MQ=255)
                                                                                               aCGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                      <  1:502212/140‑1 (MQ=255)
                                                                                                   tGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCATGTGAAGCTGATGTGACCTGTCCTCAGAACGAGAGAAt                  >  1:359547/1‑140 (MQ=255)
                                                                                                   tGCCGGGCGTGGGTTGTTGGCGACCGGCGCTGTACCAGAGCCCCCCCAAAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAt                  <  2:359403/140‑1 (MQ=255)
                                                                                                           gTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTGAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa          <  2:212033/140‑1 (MQ=255)
                                                                                                           gTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa          <  2:125723/140‑1 (MQ=255)
                                                                                                              ggTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc       <  2:653198/140‑1 (MQ=255)
                                                                                                                   ttgGCGTCCGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa  >  1:797395/1‑140 (MQ=255)
                                                                                                                          ccGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa          <  2:745171/125‑1 (MQ=255)
                                                                                                                          ccGGCGCTTTACCAGAGCCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa          >  1:745171/1‑125 (MQ=255)
                                                                                                                                           |                                                                                                                    
AGCCAGATGTACGCCACCGCGCAATCTCGTGAAACCGGAATTGCCTCCGTTAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAA  >  NC_000913/2165429‑2165682

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: