Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,285,716	T→C	V193A (GTC→GCC)	narJ →	nitrate reductase 1 molybdenum cofactor assembly chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,285,716	0	T	C	100.0%	113.7 / NA	33	V193A (GTC→GCC)	narJ	nitrate reductase 1 molybdenum cofactor assembly chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/18);  total (15/18)

CAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCGGAGGACAGCACTAATGC  >  NC_000913/1285583‑1285852                                                                                                                                      |                                                                                                                                         cAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGga                                                                                                                                    >  2:547142/1‑140 (MQ=255)    cGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGgaaga                                                                                                                                 >  2:615777/1‑140 (MQ=255)     gTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGgaagaa                                                                                                                                >  1:690918/1‑140 (MQ=255)         aaGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGc                                                                                                                            <  2:156030/140‑1 (MQ=255)         aaGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGc                                                                                                                            <  2:867597/140‑1 (MQ=255)         aaGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGc                                                                                                                            <  1:30495/140‑1 (MQ=255)                       tGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCttt                                                                                                              >  2:248570/1‑140 (MQ=255)                       tGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCttt                                                                                                              >  1:577196/1‑140 (MQ=255)                        gCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTg                                                                                                             <  1:1045518/140‑1 (MQ=255)                        gCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTg                                                                                                             <  1:269635/140‑1 (MQ=255)                                    gctgctgAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCt                                                                                                 >  1:766680/1‑140 (MQ=255)                                             aCTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGAtt                                                                                        <  1:820735/140‑1 (MQ=255)                                                    aaTACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAAt                                                                                 <  2:748241/140‑1 (MQ=255)                                                      tACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCa                                                                               <  1:547142/140‑1 (MQ=255)                                                                 aCAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAg                                                                    >  1:1092927/1‑140 (MQ=255)                                                                    gCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGt                                                                 <  2:577196/140‑1 (MQ=255)                                                                          aaGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGGGATTCAGCAATCACTGCGCATCAGCGTCGCttt                                                           >  2:988049/1‑140 (MQ=255)                                                                                  gAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGc                                                   >  2:957936/1‑140 (MQ=255)                                                                                  gAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGc                                                   >  2:811793/1‑140 (MQ=255)                                                                                   aaaaaaTTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGcc                                                  <  2:706033/140‑1 (MQ=255)                                                                                     aaaaTTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGt                                                <  1:615777/140‑1 (MQ=255)                                                                                             cgacgaAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGcgccgc                                        <  2:1124035/140‑1 (MQ=255)                                                                                                   aGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATc                                  <  2:573407/140‑1 (MQ=255)                                                                                                               tACGCCGCAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATcacca                      >  2:352698/1‑140 (MQ=255)                                                                                                                      cAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCggagg               >  2:1180254/1‑140 (MQ=255)                                                                                                                      cAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCggagg               >  1:540213/1‑140 (MQ=255)                                                                                                                      cAGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCggagg               >  1:295766/1‑140 (MQ=255)                                                                                                                       aGGCGCTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCggagga              <  2:295766/140‑1 (MQ=255)                                                                                                                          cgcTGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCGGAGGACAg           <  2:231667/140‑1 (MQ=255)                                                                                                                             tGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATcac                        >  2:182346/1‑124 (MQ=255)                                                                                                                             tGGATGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATcac                        <  1:182346/124‑1 (MQ=255)                                                                                                                              ggATGCAGACTGGGCAGACGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCGGAGGACAGCACt       >  2:1178928/1‑140 (MQ=255)                                                                                                                                aTGCTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCGGAGGACAGCACTaa     <  1:1180254/140‑1 (MQ=255)                                                                                                                                   cTGCCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCGGAGGACAGCACTAATGc  <  1:811793/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         CAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTGAATATCACCACCGGAGGACAGCACTAATGC  >  NC_000913/1285583‑1285852

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑