Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,786,070	C→T	intergenic (‑159/+116)	kbl ← / ← yibB	glycine C‑acetyltransferase/protein YibB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,786,070	0	C	T	100.0%	132.0 / NA	33	intergenic (‑159/+116)	kbl/yibB	glycine C‑acetyltransferase/protein YibB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/15);  total (18/15)

CAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTT  >  NZ_CP009273/3785938‑3786195                                                                                                                                     |                                                                                                                              cAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATTGTAGAAAAACTGTACAGTTTTGGTTGGGc                                                                                                                          <  2:536714/138‑1 (MQ=255) cAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGc                                                                                                                          <  1:550985/138‑1 (MQ=255)       atatTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAatga                                                                                                                   <  1:519844/139‑1 (MQ=255)            aCCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATgatgat                                                                                                               <  2:419960/138‑1 (MQ=255)            aCCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATgatgat                                                                                                               <  2:450902/138‑1 (MQ=255)              ccAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGAtgattga                                                                                                            <  2:154663/139‑1 (MQ=255)                              cgacgaGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGata                                                                                            >  2:226394/1‑139 (MQ=255)                                    gaTAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATaaa                                                                                        >  2:441936/1‑137 (MQ=255)                                                    ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATgatgat                                                                                                               >  2:548875/1‑98 (MQ=255)                                                    ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATgatgat                                                                                                               <  1:548875/98‑1 (MQ=255)                                                            aaTTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatcat                                                              >  1:1070/1‑139 (MQ=255)                                                                            cGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCa                                                                                                  >  2:131675/1‑87 (MQ=255)                                                                            cGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCa                                                                                                  <  1:131675/87‑1 (MQ=255)                                                                                      gATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCaaaa                                    >  2:58801/1‑139 (MQ=255)                                                                                       aTGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCaaaa                                    <  1:226394/138‑1 (MQ=255)                                                                                       aTGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAAc                                   >  1:289793/1‑139 (MQ=255)                                                                                             ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAcc                             >  1:349864/1‑139 (MQ=255)                                                                                              tcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTa                               <  2:550379/136‑1 (MQ=255)                                                                                                aaCACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAcc                             <  2:349864/136‑1 (MQ=255)                                                                                                aaCACCTACAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAcc                             <  2:96117/136‑1 (MQ=255)                                                                                                       aCAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTAt                                                     >  2:35671/1‑105 (MQ=255)                                                                                                       aCAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTAt                                                     <  1:35671/105‑1 (MQ=255)                                                                                                        cAGATAGTAGAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                  >  1:174858/1‑139 (MQ=255)                                                                                                             agtagAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa             >  1:71307/1‑139 (MQ=255)                                                                                                             agtagAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaaa             >  1:308812/1‑139 (MQ=255)                                                                                                               tagAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTa                                                                                                       <  1:49724/47‑1 (MQ=255)                                                                                                               tagAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTa                                                                                                       >  2:49724/1‑47 (MQ=255)                                                                                                               tagAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAcc           >  2:193555/1‑139 (MQ=255)                                                                                                               tagAAAAACTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAAcc           >  2:13003/1‑139 (MQ=255)                                                                                                                     aaCTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAcc                             <  2:227093/115‑1 (MQ=255)                                                                                                                     aaCTGTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTAcc                             >  1:227093/1‑115 (MQ=255)                                                                                                                         gTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt  >  1:623690/1‑138 (MQ=255)                                                                                                                         gTACAGTTTTGGTTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCtt  >  1:31671/1‑138 (MQ=255)                                                                                                                                     |                                                                                                                              CAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTT  >  NZ_CP009273/3785938‑3786195

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