Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,662,350 C→T P599S (CCC→TCC)  ynfF → S‑ and N‑oxide reductase, A subunit, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,662,3500CT97.5% 125.6 / ‑3.9 40P599S (CCC→TCC) ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/20);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

AAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662212‑1662493
                                                                                                                                          |                                                                                                                                               
aaaCCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                         >  2:123274/1‑147 (MQ=255)
 aaCCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      <  1:63208/149‑1 (MQ=255)
    ccATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      >  4:48425/1‑146 (MQ=255)
           cTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGaaa                                                                                                                             <  4:18899/148‑1 (MQ=255)
                    gAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGa                                                                                                                   >  4:99276/1‑149 (MQ=255)
                     aGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGAt                                                                                                                  <  1:123274/149‑1 (MQ=255)
                           gTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTaa                                                                                                            <  3:48425/149‑1 (MQ=255)
                                cAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                         <  4:62377/115‑1 (MQ=255)
                                       ttAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGccc                                                                                                >  7:142067/1‑149 (MQ=255)
                                           gACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaa                                                                                            <  8:142067/149‑1 (MQ=255)
                                                  cGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACt                                                                                     <  8:78026/149‑1 (MQ=255)
                                                       aTCAAACCTTAACTGAAGGTCGCAGTCAGCATGAATGGATCAAAAATCTCCATGCGAAAACGAAGGAACGTAACCCAGAGATCTCCCACTAAGATGAGAAGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACgac                                                                                >  8:18528/1‑147 (MQ=255)
                                                        tCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaagaa                                                                                         <  1:141515/139‑1 (MQ=255)
                                                        tCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaagaa                                                                                         >  2:141515/1‑139 (MQ=255)
                                                          aaaCCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCt                                                                             >  2:162943/1‑149 (MQ=255)
                                                           aaCCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCtt                                                                            >  5:15214/1‑149 (MQ=255)
                                                               tttACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCcg                                                                        >  6:9168/1‑149 (MQ=255)
                                                                         gTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCg                                                                                                                      >  2:52148/1‑93 (MQ=255)
                                                                             cAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTg                                                             >  4:1324/1‑146 (MQ=255)
                                                                              aGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGATGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGa                                                         <  3:99276/149‑1 (MQ=255)
                                                                              aGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGaa                                                           >  2:159779/1‑147 (MQ=255)
                                                                              aGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGaa                                                           <  1:159779/147‑1 (MQ=255)
                                                                                     aTGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCcac                                                    >  3:51747/1‑147 (MQ=255)
                                                                                     aTGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCcac                                                    <  4:51747/147‑1 (MQ=255)
                                                                                      tGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCTACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCGTTCCGCGCATTCCGTGAAGATCCACAgg                                                 >  7:109645/1‑149 (MQ=255)
                                                                                          tGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAgg                                                 <  4:11144/145‑1 (MQ=255)
                                                                                          tGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAgg                                                 >  3:11144/1‑145 (MQ=255)
                                                                                                aaaTATCTCCATGCTAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGtt                                       <  3:1324/149‑1 (MQ=255)
                                                                                                aaaTATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGtt                                       <  5:9168/149‑1 (MQ=255)
                                                                                                          aTGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCt                             <  1:162943/149‑1 (MQ=255)
                                                                                                               aaaaCGAAGGAACGTAACCCTGAGATGTCCGACTAAGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGTAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCggg                        <  7:18528/149‑1 (MQ=255)
                                                                                                                 aaCGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCgggg                       >  7:23305/1‑148 (MQ=255)
                                                                                                                    gAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGaaaa                   >  4:109618/1‑149 (MQ=255)
                                                                                                                      aGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                 <  3:109618/149‑1 (MQ=255)
                                                                                                                         aaCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGCGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTCAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGa               >  4:43420/1‑148 (MQ=255)
                                                                                                                             tAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTcc                                                                <  4:94742/95‑1 (MQ=255)
                                                                                                                                  cTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGTAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAt            <  2:152301/142‑1 (MQ=255)
                                                                                                                                  cTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAt            >  1:152301/1‑142 (MQ=255)
                                                                                                                                   tGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCt     >  5:76238/1‑148 (MQ=255)
                                                                                                                                     agaTGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGaa  <  8:23305/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
AAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662212‑1662493

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