Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,165,590 Δ1 bp intergenic (+69/‑78) yegP → / → yegQ UPF0339 family protein/putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,5880C.97.8% 202.0 / 0.1 47intergenic (+67/‑80)yegP/yegQUPF0339 family protein/putative peptidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (21/24);  total (21/26)
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.54e-01

GCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAA  >  NC_000913/2165490‑2165682
                                                                                                  |                                                                                              
gCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATgcc                                                                                               >  1:289678/1‑100 (MQ=255)
  acaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATgccgc                                                                                             <  2:116185/100‑3 (MQ=255)
  acaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATgccgc                                                                                             <  1:402941/100‑3 (MQ=255)
  acaAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATgccg                                                                                              <  1:122203/99‑2 (MQ=255)
     aGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCt                                                                                         >  1:983185/1‑100 (MQ=255)
         aGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTaaa                                                                                     <  1:840721/100‑1 (MQ=255)
          gACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAg                                                                                    <  1:269608/100‑1 (MQ=255)
            cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCGGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTg                                                                                  >  1:598868/1‑100 (MQ=255)
            cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTg                                                                                  >  2:146238/1‑100 (MQ=255)
              gTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTggg                                                                                >  2:561833/1‑100 (MQ=255)
               tGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTgggg                                                                               >  2:1268856/1‑100 (MQ=255)
                gAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTggggg                                                                              >  1:173350/1‑100 (MQ=255)
                      cAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTcc                                                                        <  1:9010/100‑1 (MQ=255)
                        aTACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTcccc                                                                      <  1:159782/100‑1 (MQ=255)
                         tACGTAATGAAGAGGCCGGGCGTGCGTTGTTGGCGTCCGGCGCTTGACCAGAGCAGCCCACAGCGCTACAGTGCCGCCCTTAAAGTGGGGGCACTcc                                                                        <  1:635295/97‑1 (MQ=255)
                          aCGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTa                                                                    <  2:974341/100‑1 (MQ=255)
                                 tACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCtt                                                             <  1:831015/100‑1 (MQ=255)
                                  aCGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtc                                                            >  1:1300950/1‑100 (MQ=255)
                                     aTGCCGGACGTGGGTTGTTGGCGTCCGGCGCTTTACCGGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatc                                                         <  2:1238338/100‑1 (MQ=255)
                                      tGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatca                                                        <  1:300946/100‑1 (MQ=255)
                                        ccGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                      >  2:763630/1‑100 (MQ=255)
                                          gggCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGt                                                     <  1:79986/100‑1 (MQ=255)
                                          gggCGTGGGTTGTTGGCGTCCGGAGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTg                                                    >  2:1189372/1‑100 (MQ=255)
                                           ggCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCTCCCACAGCGCTACAAGGCCGCCCGTAAGGTGGGGGCACTCCGCTAGCCGCGTCATCAGGTg                                                    <  2:386375/99‑1 (MQ=255)
                                             cGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                 <  2:1300939/100‑1 (MQ=255)
                                             cGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAg                                                 <  1:357325/100‑1 (MQ=255)
                                               tGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCg                                               <  2:1078904/100‑1 (MQ=255)
                                                gggTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCgg                                              <  2:1242188/100‑1 (MQ=255)
                                                    tgttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCt                                          <  1:63453/100‑1 (MQ=255)
                                                       tgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                       >  2:256655/1‑100 (MQ=255)
                                                       tgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                       >  2:1043874/1‑100 (MQ=255)
                                                        ggCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGCAGCGGATCTGagc                                      <  1:884412/100‑3 (MQ=255)
                                                        ggCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAcc                                      <  1:763635/100‑1 (MQ=255)
                                                               ggCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatc                               <  1:1093016/100‑1 (MQ=255)
                                                               ggCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatc                               >  1:691155/1‑100 (MQ=255)
                                                                   cTTTACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGaa                           >  2:183528/1‑100 (MQ=255)
                                                                    tttACCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAAc                          >  2:1058774/1‑100 (MQ=255)
                                                                       aCCAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgag                       >  1:635841/1‑100 (MQ=255)
                                                                        ccAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                      >  1:631647/1‑100 (MQ=255)
                                                                        ccAGAGCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgag                       >  1:701384/1‑99 (MQ=255)
                                                                           gagCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGaa                   >  1:1305751/1‑100 (MQ=255)
                                                                           gagCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCCGACCTGTCATCAGAACGAGAGaa                   <  1:451892/100‑1 (MQ=255)
                                                                            agCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAt                  <  1:1058784/100‑1 (MQ=255)
                                                                             gCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCGGAACGAGAGAAtt                 >  1:423057/1‑100 (MQ=255)
                                                                             gCCCCCCACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAtt                 <  2:826550/100‑1 (MQ=255)
                                                                              ccccccACAGCGCTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTa                <  2:691151/100‑1 (MQ=255)
                                                                                       gcgcTACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc       >  1:892535/1‑100 (MQ=255)
                                                                                           tACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGa   >  1:794011/1‑100 (MQ=255)
                                                                                           tACAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGa   <  2:289675/100‑1 (MQ=255)
                                                                                            aCAATGCCGCCCTTAGAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa  <  1:1281986/100‑1 (MQ=255)
                                                                                            aCAATGCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa  <  1:1281644/100‑1 (MQ=255)
                                                                                                  |                                                                                              
GCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAA  >  NC_000913/2165490‑2165682

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