Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	577,085	G→A	intergenic (‑260/‑313)	nmpC ← / → essD	DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage/DLP12 prophage; putative phage lysis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	577,085	0	G	A	100.0%	203.6 / NA	52	intergenic (‑260/‑313)	nmpC/essD	DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage/DLP12 prophage; putative phage lysis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/24);  total (28/24)

GTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGGCTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTG  >  NC_000913/576937‑577232                                                                                                                                                     |                                                                                                                                                    gtgAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACt                                                                                                                                                   >  2:577508/1‑151 (MQ=255)   gAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTTATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTaa                                                                                                                                                 <  1:126810/151‑1 (MQ=255)        aGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAAt                                                                                                                                             >  2:11857/1‑150 (MQ=255)          aTCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATaa                                                                                                                                           >  2:4878/1‑150 (MQ=255)          aTCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATaa                                                                                                                                           <  1:4878/150‑1 (MQ=255)           tctcTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAac                                                                                                                                          >  1:424532/1‑150 (MQ=255)             tctATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACt                                                                                                                                       >  2:647813/1‑151 (MQ=255)             tctATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACt                                                                                                                                       >  1:17738/1‑151 (MQ=255)                aTTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTg                                                                                                                                      <  1:546433/149‑1 (MQ=255)                          gAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAattatt                                                                                                                             <  1:560526/148‑1 (MQ=255)                          gAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAattat                                                                                                                              <  2:75166/147‑1 (MQ=255)                           aaCTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTAtttttt                                                                                                                         <  2:102891/151‑1 (MQ=255)                                  ttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAACCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTAt                                                                                                                   <  2:13874/150‑1 (MQ=255)                                      aTCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTAttt                                                                                                                 <  2:333814/148‑1 (MQ=255)                                           aTCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgttg                                                                                                             >  2:404981/1‑147 (MQ=255)                                                aTTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGc                                                                                                     >  1:96162/1‑150 (MQ=255)                                                                tGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAgtgg                                                                                    >  2:471312/1‑151 (MQ=255)                                                                 gAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAgtgg                                                                                    <  1:225494/150‑1 (MQ=255)                                                                   ataaaaataaaTCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGgtgt                                                                                 >  2:249144/1‑151 (MQ=255)                                                                      aaaataaaTCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGtttt                                                                              >  2:375161/1‑151 (MQ=255)                                                                      aaaataaaTCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGtttt                                                                              >  1:153552/1‑151 (MQ=255)                                                                           aaatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAAt                                                                          <  2:651835/150‑1 (MQ=255)                                                                                 aTTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGt                                                                   >  2:211315/1‑151 (MQ=255)                                                                                   tAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGt                                                                   >  1:664012/1‑149 (MQ=255)                                                                                   tAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGt                                                                   >  1:44881/1‑149 (MQ=255)                                                                                    aGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAg                                                                 >  1:79833/1‑150 (MQ=255)                                                                                         aTTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtct                                                           <  1:311167/151‑1 (MQ=255)                                                                                            tataTTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtctctc                                                        >  2:197296/1‑151 (MQ=255)                                                                                            tataTTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtctct                                                         <  2:114675/150‑1 (MQ=255)                                                                                             atatTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtctctc                                                        >  1:632666/1‑150 (MQ=255)                                                                                                 tAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAat                                                    <  2:261654/150‑1 (MQ=255)                                                                                                  aaCGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATCTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAatat                                                  <  1:642550/151‑1 (MQ=255)                                                                                                    cGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATc                                                 >  1:55779/1‑150 (MQ=255)                                                                                                    cGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATc                                                 >  1:427559/1‑150 (MQ=255)                                                                                                         gTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCAttt                                             <  1:134955/149‑1 (MQ=255)                                                                                                         gTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAg                                           <  1:211314/151‑1 (MQ=255)                                                                                                         gTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAg                                           <  2:664023/151‑1 (MQ=255)                                                                                                          ttatttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAg                                           >  1:430345/1‑150 (MQ=255)                                                                                                            atttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTgg                                        <  1:331186/151‑1 (MQ=255)                                                                                                            atttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTgg                                        <  1:212802/151‑1 (MQ=255)                                                                                                                ataAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGttt                                     >  2:391122/1‑150 (MQ=255)                                                                                                                  aaGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTtac                                   <  1:647180/150‑1 (MQ=255)                                                                                                                   aGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCaa                                                      <  1:7429/130‑1 (MQ=255)                                                                                                                   aGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCaa                                                      >  2:7428/1‑130 (MQ=255)                                                                                                                   aGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTtacgt                                 <  1:80935/151‑1 (MQ=255)                                                                                                                       ctctATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGAATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTacgtacgt                             >  1:554442/1‑151 (MQ=255)                                                                                                                          tataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTAcgtacgta                            >  1:493502/1‑149 (MQ=255)                                                                                                                          tataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTAcgtacgta                            <  2:493504/149‑1 (MQ=255)                                                                                                                          tataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaa                           >  1:254660/1‑150 (MQ=255)                                                                                                                             aaTTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaaaaa                       <  1:524821/151‑1 (MQ=255)                                                                                                                              aTTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAAt                      >  2:22655/1‑151 (MQ=255)                                                                                                                                                  ccGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTg  >  1:312427/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGGCTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTG  >  NC_000913/576937‑577232

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