Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,912,311 G→A 7.6% A89V (GCG→GTG)  htpX ← putative endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,3110GA 7.6% 164.2 / 3.5 53A89V (GCG→GTG) htpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (21/27);  new base A (2/2);  total (23/30)
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.23e-01

TTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCA  >  NC_000913/1912168‑1912458
                                                                                                                                               |                                                                                                                                                   
ttACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTg                                                                                                                                                 <  2:133095/148‑1 (MQ=255)
  aCCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGa                                                                                                                                                <  1:317309/147‑1 (MQ=255)
       cTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGcc                                                                                                                                          <  8:307099/148‑1 (MQ=255)
                     ggCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGt                                                                                                                            <  2:2000/148‑1 (MQ=255)
                       cTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTAt                                                                                                                          >  3:17656/1‑148 (MQ=255)
                       cTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTAt                                                                                                                          >  3:222969/1‑148 (MQ=255)
                         cATGTTCTGCAGCAAACCTGTGCTGACAGCAACAAGAGAGGCATAACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGAGGAATAGCTATCCCCGCCTGACGAGCCTGGGTTGCTACAGTAtt                                                                                                                         >  6:92147/1‑147 (MQ=255)
                            gTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCt                                                                                                                                 >  4:63564/1‑136 (MQ=255)
                            gTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCt                                                                                                                                 <  3:63564/136‑1 (MQ=255)
                               cTGCAGCACACCGGTTCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGGCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCTCCTGACGAGCCTGGGTTGCTACAGTATTGACCAg                                                                                                                   <  5:92147/147‑1 (MQ=255)
                                   agcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAAc                                                                                                              >  8:224687/1‑148 (MQ=255)
                                   agcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAAc                                                                                                              >  5:83742/1‑148 (MQ=255)
                                    gcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCACCTGACGAGCCTGGGTTGCTACAGt                                                                                                                            <  5:43124/133‑1 (MQ=255)
                                    gcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCACCTGACGAGCCTGGGTTGCTACAGt                                                                                                                            <  5:38684/133‑1 (MQ=255)
                                    gcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCACCTGACGAGCCTGGGTTGCTACAGt                                                                                                                            >  6:38684/1‑133 (MQ=255)
                                    gcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCACCTGACGAGCCTGGGTTGCTACAGt                                                                                                                            >  6:43124/1‑133 (MQ=255)
                                       aaCCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTa                                                                                                                           <  7:254668/131‑1 (MQ=255)
                                       aaCCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTc                                                                                                          <  1:144688/148‑1 (MQ=255)
                                                   cAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTAcg                                                                                              <  4:119213/148‑1 (MQ=255)
                                                        aCCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGTTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGACTGGGTTTATACAGTATTGACCAGCCAACGTTCCATTTCGTTACGCGGtt                                                                                         >  8:167978/1‑148 (MQ=255)
                                                             agagGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTTGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCg                                                                                    <  7:167978/148‑1 (MQ=255)
                                                                ggCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATc                                                                                 <  3:118195/148‑1 (MQ=255)
                                                                     cACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTc                                                                            >  4:190385/1‑148 (MQ=255)
                                                                     cACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTc                                                                            >  7:223259/1‑148 (MQ=255)
                                                                            cgcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGcc                                                                                                                                          >  3:191110/1‑79 (MQ=255)
                                                                            cgcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGcc                                                                                                                                          <  4:191110/79‑1 (MQ=255)
                                                                            cgcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCa                                                                     <  5:59274/148‑1 (MQ=255)
                                                                            cgcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCa                                                                     <  4:17656/148‑1 (MQ=255)
                                                                              cACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCAATTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAAc                                                                   <  4:133137/148‑1 (MQ=255)
                                                                                ccGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTc                                                                                                          >  8:188604/1‑107 (MQ=255)
                                                                                ccGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTc                                                                                                          <  7:188604/107‑1 (MQ=255)
                                                                                   gTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGa                                                                                                                       >  5:111724/1‑91 (MQ=255)
                                                                                   gTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGa                                                                                                                       <  6:111724/91‑1 (MQ=255)
                                                                                       cAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTaa                                                          >  7:42185/1‑148 (MQ=255)
                                                                                       cAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTa                                                           <  7:224687/147‑1 (MQ=255)
                                                                                         aaaGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt                                                            <  2:117244/144‑1 (MQ=255)
                                                                                         aaaGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt                                                            >  1:117244/1‑144 (MQ=255)
                                                                                             cGTTGATGTCCGGCGCATGGTAGATAGTCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCAAACAGATCGTAATGccat                                                    >  6:62183/1‑148 (MQ=255)
                                                                                                        ggCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACa                                         <  4:222969/148‑1 (MQ=255)
                                                                                                               ggTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGt                                                                                                            >  7:260977/1‑74 (MQ=255)
                                                                                                               ggTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGt                                                                                                            <  8:260977/74‑1 (MQ=255)
                                                                                                               ggTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACtt                                                                             <  8:55926/105‑1 (MQ=255)
                                                                                                               ggTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACtt                                                                             >  7:55926/1‑105 (MQ=255)
                                                                                                                    aTAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGaaa                             >  4:169848/1‑148 (MQ=255)
                                                                                                                      aGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAAc                            >  6:251155/1‑147 (MQ=255)
                                                                                                                       gCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACg                           <  5:251155/147‑1 (MQ=255)
                                                                                                                       gCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGa                          <  8:223259/148‑1 (MQ=255)
                                                                                                                        ccACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACg                           >  1:20206/1‑146 (MQ=255)
                                                                                                                         cACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAg                        <  3:190385/148‑1 (MQ=255)
                                                                                                                          aCTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAgg                       <  5:62183/148‑1 (MQ=255)
                                                                                                                                         aTCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCACCAGAGCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGa        <  1:239685/148‑1 (MQ=255)
                                                                                                                                            cccGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGCTCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACa     <  6:83742/148‑1 (MQ=255)
                                                                                                                                             ccGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAg    >  8:28647/1‑148 (MQ=255)
                                                                                                                                               gCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCa  <  3:169848/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCA  >  NC_000913/1912168‑1912458

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