Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,106,201 C→T Q128* (CAG→TAG)  cpxP → periplasmic protein CpxP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,106,2010CT83.3% 85.1 / 7.2 36Q128* (CAG→TAG) cpxPperiplasmic protein CpxP
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (14/16);  total (17/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 = 1.00e+00

AGAGCCATATGTTCGACGGCATAAGTTTAACCGAACATCAGCGTCAGCAGATGCGAGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTA  >  NC_000913/4105944‑4106469
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                            
agagCCATATGTTCGACGGCATAAGTTTAACCGAACATCAGCGTCAGCAGATGCGAGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAAc                                                                                                                                                                                                                                                      >  1:323600/1‑282 (MQ=255)
  agCCATATGTTCGACGGCATAAGTTTAACCGAACATCAGCGTCAGCAGATGCGCGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGTAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAGAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACAt                                                                                                                                                                                                                                                    >  1:393900/1‑282 (MQ=255)
        aagaggccagccATAACATTAACCGAACATCAGCATCAGCGGATGCGGGAGCGTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGAGAATGAGCAAATTGCTCGTCAGGTGGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACgaga                                                                                                                                                                                                                                                         <  2:52359/260‑1 (MQ=255)
                    taaacttaaaCCGAACATCAGGCGCAGCAGAGGGGGGATATAAAGAAACAGGCCCGGCACGAAAAGCCGCCTGTTAAAGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAGAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACAAAATGTATCGCATGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGaa                                                                                                                                                                                                                                          <  2:54649/266‑1 (MQ=255)
                          taaaCCGAACATCAGCGTCAGCAGATGCGAGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATg                                                                                                                                                                                                                                        <  1:112850/268‑1 (MQ=255)
                             aCCGAACATCAGCGTCAGCAGATGCGAGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAAc                                                                                                                                                                                                                                                      >  1:52359/1‑253 (MQ=255)
                                                  gaggcgatcctgaTGAAACAGGCCAGGAACTAACAGCCTACTGATAAGGTTGTCGAACAGGAGAGAATGCAACGACTTGCAGCCGGAGAAAATTATGATAAAAACGCTGGGGGCTCACAGGCAGAAAAAAGGGCGAATGAGCCAATTGCTGGTCAGGTTAAGAAGGGAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCTACACCGAAGGGAGCAGTTGcg                                                                                                                                                                                                                             <  2:70191/245‑1 (MQ=255)
                                                    gCGAGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTc                                                                                                                                                                                                  >  1:133019/1‑282 (MQ=255)
                                                            ttATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTACCGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGcgtga                                                                                                                                                                                                                          >  1:70191/1‑250 (MQ=255)
                                                                    cAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTg                                                                                                                                                                                  >  1:382487/1‑282 (MQ=255)
                                                                    cAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAAAGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCGACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTg                                                                                                                                                                                  >  1:30588/1‑282 (MQ=255)
                                                                                                                           gCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGacac                                                                                                                           <  1:282529/282‑1 (MQ=255)
                                                                                                                           gCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGacac                                                                                                                           <  1:150413/282‑1 (MQ=255)
                                                                                                                               tGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCa                                                                                                                         >  1:274073/1‑280 (MQ=255)
                                                                                                                                           aaaTTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGCAGCAAGc                                                                                                                                                                                                                                                                       >  2:173581/1‑126 (MQ=255)
                                                                                                                                           aaaTTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGCAGCAAGc                                                                                                                                                                                                                                                                       <  1:173581/126‑1 (MQ=255)
                                                                                                                                                             tgggaggggacagacagAAAAAAGGTCGAAGTAGCAAATTGCTCGTCAAGTTGGGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGTAGTAGCAAGCGGTTTTAAACGAGAAACATCGACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATcc                                                                                                                      <  2:30588/244‑1 (MQ=255)
                                                                                                                                                                                    tGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACAt                                                                                                                                                                                                                                                    >  1:169789/1‑104 (MQ=255)
                                                                                                                                                                                    tGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACAt                                                                                                                                                                                                                                                    <  2:169789/104‑1 (MQ=255)
                                                                                                                                                                                    tGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCa                                                                                                                                                                                                 >  1:195878/1‑155 (MQ=255)
                                                                                                                                                                                    tGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCa                                                                                                                                                                                                 <  2:195878/155‑1 (MQ=255)
                                                                                                                                                                                                              gTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTATACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGt                                                                                                                                                                                             >  2:134440/1‑133 (MQ=255)
                                                                                                                                                                                                              gTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGt                                                                                                                                                                                             <  1:134440/133‑1 (MQ=255)
                                                                                                                                                                                                                     tGGAAAACGTCCGCAACCTAATGTATCGCCTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGtt                                                                                                                                                               <  2:46447/156‑1 (MQ=255)
                                                                                                                                                                                                                        cAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACGACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTccc                                                                                                                                                            >  1:46447/1‑156 (MQ=255)
                                                                                                                                                                                                                          aacccccccaacaaAATGAACTGCCGATTACCGCCGGAGCAGAAAGCGGTTTTAAACGGGAAACATCAAAAACGAATGGGGCAGTTGCGTGACGGGACGCAATGGCAAAAAAGGTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACt                             <  2:382487/273‑1 (MQ=255)
                                                                                                                                                                                                                             gTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTa                                                             <  1:34249/246‑1 (MQ=255)
                                                                                                                                                                                                                                             aTCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAg                                                                                                                                                                                                                                     >  1:223854/1‑62 (MQ=255)
                                                                                                                                                                                                                                             aTCGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAg                                                                                                                                                                                                                                     <  2:223854/62‑1 (MQ=255)
                                                                                                                                                                                                                                               cGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGGGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGagtagt                                                                                                                                                                            <  1:64147/117‑1 (MQ=255)
                                                                                                                                                                                                                                               cGCCTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGagtagt                                                                                                                                                                            >  2:64147/1‑117 (MQ=255)
                                                                                                                                                                                                                                                 ccTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                                    >  1:104605/1‑219 (MQ=255)
                                                                                                                                                                                                                                                 ccTGTTAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                                    <  2:104605/219‑1 (MQ=255)
                                                                                                                                                                                                                                                    tgtAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTa  <  1:94263/280‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCa                                                                                                                         <  2:56955/162‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttAACGCCGGAGTAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCa                                                                                                                         >  1:56955/1‑162 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                            
AGAGCCATATGTTCGACGGCATAAGTTTAACCGAACATCAGCGTCAGCAGATGCGAGATCTTATGCAACAGGCCCGGCACGAACAGCCTCCTGTTAATGTTAGCGAACTGGAGACAATGCATCGCCTTGTCACCGCAGAAAATTTTGATGAAAACGCTGTGCGCGCACAGGCAGAAAAAATGGCGAATGAGCAAATTGCTCGTCAGGTTGAGATGGCAAAAGTCCGCAACCAAATGTATCGCCTGTTAACGCCGGAGCAGCAAGCGGTTTTAAACGAGAAACATCAACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGATAGCGGGAGTATTTA  >  NC_000913/4105944‑4106469

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