Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,497,014	C→T	D4D (GAC→GAT)	nirC →	nitrite transporter NirC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,497,014	0	C	T	100.0%	159.6 / NA	46	D4D (GAC→GAT)	nirC	nitrite transporter NirC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/25);  total (21/25)

TGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCG  >  NC_000913/3496878‑3497147                                                                                                                                         |                                                                                                                                      tGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACt                                                                                                                                    <  2:464916/140‑1 (MQ=255) tGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACt                                                                                                                                    <  1:997527/140‑1 (MQ=255) tGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACt                                                                                                                                    <  1:286572/140‑1 (MQ=255)     ttAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATta                                                                                                                                >  2:977206/1‑140 (MQ=255)       aaCGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATtaat                                                                                                                              >  1:213525/1‑140 (MQ=255)                     tGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTaa                                                                                                                <  1:370748/140‑1 (MQ=255)                          ccccTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTg                                                                                                           >  1:899442/1‑140 (MQ=255)                            ccTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTgcg                                                                                                         >  1:411583/1‑140 (MQ=255)                                tttGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCa                                                                                                     <  1:190010/140‑1 (MQ=255)                                                   atCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCaaa                                                                                  >  2:576155/1‑140 (MQ=255)                                                   atCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCaaa                                                                                  >  2:62801/1‑140 (MQ=255)                                                   atCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCaaa                                                                                  >  2:1003103/1‑140 (MQ=255)                                                   atCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCaaa                                                                                  >  2:1027803/1‑140 (MQ=255)                                                   atCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCaaa                                                                                  >  2:356272/1‑140 (MQ=255)                                                   atCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCaaa                                                                                  >  1:468639/1‑140 (MQ=255)                                                      ttCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATaa                                                                               <  2:1008440/140‑1 (MQ=255)                                                             tatTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTc                                                                        <  2:1132746/140‑1 (MQ=255)                                                              atTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCg                                                                       <  1:364501/140‑1 (MQ=255)                                                                     gACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                                >  1:389719/1‑140 (MQ=255)                                                                     gACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                                >  1:662017/1‑140 (MQ=255)                                                                     gACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                                >  1:289839/1‑140 (MQ=255)                                                                       cATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATa                                                                                >  2:13065/1‑122 (MQ=255)                                                                       cATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATa                                                                                <  1:13065/122‑1 (MQ=255)                                                                       cATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTggg                                                              >  1:982618/1‑140 (MQ=255)                                                                        aTTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGt                                                             <  1:576155/140‑1 (MQ=255)                                                                          tttttAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCa                                                           >  1:1188262/1‑140 (MQ=255)                                                                           ttttAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAg                                                          <  1:356272/140‑1 (MQ=255)                                                                                    ttttGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCAt                                                 >  1:1152773/1‑140 (MQ=255)                                                                                     tttGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATg                                                >  1:238300/1‑140 (MQ=255)                                                                                       tGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggc                                              <  2:899442/140‑1 (MQ=255)                                                                                          tttATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcggg                                           <  1:977206/140‑1 (MQ=255)                                                                                            tATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgggcg                                         <  1:573668/140‑1 (MQ=255)                                                                                             aTTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATTCTATTAATAAGGGGGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTTGGTCAGCTCCGCCATGGCGGgcgc                                        <  1:1223115/140‑1 (MQ=255)                                                                                             aTTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGgcgc                                        >  2:531510/1‑140 (MQ=255)                                                                                                   tGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                                  <  2:213525/140‑1 (MQ=255)                                                                                                   tGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                                  <  2:944642/140‑1 (MQ=255)                                                                                                   tGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                                  <  2:982618/140‑1 (MQ=255)                                                                                                        ttATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATAGCGGGCGCGTATGTGGGTc                             <  2:411583/140‑1 (MQ=255)                                                                                                         tATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCt                            <  1:360317/140‑1 (MQ=255)                                                                                                         tATTTTTTAAAGGATAATCAAATGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCt                            <  2:289839/140‑1 (MQ=255)                                                                                                                              aTGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAc       <  2:678613/140‑1 (MQ=255)                                                                                                                              aTGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAc       <  2:19961/140‑1 (MQ=255)                                                                                                                               tGTTTACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACg      <  2:468639/140‑1 (MQ=255)                                                                                                                                 tttACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCt    <  1:1003103/140‑1 (MQ=255)                                                                                                                                   tACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCg  >  2:1072922/1‑140 (MQ=255)                                                                                                                                   tACAGATACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGAACATCCTGATTTTCACGCTCg  >  2:898594/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      TGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCG  >  NC_000913/3496878‑3497147

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: ‑