Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,190 T→A 19.1% L18I (TTA→ATA)  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,1900TA19.1% 119.4 / 22.9 47L18I (TTA→ATA) narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (5/4);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGA  >  NC_000913/1285044‑1285329
                                                                                                                                                  |                                                                                                                                            
ggTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCtt                                                                                                                                             <  3:315497/148‑1 (MQ=255)
   gATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                                                                                                          <  1:658601/148‑1 (MQ=255)
    aTGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                                                                                                         <  4:953960/148‑1 (MQ=255)
     tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCTATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTa                                                                                                                                            <  4:649941/144‑1 (MQ=255)
     tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTa                                                                                                                                            >  3:649941/1‑144 (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                                                                                                                                           >  1:739463/1‑116 (MQ=255)
                                  aTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCACCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           <  2:739463/116‑1 (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)
                                                                           aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑t                                                                 >  2:493365/1‑148 (MQ=255)
                                                                            gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑Tg                                                                <  1:493365/148‑1 (MQ=255)
                                                                                 aCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           >  3:186404/1‑69 (MQ=255)
                                                                                 aCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                           <  4:186404/69‑1 (MQ=255)
                                                                                    gAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCAtttt                                                        <  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)
                                                                                                                               tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCAtt                                                          >  2:425247/1‑103 (MQ=255)
                                                                                                                               tGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCAtt                                                          <  1:425247/103‑1 (MQ=255)
                                                                                                                                  atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGa                                               >  3:761129/1‑111 (MQ=255)
                                                                                                                                  atatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGa                                               <  4:761129/111‑1 (MQ=255)
                                                                                                                                   tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt                                                           <  4:1061/98‑1 (MQ=255)
                                                                                                                                   tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt                                                           >  3:1061/1‑98 (MQ=255)
                                                                                                                                          aTGCTGCCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGa  >  3:1015131/1‑148 (MQ=255)
                                                                                                                                               gcCTTATGGCAGCATCAAGAAGAGATGTTTGAGGCGATTGCCGCGTCGAAATATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATgg                                 >  4:717623/1‑112 (MQ=255)
                                                                                                                                               gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATgg                                 <  3:717623/112‑1 (MQ=255)
                                                                                                                                               gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAtc                                                          <  4:380005/87‑2 (MQ=255)
                                                                                                                                               gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAtc                                                          >  3:380005/1‑86 (MQ=255)
                                                                                                                                                 cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt                          <  4:108979/117‑1 (MQ=255)
                                                                                                                                                 cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCt                          >  3:108979/1‑117 (MQ=255)
                                                                                                                                                  ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt                                                           <  2:917031/83‑1 (MQ=255)
                                                                                                                                                  ttATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCA‑GCGCATGGGCAt                                                           >  1:917031/1‑83 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGC‑TGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGA  >  NC_000913/1285044‑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: ‑