Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,493,804	T→C	intergenic (+55/‑207)	tsgA → / → nirB	putative transporter TsgA/nitrite reductase catalytic subunit NirB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,493,804	0	T	C	100.0%	100.6 / NA	29	intergenic (+55/‑207)	tsgA/nirB	putative transporter TsgA/nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/21);  total (8/21)

CATAACACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCTACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTG  >  NC_000913/3493726‑3493879                                                                               |                                                                            cATAACACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCAc                                                                           <  1:667356/81‑1 (MQ=255) cATAACACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCAc                                                                           <  1:667355/81‑1 (MQ=255) cATAACACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCCTCCCCCGCCAc                                                                           <  2:667354/81‑1 (MQ=255)  aTAACACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACg                                                                          >  1:667383/1‑81 (MQ=255)    aaCACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTg                                                                        >  1:667385/1‑81 (MQ=255)        cccTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATc                                                                    <  1:667357/81‑1 (MQ=255)        cccTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATc                                                                    <  1:667358/81‑1 (MQ=255)          cTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTa                                                                  >  1:667394/1‑81 (MQ=255)               ctctCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCt                                                             <  1:667359/81‑1 (MQ=255)                        aTTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCg                                                    <  2:667360/81‑1 (MQ=255)                           gCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCt                                                  <  2:667361/80‑1 (MQ=255)                            cTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTc                                                 <  2:667363/80‑1 (MQ=255)                            cTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTc                                                 <  1:667364/80‑1 (MQ=255)                            cTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCa                                                <  2:667362/81‑1 (MQ=255)                              cATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCAtt                                              >  2:667391/1‑81 (MQ=255)                               aTGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCAt                                               <  2:667365/79‑1 (MQ=255)                                 gCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTg                                           >  1:667389/1‑81 (MQ=255)                                  ccGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTg                                           <  1:667366/80‑1 (MQ=255)                                    ggACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTc                                         >  1:667395/1‑80 (MQ=255)                                                    gtcCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTCTTTTTTGCACAAACa                        >  1:667398/1‑81 (MQ=255)                                                    gtcCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAAc                         <  1:667367/80‑1 (MQ=255)                                                    gtcCGGCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAAc                         <  2:667368/80‑1 (MQ=255)                                                        ggCCTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGaa                    <  2:667369/81‑1 (MQ=255)                                                          ccTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAat                  <  2:667371/81‑1 (MQ=255)                                                          ccTTTTCCTCTCTTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAat                  <  1:667370/81‑1 (MQ=255)                                                                   ctctTCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAATATCAGACaa         >  2:667390/1‑81 (MQ=255)                                                                      ttCCCCCGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAATATCAGACaa         <  1:667372/78‑1 (MQ=255)                                                                        cccccGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAATATCAGACAATTc      <  1:667373/79‑1 (MQ=255)                                                                          cccGCCACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTg  <  2:667374/81‑1 (MQ=255)                                                                               |                                                                            CATAACACCCTGACCTCTCATTAATTGCTCATGCCGGACGGCACTATCGTCGTCCGGCCTTTTCCTCTCTTCCCCCGCTACGTGCATCTATTTCTATAAACCCGCTCATTTTGTCTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTG  >  NC_000913/3493726‑3493879

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