Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,236,848	(T)5→6	coding (394/1533 nt)	ycgB ←	PF04293 family protein YcgB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,236,843	1	.	T	100.0%	69.1 / NA	20	N133K (AAC→AAA)	ycgB	PF04293 family protein YcgB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (10/10);  total (10/10)

ACTCGGTAATATATTTACGGGCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTG‑TTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATGGTGTTCTCT  >  NC_000913/1236764‑1236915                                                                                 |                                                                         aCTCGGTAATATATTTACGGGCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTG‑t                                                                         >  1:242277/1‑81 (MQ=255)   tCGGTAATATATTTACGGGCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTG‑ttt                                                                       >  1:242278/1‑81 (MQ=255)            tatTTACGGGCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTgaagaa                                                               <  2:242257/81‑1 (MQ=255)                  cGGGCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGt                                                          >  1:242292/1‑80 (MQ=255)                     gCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATg                                                      <  1:242258/81‑1 (MQ=255)                          aaTCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCAt                                                 <  1:242259/81‑1 (MQ=255)                                aTAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAgg                                           >  1:242283/1‑81 (MQ=255)                                  aaTCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGGAGAAAGAGTTATGCCCATAGCAGGCa                                         >  2:242317/1‑81 (MQ=255)                                          aTCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCAt                                 >  1:242285/1‑81 (MQ=255)                                           tCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCAt                                 >  1:242276/1‑80 (MQ=255)                                               aCTGGCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCa                            <  2:242260/81‑1 (MQ=255)                                                  ggCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAgcg                         >  1:242281/1‑81 (MQ=255)                                                   gCGTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAgcgc                        >  1:242284/1‑81 (MQ=255)                                                     gTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCtt                      <  2:242261/81‑1 (MQ=255)                                                     gTCGGTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCtt                      >  2:242293/1‑81 (MQ=255)                                                        ggTCCAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCa                   >  1:242288/1‑81 (MQ=255)                                                            cAGCTACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGt               <  1:242262/81‑1 (MQ=255)                                                                tACGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATg           >  2:242294/1‑81 (MQ=255)                                                                  cGGAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATg           <  2:242263/79‑1 (MQ=255)                                                                    gAATAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATGgtgt       <  1:242264/81‑1 (MQ=255)                                                                     aaTAAGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATGGTGtt      <  2:242265/81‑1 (MQ=255)                                                                         aGTAATTGTTTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATGGTGTTctct  <  1:242266/81‑1 (MQ=255)                                                                                 |                                                                         ACTCGGTAATATATTTACGGGCAAAAATCAGATAATCGACAATCGAACTGGCGTCGGTCCAGCTACGGAATAAGTAATTG‑TTTTTGAAGAAAGAGTTATGCCCATAGCAGGCATGAGCCATCACCAGCGCTTGCATGGTAATGGTGTTCTCT  >  NC_000913/1236764‑1236915

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