Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,193 T→A 20.9% W19R (TGG→AGG)  narJ → molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,285,1930TA20.9% 103.7 / 22.6 43W19R (TGG→AGG) narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base A (5/4);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGA  >  NC_000913/1285047‑1285329
                                                                                                                                                  |                                                                                                                                         
gATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                                                                                                          <  1:658601/148‑1 (MQ=255)
 aTGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                                                                                                         <  4:953960/148‑1 (MQ=255)
           cGGTTCATATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgc                                                                                                                                        <  2:852697/139‑1 (MQ=255)
           cGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGGAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGACTTATGgc                                                                                                                                        >  1:852697/1‑139 (MQ=255)
                        aaaTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGgcag                                                                                                                                      <  3:511951/128‑1 (MQ=255)
                        aaaTTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGgcag                                                                                                                                      >  4:511951/1‑128 (MQ=255)
                            tCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGgcagc                                                                                                                                     >  4:892583/1‑125 (MQ=255)
                            tCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGgcagc                                                                                                                                     <  3:892583/125‑1 (MQ=255)
                               aTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCACCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           <  2:739463/116‑1 (MQ=255)
                               aTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCACCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           >  1:739463/1‑116 (MQ=255)
                                        aCAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCg                                                                                                  >  4:20995/1‑148 (MQ=255)
                                         cAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGa                                                                                                 >  4:794859/1‑148 (MQ=255)
                                                   aTCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGcc                                                                                       >  4:1036419/1‑148 (MQ=255)
                                                    tCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                   >  2:113497/1‑103 (MQ=255)
                                                    tCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                   <  3:782894/103‑1 (MQ=255)
                                                    tCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCACCCATGATCGATCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                   <  1:113497/103‑1 (MQ=255)
                                                    tCGATACCATCGAGGTGACCAGCAAAGCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                   >  4:782894/1‑103 (MQ=255)
                                                      gATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCaaa                                                                                    <  3:794859/148‑1 (MQ=255)
                                                         gCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                 <  3:20995/148‑1 (MQ=255)
                                                                      ccAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTg                                                                                                                   >  1:1047312/1‑101 (MQ=255)
                                                                      ccAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTg                                                                                                                   <  2:1047312/101‑1 (MQ=255)
                                                                        aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCt                                                                 >  2:493365/1‑148 (MQ=255)
                                                                         gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTg                                                                <  1:493365/148‑1 (MQ=255)
                                                                              aCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           >  3:186404/1‑69 (MQ=255)
                                                                              aCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           <  4:186404/69‑1 (MQ=255)
                                                                                 gAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtttt                                                        <  3:1036419/148‑1 (MQ=255)
                                                                                                         tGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgcag                                                                                                                                      >  1:429800/1‑47 (MQ=255)
                                                                                                         tGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgcag                                                                                                                                      <  2:429800/47‑1 (MQ=255)
                                                                                                                  cgcgTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                 <  1:688932/91‑1 (MQ=255)
                                                                                                                  cgcgTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                 >  2:688932/1‑91 (MQ=255)
                                                                                                                       ctcCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTg                                                                                         >  1:1041767/1‑78 (MQ=255)
                                                                                                                       ctcCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTg                                                                                         <  2:1041767/78‑1 (MQ=255)
                                                                                                                            tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtt                                                          >  2:425247/1‑103 (MQ=255)
                                                                                                                            tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCAtt                                                          <  1:425247/103‑1 (MQ=255)
                                                                                                                               atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGa                                               >  3:761129/1‑111 (MQ=255)
                                                                                                                               atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGa                                               <  4:761129/111‑1 (MQ=255)
                                                                                                                                tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                           <  4:1061/98‑1 (MQ=255)
                                                                                                                                tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                           >  3:1061/1‑98 (MQ=255)
                                                                                                                                       aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa  >  3:1015131/1‑148 (MQ=255)
                                                                                                                                            gcCTTATGGCAGCATCAAGAAGAGATGTTTGAGGCGATTGCCGCGTCGAAATATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATgg                                 >  4:717623/1‑112 (MQ=255)
                                                                                                                                            gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATgg                                 <  3:717623/112‑1 (MQ=255)
                                                                                                                                            gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAtc                                                          >  3:380005/1‑86 (MQ=255)
                                                                                                                                            gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAtc                                                          <  4:380005/87‑2 (MQ=255)
                                                                                                                                              cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt                          <  4:108979/117‑1 (MQ=255)
                                                                                                                                              cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt                          >  3:108979/1‑117 (MQ=255)
                                                                                                                                               ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                           >  1:917031/1‑83 (MQ=255)
                                                                                                                                               ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAt                                                           <  2:917031/83‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
GATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGA  >  NC_000913/1285047‑1285329

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: