Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,548,360	C→T	D320D (GAC→GAT)	fdnG →	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,548,360	0	C	T	100.0%	164.7 / NA	51	D320D (GAC→GAT)	fdnG	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/27);  total (24/27)

GTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGG  >  NC_000913/1548233‑1548483                                                                                                                                |                                                                                                                            gTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAg                                                                                                                  <  2:1103424/139‑1 (MQ=255)           tACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacgac                                                                                                        <  2:488800/139‑1 (MQ=255)            aCGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacgacg                                                                                                       >  1:342249/1‑139 (MQ=255)                 ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                  >  1:1020513/1‑139 (MQ=255)                 ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                  >  2:407029/1‑139 (MQ=255)                 ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                  >  2:275581/1‑139 (MQ=255)                 ccTGTCTGACGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTACGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                  <  2:913234/139‑1 (MQ=255)                   tGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACGAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaaa                                                                                                >  1:941866/1‑139 (MQ=255)                       tGGCGTTTTGTGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                            <  1:275581/139‑1 (MQ=255)                          cGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCac                                                                                         <  2:941866/139‑2 (MQ=255)                                   cTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                <  2:1034102/139‑1 (MQ=255)                                   cTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                <  2:1020513/139‑1 (MQ=255)                                      ccTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGt                                                                                                                          <  1:1122461/94‑1 (MQ=255)                                      ccTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGt                                                                                                                          >  2:1122461/1‑94 (MQ=255)                                         gATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            <  2:606455/105‑1 (MQ=255)                                         gATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            >  1:606455/1‑105 (MQ=255)                                              aaaaCAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGaa                                                                     >  1:12811/1‑139 (MQ=255)                                                aacaatAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACt                                                                   >  2:452417/1‑139 (MQ=255)                                                aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacg                                                                                                          >  1:113212/1‑100 (MQ=255)                                                aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacg                                                                                                          <  2:113212/100‑1 (MQ=255)                                                    acaaAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCa                                                               >  1:343399/1‑139 (MQ=255)                                                         aTCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            <  1:203942/89‑1 (MQ=255)                                                         aTCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            >  2:203942/1‑89 (MQ=255)                                                              cGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaa                                                     <  1:452417/139‑1 (MQ=255)                                                                                ttACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            >  1:1376567/1‑66 (MQ=255)                                                                                ttACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            <  2:1376567/66‑1 (MQ=255)                                                                                        aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                         >  2:531479/1‑77 (MQ=255)                                                                                        aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                         <  1:531479/77‑1 (MQ=255)                                                                                           ccAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACt                        >  1:1056336/1‑139 (MQ=255)                                                                                              gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCa                                                                                          >  2:751457/1‑70 (MQ=255)                                                                                              gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCa                                                                                          <  1:751457/70‑1 (MQ=255)                                                                                              gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATAc                                                                                      <  1:1441065/74‑1 (MQ=255)                                                                                              gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATAc                                                                                      >  2:1441065/1‑74 (MQ=255)                                                                                                       tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg            <  2:12811/139‑1 (MQ=255)                                                                                                       tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg            <  2:1056336/139‑1 (MQ=255)                                                                                                        gCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                >  2:1297880/1‑70 (MQ=255)                                                                                                        gCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                <  1:1297880/70‑1 (MQ=255)                                                                                                        gCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCgcgc           <  1:455582/139‑1 (MQ=255)                                                                                                           tgatgaTTTTGCTTTGGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcg                                    >  2:1008372/1‑111 (MQ=255)                                                                                                           tgatgaTTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            >  2:808382/1‑39 (MQ=255)                                                                                                           tgatgaTTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                            <  1:808382/39‑1 (MQ=255)                                                                                                           tgatgaTTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCCGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcg                                    <  1:1008372/111‑1 (MQ=255)                                                                                                             atgatTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACt                    >  1:38851/1‑125 (MQ=255)                                                                                                             atgatTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACt                    <  2:38851/125‑1 (MQ=255)                                                                                                                atTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                           <  2:781125/67‑1 (MQ=255)                                                                                                                atTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                           >  1:781125/1‑67 (MQ=255)                                                                                                                 ttttGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTgg  <  1:407029/139‑1 (MQ=255)                                                                                                                      cTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaa                                                                                                 <  1:357169/39‑1 (MQ=255)                                                                                                                      cTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaa                                                                                                 >  2:357169/1‑39 (MQ=255)                                                                                                                          cGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                         >  1:211221/1‑43 (MQ=255)                                                                                                                          cGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                         <  2:211221/43‑1 (MQ=255)                                                                                                                                |                                                                                                                            GTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGG  >  NC_000913/1548233‑1548483

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