Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,163,286	(A)5→4	intergenic (‑16/‑353)	trmA ← / → btuB	tRNA m(5)U54 methyltransferase/cobalamin/cobinamide outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,163,282	0	A	.	100.0%	150.6 / NA	36	intergenic (‑12/‑357)	trmA/btuB	tRNA m(5)U54 methyltransferase/cobalamin/cobinamide outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (23/13);  total (23/13)

ACCACTTTTTCGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTGAAAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGC  >  NC_000913/4163212‑4163356                                                                       |                                                                           aCCACTTTTTCGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGt                                                                 <  2:846328/81‑1 (MQ=255)      tttttCGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTAcgcg                                                            >  1:846394/1‑81 (MQ=255)        tttCGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACgcggc                                                          >  1:846329/1‑81 (MQ=255)          tCGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGa                                                        >  1:846361/1‑81 (MQ=255)           cGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGa                                                        <  2:846330/80‑1 (MQ=255)           cGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGAt                                                       >  2:846366/1‑81 (MQ=255)            ggCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGAtt                                                      >  1:846385/1‑81 (MQ=255)            ggCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGAtt                                                      <  2:846331/81‑1 (MQ=255)               tAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTg                                                     >  1:846390/1‑79 (MQ=255)                aaCTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAg                                                  >  2:846373/1‑81 (MQ=255)                  cTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAgg                                                 >  1:846371/1‑80 (MQ=255)                     cgcTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGAtt                                             <  2:846332/81‑1 (MQ=255)                       cTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGc                                           >  1:846384/1‑81 (MQ=255)                         tCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCt                                          >  2:846381/1‑80 (MQ=255)                          cATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCtca                                        >  2:846377/1‑81 (MQ=255)                           aTACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCtcat                                       >  2:846374/1‑81 (MQ=255)                           aTACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCtcat                                       <  1:846333/81‑1 (MQ=255)                                gttctgttGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGAt                                  >  2:846378/1‑81 (MQ=255)                                   ctgttGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTc                               >  2:846382/1‑81 (MQ=255)                                      ttGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTcc                              <  2:846336/79‑1 (MQ=255)                                      ttGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTcc                              <  1:846335/79‑1 (MQ=255)                                      ttGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAg                            <  2:846334/81‑1 (MQ=255)                                       tGGAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTc                               >  2:846364/1‑77 (MQ=255)                                         gAAGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGAt                          <  1:846337/80‑1 (MQ=255)                                          aaGGTGTTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCt                        >  2:846380/1‑81 (MQ=255)                                             gtgtTCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTtga                     >  1:846379/1‑81 (MQ=255)                                                ttCGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTtgatg                   <  1:846338/80‑1 (MQ=255)                                                  cGGGGGTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGaa                 <  1:846339/80‑1 (MQ=255)                                                   gggggTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAAtt               >  2:846403/1‑81 (MQ=255)                                                   gggggTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAAtt               <  2:846340/81‑1 (MQ=255)                                                   gggggTCATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAAGGATTGCTCATCAGATGTCCAGATCTTGATGAAtt               >  2:846386/1‑81 (MQ=255)                                                         cATGATGTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTAtt         >  2:846393/1‑81 (MQ=255)                                                          atgatgTATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTAttt        >  1:846388/1‑81 (MQ=255)                                                                tATCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAg   <  1:846341/80‑1 (MQ=255)                                                                 aTCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGc  >  2:846376/1‑80 (MQ=255)                                                                 aTCCTG‑AAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGc  >  2:846401/1‑80 (MQ=255)                                                                       |                                                                           ACCACTTTTTCGGCTAACTGCGCTTCATACTGTTCTGTTGGAAGGTGTTCGGGGGTCATGATGTATCCTGAAAAATTAAAGTACGCGGCGATTGTAGGGATTGCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGC  >  NC_000913/4163212‑4163356

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