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,3500CT96.0% 161.0 / ‑1.9 50P599S (CCC→TCC) ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (25/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

ATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAAC  >  NC_000913/1662218‑1662494
                                                                                                                                    |                                                                                                                                                
aTCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGgaga                                                                                                                                   >  3:69542/1‑148 (MQ=255)
aTCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGAt                                                                                                                                  >  1:167690/1‑149 (MQ=255)
            cTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCgg                                                                                                                      <  3:77722/149‑1 (MQ=255)
              gAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGa                                                                                                                    >  2:182567/1‑149 (MQ=255)
                 tGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCt                                                                                                                 <  3:30553/149‑1 (MQ=255)
                   aaGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGAt                                                                                                                   <  2:116372/145‑1 (MQ=255)
                   aaGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGAt                                                                                                                   >  1:116372/1‑145 (MQ=255)
                   aaGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCttt                                                                                                               <  1:182567/149‑1 (MQ=255)
                          cAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGaaa                                                                                                         >  8:119448/1‑148 (MQ=255)
                               gCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGc                                                                                                   <  7:119448/149‑1 (MQ=255)
                               gCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGc                                                                                                   >  5:173023/1‑149 (MQ=255)
                               gCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACACTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGc                                                                                                   >  2:22388/1‑149 (MQ=255)
                                   aGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCg                                                                                                >  4:31626/1‑148 (MQ=255)
                                     gACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaa                                                                                             <  6:173023/149‑1 (MQ=255)
                                     gACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGgaa                                                                                             <  5:107790/149‑1 (MQ=255)
                                          gACGTTTATCAAACCTTTACTGAAGGTCTCAGTCAGCATGAATGGATGAAATATCTCCATGCGAAAACGAAGGAAGGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAaca                                                                                        <  7:19092/149‑1 (MQ=255)
                                          gACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGa                                                                                                                                              >  7:118924/1‑95 (MQ=255)
                                          gACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGa                                                                                                                                              <  8:118924/95‑1 (MQ=255)
                                             gTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTa                                                                                     >  8:76106/1‑149 (MQ=255)
                                                 aTCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACTCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGt                                                                                  <  1:22388/148‑1 (MQ=255)
                                                 aTCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTc                                                                                 >  7:27441/1‑149 (MQ=255)
                                                  tCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCg                                                                                <  8:27441/149‑1 (MQ=255)
                                                         tttACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCcg                                                                         >  8:46174/1‑149 (MQ=255)
                                                             cTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGTAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCa                                                                     <  7:19110/149‑1 (MQ=255)
                                                               gAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGCCCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCAtt                                                                   >  6:18534/1‑149 (MQ=255)
                                                                   gTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGt                                                               >  3:29858/1‑149 (MQ=255)
                                                                      gCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTg                                                              >  4:153249/1‑147 (MQ=255)
                                                                       cAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                          <  4:112734/70‑1 (MQ=255)
                                                                       cAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTAc                                                                                                                                          >  3:112734/1‑70 (MQ=255)
                                                                            aGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATcc                                                       <  7:76106/148‑1 (MQ=255)
                                                                            aGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCca                                                      <  4:29858/149‑1 (MQ=255)
                                                                            aGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCca                                                      >  6:186283/1‑149 (MQ=255)
                                                                              cATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCcaca                                                    <  7:46174/149‑1 (MQ=255)
                                                                                        tCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAAc                                             <  1:49673/146‑1 (MQ=255)
                                                                                        tCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAAc                                             >  2:49673/1‑146 (MQ=255)
                                                                                            atatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                  <  3:33033/89‑1 (MQ=255)
                                                                                            atatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGcc                                                                                                  >  4:33033/1‑89 (MQ=255)
                                                                                              atCTCCATGCGAAAACGAATGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGAt                                                                                                                   <  5:63145/70‑1 (MQ=255)
                                                                                              atCTCCATGCGAAAACGAATGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGAt                                                                                                                   >  6:63145/1‑70 (MQ=255)
                                                                                                    aTGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCt                              <  3:31626/149‑1 (MQ=255)
                                                                                                        gAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAac                                                                                         >  3:123574/1‑86 (MQ=255)
                                                                                                        gAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAac                                                                                         <  4:123574/86‑1 (MQ=255)
                                                                                                        gAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCAGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCgg                          >  2:35290/1‑149 (MQ=255)
                                                                                                         aaaaCGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCggg                         >  1:202923/1‑149 (MQ=255)
                                                                                                                aGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                  <  2:167690/149‑1 (MQ=255)
                                                                                                                 ggAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCg                 <  5:186283/149‑1 (MQ=255)
                                                                                                                          cccTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACCGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTAtt        >  8:128104/1‑149 (MQ=255)
                                                                                                                           ccTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTc       <  3:153249/149‑1 (MQ=255)
                                                                                                                            cTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCt      >  7:38295/1‑149 (MQ=255)
                                                                                                                                 aTGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAAc  >  4:99270/1‑148 (MQ=255)
                                                                                                                                    |                                                                                                                                                
ATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAAC  >  NC_000913/1662218‑1662494

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: ‑