Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	287,977	A→G	M397V (ATG→GTG)	yagH →	CP4‑6 prophage, putative xylosidase/arabinosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	287,977	0	A	G	100.0%	109.9 / NA	33	M397V (ATG→GTG)	yagH	CP4‑6 prophage, putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/19);  total (14/19)

GCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAGTTATCCAGCTCGACC  >  NC_000913/287866‑288109                                                                                                                |                                                                                                                                     gcgcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagc                                                                                                           <  2:333805/139‑1 (MQ=255)  cgcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                                          <  2:254110/139‑1 (MQ=255)          ggCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgg                                                                                                  >  1:314705/1‑139 (MQ=255)                ttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                                            <  2:314053/139‑1 (MQ=255)                  aCGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAcc                                                                                          <  1:604077/139‑1 (MQ=255)                   cGGCTCTATGGCGACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGTACTTCCAGCAGAGCGCGGGGCTGACCt                                                                                         <  2:459448/139‑1 (MQ=255)                   cGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                                         <  2:622306/139‑1 (MQ=255)                   cGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                                         <  1:662503/139‑1 (MQ=255)                     gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                                           <  1:311474/135‑1 (MQ=255)                     gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                                           >  2:311474/1‑135 (MQ=255)                     gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                                                       >  1:40480/1‑139 (MQ=255)                         tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                                    >  2:644852/1‑106 (MQ=255)                         tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                                    <  1:644852/106‑1 (MQ=255)                                aCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAg                                                                            >  2:770921/1‑139 (MQ=255)                                   aCTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGTAGAGCGCGGGGCTGACCTGCTACTACAACAGCaa                                                                         >  1:353854/1‑139 (MQ=255)                                            aTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                                            <  1:490627/111‑1 (MQ=255)                                            aTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                                            >  2:490627/1‑111 (MQ=255)                                             ttCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGc                                                               <  1:601486/139‑1 (MQ=255)                                               cGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTa                                                             <  1:211552/139‑1 (MQ=255)                                                                aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                                          <  1:315030/77‑1 (MQ=255)                                                                aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                                          >  2:315030/1‑77 (MQ=255)                                                                aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGGTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACg                                            >  2:364369/1‑139 (MQ=255)                                                                     ggCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag                                       <  1:364369/139‑1 (MQ=255)                                                                        gcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa                                    <  2:302087/139‑1 (MQ=255)                                                                           ccgcTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                                         >  2:38991/1‑83 (MQ=255)                                                                           ccgcTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                                         <  1:38991/83‑1 (MQ=255)                                                                                          cGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGTAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCaaa                  <  2:353854/139‑1 (MQ=255)                                                                                           gCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAg                 <  1:301780/139‑1 (MQ=255)                                                                                           gCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAg                 >  1:131866/1‑139 (MQ=255)                                                                                                ccGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCTCTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAGTTATc            >  2:713267/1‑139 (MQ=255)                                                                                                         aaCGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTaca                                                                                <  2:106040/62‑1 (MQ=255)                                                                                                         aaCGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTaca                                                                                >  1:106040/1‑62 (MQ=255)                                                                                                          aCGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAGTTATCCAGCTCGAcc  >  1:618519/1‑139 (MQ=255)                                                                                                                |                                                                                                                                     GCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAGTTATCCAGCTCGACC  >  NC_000913/287866‑288109

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