Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,283,686	(A)5→4	coding (83/1539 nt)	narH →	nitrate reductase A subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,283,682	0	A	.	100.0%	212.7 / NA	50	coding (79/1539 nt)	narH	nitrate reductase A subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (23/27);  total (23/27)

GAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAAC  >  NC_000913/1283550‑1283816                                                                                                                                     |                                                                                                                                       gaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaaacgtc                                                                                                                                 >  2:613992/1‑136 (MQ=255)     aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                                                                                                      <  2:503750/131‑1 (MQ=255)     aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                                                                                                      >  1:503750/1‑131 (MQ=255)      ttGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGAc                                                                                                                           <  2:424592/140‑1 (MQ=255)      ttGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGAc                                                                                                                           <  2:123562/140‑1 (MQ=255)         aCTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAg                                                                                                                        >  2:918958/1‑140 (MQ=255)         aCTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAg                                                                                                                        >  2:193738/1‑140 (MQ=255)              ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTg                                                                                                                   <  1:193738/140‑1 (MQ=255)              ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTg                                                                                                                   <  1:996586/140‑1 (MQ=255)              ttAGATGGCGAAGACAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTg                                                                                                                   <  2:318088/140‑1 (MQ=255)               tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGa                                                                                                                  >  2:652290/1‑140 (MQ=255)                   tGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGa                                                                                                                            <  2:432296/126‑1 (MQ=255)                   tGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGa                                                                                                                            >  1:432296/1‑126 (MQ=255)                    ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCg                                                                                                             <  2:956433/140‑1 (MQ=255)                    ggCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCg                                                                                                             <  1:652290/140‑1 (MQ=255)                       gAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTg                                                                                                                   >  1:345810/1‑131 (MQ=255)                       gAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTg                                                                                                                   <  2:345810/131‑1 (MQ=255)                        aaGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGa                                                                                                         <  1:821783/140‑1 (MQ=255)                               tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCg                                                                                                  >  1:788864/1‑140 (MQ=255)                               tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCg                                                                                                  <  1:91354/140‑1 (MQ=255)                               tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCg                                                                                                  <  1:558254/140‑1 (MQ=255)                               tGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGCCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAGTATGCg                                                                                                  <  2:1015935/140‑1 (MQ=255)                                         cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTATCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                                                        >  1:385570/1‑140 (MQ=255)                                         cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                                                        >  1:664009/1‑140 (MQ=255)                                         cAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                                                        <  1:613992/140‑1 (MQ=255)                                               agCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTgg                                                                                  >  1:1079859/1‑140 (MQ=255)                                               agCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTgg                                                                                  >  2:4392/1‑140 (MQ=255)                                                gCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                                 >  1:403233/1‑140 (MQ=255)                                                 cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                                                                                >  2:60039/1‑140 (MQ=255)                                                 cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                                                                                >  1:948622/1‑140 (MQ=255)                                                 cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                                                                                >  1:14978/1‑140 (MQ=255)                                                 cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                                                                                >  1:272030/1‑140 (MQ=255)                                                        tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                                                                         <  1:153167/140‑1 (MQ=255)                                                        tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                                                                         <  2:948622/140‑1 (MQ=255)                                                        tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                                                                         <  1:610249/140‑1 (MQ=255)                                                            aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                                                     <  2:78903/140‑1 (MQ=255)                                                                 gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGTAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAgg                                                                >  2:1020507/1‑140 (MQ=255)                                                                     acaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                                                            <  2:619285/140‑1 (MQ=255)                                                                         gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCg                                                        <  2:664009/140‑1 (MQ=255)                                                                         gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCCTCCCg                                                        <  2:385570/140‑1 (MQ=255)                                                                          tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCga                                                       <  2:1079859/140‑1 (MQ=255)                                                                                  gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTgg                                               <  2:14978/140‑1 (MQ=255)                                                                                         aTCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAAcc                                        <  1:1078065/140‑1 (MQ=255)                                                                                             cGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGa                                    >  1:58164/1‑140 (MQ=255)                                                                                             cGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGa                                    >  1:939157/1‑140 (MQ=255)                                                                                              gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaa                                   <  2:637265/140‑1 (MQ=255)                                                                                                  aGTGCATCGGCTGCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAAt                               >  2:144395/1‑140 (MQ=255)                                                                                                              gCCACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCt                   >  2:585360/1‑140 (MQ=255)                                                                                                               ccACACCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTg                  <  2:788864/140‑1 (MQ=255)                                                                                                                 acacCTGTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGa                >  2:246340/1‑140 (MQ=255)                                                                                                                       gTTCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTa          <  2:418757/140‑1 (MQ=255)                                                                                                                         tCAGTTACCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTaaa        >  2:877933/1‑140 (MQ=255)                                                                                                                               aCCTGT‑AAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAAc  >  2:443472/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                       GAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGATCCGTAAAATCAAC  >  NC_000913/1283550‑1283816

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