Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,204,160	G→A	*113* (TAG→TAA)	ymfM →	e14 prophage, uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,204,160	0	G	A	100.0%	190.4 / NA	59	*113* (TAG→TAA)	ymfM	e14 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/31);  total (28/31)

ACTACGGAGTGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAGAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGATTTCAGTAATGACGAAGGTGTGTGCTGGCCA  >  NC_000913/1204032‑1204292                                                                                                                                 |                                                                                                                                     aCTACGGAGTGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAggt                                                                                                                                  >  2:1340108/1‑133 (MQ=255) aCTACGGAGTGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAggt                                                                                                                                  <  1:1340108/133‑1 (MQ=255)          tGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTa                                                                                                                   >  2:692774/1‑139 (MQ=255)          tGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTa                                                                                                                   >  1:1187843/1‑139 (MQ=255)          tGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTa                                                                                                                   >  2:44302/1‑139 (MQ=255)                 cAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAAc                                                                                                            <  2:1015937/139‑1 (MQ=255)                   gagTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCg                                                                                                          <  2:1187843/139‑1 (MQ=255)                    agTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCgg                                                                                                         >  1:358993/1‑139 (MQ=255)                    agTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCgg                                                                                                         >  2:1146107/1‑139 (MQ=255)                    agTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCgg                                                                                                         <  2:1326742/139‑1 (MQ=255)                      tCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCt                                                                                                       <  1:44302/139‑1 (MQ=255)                      tCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCt                                                                                                       <  1:737916/139‑1 (MQ=255)                                        aaaCACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcg                                                                                     <  2:358993/139‑1 (MQ=255)                                              taCGGTTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTATCCGTCTATGTATGGGATGGTTGCGCTGCAt                                                                               <  2:84828/138‑1 (MQ=255)                                                            cTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCAt                                                                                                                           <  2:370612/81‑1 (MQ=255)                                                            cTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCAt                                                                                                                           >  1:370612/1‑81 (MQ=255)                                                              ccGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCa                                                                         >  2:756316/1‑129 (MQ=255)                                                              ccGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCa                                                                         <  1:756316/129‑1 (MQ=255)                                                                gcgAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCa                                                             >  2:1329691/1‑139 (MQ=255)                                                                 cgAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAg                                                            <  1:1329691/139‑1 (MQ=255)                                                                     ggTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCAt                                                                                                                           >  1:331855/1‑72 (MQ=255)                                                                     ggTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCAt                                                                                                                           <  2:331855/72‑1 (MQ=255)                                                                                tGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGc                                                                                                                       <  1:583934/65‑1 (MQ=255)                                                                                tGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGc                                                                                                                       >  2:583934/1‑65 (MQ=255)                                                                                tGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcgc                                                                                    <  2:236339/100‑1 (MQ=255)                                                                                tGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcgc                                                                                    >  1:236339/1‑100 (MQ=255)                                                                                       ttCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCta                                                                                                      <  1:471804/75‑1 (MQ=255)                                                                                       ttCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCta                                                                                                      >  2:471804/1‑75 (MQ=255)                                                                                        tCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCt                                                                                                       >  2:354620/1‑73 (MQ=255)                                                                                        tCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCt                                                                                                       <  1:354620/73‑1 (MQ=255)                                                                                         cAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGc                                                                                                                       <  1:1169/56‑1 (MQ=255)                                                                                         cAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGc                                                                                                                       >  2:1169/1‑56 (MQ=255)                                                                                         cAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAAcc                                                                                                           >  1:1082683/1‑68 (MQ=255)                                                                                         cAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAAcc                                                                                                           <  2:1082683/68‑1 (MQ=255)                                                                                            cccGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGc                                                                                                                       <  1:1037751/53‑1 (MQ=255)                                                                                            cccGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGc                                                                                                                       >  2:1037751/1‑53 (MQ=255)                                                                                             ccGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGAt                                <  1:257998/139‑1 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTa                                                                                                              <  1:1134890/58‑1 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTa                                                                                                              >  2:1134890/1‑58 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCgg                                                                                                         <  2:1167230/63‑1 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCgg                                                                                                         >  1:1167230/1‑63 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCAt                                                                               >  1:388467/1‑89 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCAt                                                                               <  2:388467/89‑1 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATg                                                                       >  2:1111657/1‑97 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATg                                                                       <  1:1111657/97‑1 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGaa                                                                     >  1:1113463/1‑99 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGaa                                                                     <  2:1113463/99‑1 (MQ=255)                                                                                                cTCGAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGATTTc                             >  2:1189351/1‑139 (MQ=255)                                                                                                   gAACAGTTTCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGATTTCAGt                          <  1:692774/139‑1 (MQ=255)                                                                                                         tttCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcg                                                                                     <  2:279691/74‑1 (MQ=255)                                                                                                         tttCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcg                                                                                     >  1:279691/1‑74 (MQ=255)                                                                                                         tttCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGATTTCAGTAATGAc                    >  1:713137/1‑139 (MQ=255)                                                                                                          ttCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcgc                                                                                    <  1:84532/74‑1 (MQ=255)                                                                                                          ttCGTCGTAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTgcgc                                                                                    >  2:84532/1‑74 (MQ=255)                                                                                                               cgtAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCAt                                                                               >  1:757250/1‑74 (MQ=255)                                                                                                               cgtAAATTCAGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCAt                                                                               <  2:757250/74‑1 (MQ=255)                                                                                                                        aGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGc                                                                                                        >  2:1273402/1‑40 (MQ=255)                                                                                                                        aGGGAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGc                                                                                                        <  1:1273402/40‑1 (MQ=255)                                                                                                                           gAAATAAAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGATTTCAGTAATGACGAAGGTGTGTGCTGGCCa  <  2:713137/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     ACTACGGAGTGCCCGTCAGAGTCATTCGTTGGGAGCCGGAAACACAACGGGTTATCTACCTCCGCGAAGGTTATGAGCATGAATGCTTCAGCCCGCTCGAACAGTTTCGTCGTAAATTCAGGGAAATAGAGGTCGGTCATGAGCACTAAATTAACCGGCTATGTATGGGATGGTTGCGCTGCATCAGGCATGAAGTTATCCAGCGTGGCAATTATGGCCCGCCTGGCTGATTTCAGTAATGACGAAGGTGTGTGCTGGCCA  >  NC_000913/1204032‑1204292

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