Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	572,378	A→G	Y4C (TAC→TGC)	ylcH →	DLP12 prophage; uncharacterized protein YlcH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	572,378	0	A	G	100.0%	170.6 / NA	48	Y4C (TAC→TGC)	ylcH	DLP12 prophage; uncharacterized protein YlcH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/25);  total (23/25)

GTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAAC  >  NC_000913/572248‑572507                                                                                                                                   |                                                                                                                                  gTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCaaa                                                                                                                          <  1:288188/140‑1 (MQ=255)       gtgtTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAActctct                                                                                                                    <  1:378244/140‑1 (MQ=255)        tgtTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAAcgctctg                                                                                                                   <  2:378789/140‑2 (MQ=255)            ttAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTcc                                                                                                               >  1:1076498/1‑140 (MQ=255)            ttAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTcc                                                                                                               >  1:1101301/1‑140 (MQ=255)                  agGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTcac                                                                                                         >  2:659957/1‑140 (MQ=255)                     aGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGa                                                                                                      <  1:778229/140‑1 (MQ=255)                     aGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGa                                                                                                      <  1:956921/140‑1 (MQ=255)                                cGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGc                                                                                           <  2:872064/140‑1 (MQ=255)                                    tttgattttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCt                                                                                       >  1:75100/1‑140 (MQ=255)                                      tgattttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGc                                                                                           <  2:8988/134‑1 (MQ=255)                                      tgattttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGc                                                                                           >  1:8988/1‑134 (MQ=255)                                      tgatcttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGc                                                                                     <  2:1076498/140‑1 (MQ=255)                                         ttttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAAt                                                                                  <  1:1151563/140‑1 (MQ=255)                                                      acacGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAg                                                                     <  1:929301/140‑1 (MQ=255)                                                            cATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTg                                                               <  1:659957/140‑1 (MQ=255)                                                              tAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGAc                                                             <  2:899276/140‑1 (MQ=255)                                                               aaTTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGAcc                                                            >  2:1122357/1‑140 (MQ=255)                                                                    tGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGAcc                                                            <  1:696681/135‑1 (MQ=255)                                                                    tGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGAcc                                                            >  2:696681/1‑135 (MQ=255)                                                                        aCCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAg                                                   >  2:16585/1‑140 (MQ=255)                                                                        aCCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAg                                                   >  1:469807/1‑140 (MQ=255)                                                                         ccGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAgg                                                  >  1:647861/1‑140 (MQ=255)                                                                         ccGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAgg                                                  >  2:265654/1‑140 (MQ=255)                                                                          cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 >  1:1108023/1‑140 (MQ=255)                                                                          cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 <  1:265654/140‑1 (MQ=255)                                                                          cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 >  2:27045/1‑140 (MQ=255)                                                                          cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 <  1:1122316/140‑1 (MQ=255)                                                                           ggAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAgggg                                                <  1:584384/140‑1 (MQ=255)                                                                            gAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGa                                               >  1:607146/1‑140 (MQ=255)                                                                                   aacaacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 >  1:111290/1‑131 (MQ=255)                                                                                   aacaacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 <  2:111290/131‑1 (MQ=255)                                                                                   aacaacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGa                                        >  1:545447/1‑140 (MQ=255)                                                                                     caacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAAc                                      <  1:27045/140‑1 (MQ=255)                                                                                        cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 <  1:426635/126‑1 (MQ=255)                                                                                        cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                                 >  2:426635/1‑126 (MQ=255)                                                                                        cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGaa                                       <  2:919892/136‑1 (MQ=255)                                                                                        cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGaa                                       >  1:919892/1‑136 (MQ=255)                                                                                                     tGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                      <  2:545447/140‑1 (MQ=255)                                                                                                     tGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                      <  1:651483/140‑1 (MQ=255)                                                                                                     tGCGCTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                      <  1:1122357/140‑1 (MQ=255)                                                                                                       cgcTAAACGGGGACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATc                    >  1:180631/1‑140 (MQ=255)                                                                                                                gggACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACAt           <  2:1101301/140‑1 (MQ=255)                                                                                                                 ggACGTTTATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATc          >  2:894406/1‑140 (MQ=255)                                                                                                                      tttATGCGCACATGCAATCCAAACTCTCTTCTCCCTTTACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACATCTCACAAGACGGCATCAAATTACATCGCGGc     >  2:1033349/1‑140 (MQ=255)                                                                                                                      tttATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGc     >  1:841740/1‑140 (MQ=255)                                                                                                                      tttATGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGc     <  1:1033349/140‑1 (MQ=255)                                                                                                                         aTGCGCACATGCAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAAc  >  1:1133745/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                  GTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAAC  >  NC_000913/572248‑572507

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