Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,951,161	A→T	intergenic (‑125/+70)	PP_3486 ← / ← PP_3487	cytochrome c/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,951,161	0	A	T	100.0%	94.8 / NA	29	intergenic (‑125/+70)	PP_3486/PP_3487	cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (19/10);  total (19/10)

TCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGG  >  NC_002947/3951024‑3951300                                                                                                                                          |                                                                                                                                            cccatcttGCGATGCCAAGCACGGAGGTAAGACATATCGCAACGCGTGTGGGAGGGGCCTTGCGTCGCGAATGGGGTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGGTGAAAGAGGGGTGCGCTGCGCGGTCCTTTCGCGACg                                                                                                                                     <  4:7208/140‑1 (MQ=25)    taGTTGCCATGCCAAGCACGCTGGAAAGCCATATCCCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAATAGGGGACCGCTGCGCGCTCCGTTCGCGACGCAAg                                                                                                                                 <  3:158515/146‑1 (MQ=38)     aGTTGCCATGCCAAGCACGCAGGTAAGCCATAGCGCAAGCCCTGTGGCAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAgg                                                                                                                                <  5:77997/147‑1 (MQ=255)          ccatgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGAGGCAGCTCAAGCCCAGAAACACGGGAGCGCTGCCCCCTCCTTTAGCGACGCAAgaacgga                                                                                                                           >  3:7208/1‑141 (MQ=37)                                                          ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCGCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCACTTCCCCAGGCGGCCCATg                                                                            >  5:9482/1‑146 (MQ=255)                                                                gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCa                                                                     >  4:234479/1‑147 (MQ=255)                                                                 tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                    >  8:111691/1‑147 (MQ=255)                                                                             tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc                                                        >  2:77745/1‑147 (MQ=255)                                                                                   gcagcCCCGGCAGTGGCAGCTCAAGCGCTTAACGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                                  <  2:54092/147‑1 (MQ=39)                                                                                     agcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCCCAGAGGCCAAATCACCATGCGGCCCATGATACTAACCGACCGGACAGCGGCtagt                                                 >  8:108714/1‑143 (MQ=38)                                                                                       ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCa                                              >  3:105554/1‑147 (MQ=255)                                                                                          cGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACa                                           >  3:254351/1‑147 (MQ=255)                                                                                           ggctggggCAGCACATGCGCATAACGCGGGGCCCGGGGCGGGCCCCTTGCGCCACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAg                                          <  5:96798/141‑1 (MQ=255)                                                                                             cAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTACTTTCGCCAAGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCACGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCt                                        >  1:72413/1‑147 (MQ=255)                                                                                              aGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCTGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTg                                       <  7:108714/147‑1 (MQ=255)                                                                                                 ggCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGCCCGGCCAGCGGCTGGGGCAACAGCTGCAt                                    >  6:243722/1‑147 (MQ=255)                                                                                                     gCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgc                                <  3:234479/147‑1 (MQ=255)                                                                                                       acAAGCGCACAACCCGGGGCGCGCGGCCGGCTCCTTGCGCCACGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTg                              <  7:231782/146‑1 (MQ=255)                                                                                                         aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTgg                             >  6:77746/1‑146 (MQ=255)                                                                                                                  aaGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgccg                   >  2:152482/1‑147 (MQ=255)                                                                                                                     agGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTg                >  5:224292/1‑147 (MQ=255)                                                                                                                     agGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCACGCGGCCCATGATACTCACCGACCGGCCAGGCGCTGGGTCAACAGCTGCATGCGCTGGGCAAAcgccgcag                >  5:224295/1‑145 (MQ=39)                                                                                                                          aGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCgggg                                                 <  4:162732/109‑1 (MQ=255)                                                                                                                          aGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTggg                                                 >  3:162732/1‑109 (MQ=255)                                                                                                                          aGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGt           >  6:43909/1‑147 (MQ=255)                                                                                                                              cTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCAt       >  6:1233/1‑147 (MQ=255)                                                                                                                              cTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCAt       >  3:61977/1‑147 (MQ=255)                                                                                                                                gcgcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATAACCCAGGCGGCCCATGATAAACACCCCCCCGCCAGAGGCTGGGTGACCAGCTGCAAGCGCTGGGCAATCGCCGCTGGCGGTGCATcc     >  1:213347/1‑147 (MQ=255)                                                                                                                                   cgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCggg  <  1:77745/147‑1 (MQ=255)                                                                                                                                          |                                                                                                                                            TCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGG  >  NC_002947/3951024‑3951300

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