Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,412,000	C→T	*412* (TAG→TAA)	intR ←	Rac prophage, integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,412,000	0	C	T	100.0%	172.6 / NA	53	*412* (TAG→TAA)	intR	Rac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/28);  total (25/28)

CGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGC  >  NC_000913/1411867‑1412134                                                                                                                                      |                                                                                                                                       cGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGa                                                                                                                                   >  2:156852/1‑139 (MQ=255)       cGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                             <  2:518842/139‑1 (MQ=255)                 cAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                                   <  2:1398823/139‑1 (MQ=255)                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCa                                                                                                                     >  1:1072403/1‑136 (MQ=255)                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCa                                                                                                                     <  2:1072403/136‑1 (MQ=255)                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                                   <  2:557182/138‑1 (MQ=255)                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                                   >  1:557182/1‑138 (MQ=255)                              aaTTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                      <  1:156852/139‑1 (MQ=255)                              aaTTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGAATTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                      <  1:665977/139‑1 (MQ=255)                                tttGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                         <  2:1301083/118‑1 (MQ=255)                                tttGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                         >  1:1301083/1‑118 (MQ=255)                                          ttGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGt                                                                                          <  2:307704/139‑1 (MQ=255)                                                   tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGt                                                                                          <  2:305272/130‑1 (MQ=255)                                                   tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGt                                                                                          >  1:305272/1‑130 (MQ=255)                                                   tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATcc                                                                                 >  1:647007/1‑139 (MQ=255)                                                      ttctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcaca                                                                              <  1:246986/139‑1 (MQ=255)                                                       tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                             >  2:738457/1‑139 (MQ=255)                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                                       >  1:138531/1‑92 (MQ=255)                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                                       <  2:138531/92‑1 (MQ=255)                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATg                                                                                                                <  2:790571/99‑1 (MQ=255)                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATg                                                                                                                >  1:790571/1‑99 (MQ=255)                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                                        <  1:165601/139‑1 (MQ=255)                                                             tAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTa                                                                       <  2:1407423/139‑1 (MQ=255)                                                                       aTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                             <  1:350745/75‑1 (MQ=255)                                                                       aTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                             >  2:350745/1‑75 (MQ=255)                                                                              tGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcac                                                                               <  2:890076/114‑1 (MQ=255)                                                                              tGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcac                                                                               >  1:890076/1‑114 (MQ=255)                                                                              tGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGAAATTACGATAACCGCGCATTCAACACGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCt                                                      >  2:316393/1‑139 (MQ=255)                                                                                   tattTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCa                                                 <  1:285052/139‑1 (MQ=255)                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                         >  1:865909/1‑65 (MQ=255)                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                         <  2:865909/65‑1 (MQ=255)                                                                                         aGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATAcc                                                                          >  2:1481765/1‑108 (MQ=255)                                                                                         aGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATAcc                                                                          <  1:1481765/108‑1 (MQ=255)                                                                                           ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                            <  2:235142/88‑1 (MQ=255)                                                                                           ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                            >  1:235142/1‑88 (MQ=255)                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGcc                                       >  1:661031/1‑139 (MQ=255)                                                                                               ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGcccc                                     >  1:197247/1‑139 (MQ=255)                                                                                                     cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGt                                                                                          >  1:554832/1‑80 (MQ=255)                                                                                                     cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGt                                                                                          <  2:554832/80‑1 (MQ=255)                                                                                                     cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                               <  2:1111473/139‑1 (MQ=255)                                                                                                        aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctg                            >  2:1364381/1‑139 (MQ=255)                                                                                                        aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctg                            <  1:1429327/139‑1 (MQ=255)                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                         <  1:1001802/41‑1 (MQ=255)                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                         >  2:1001802/1‑41 (MQ=255)                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                             >  2:396254/1‑85 (MQ=255)                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                             <  1:396254/85‑1 (MQ=255)                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                       >  2:344160/1‑139 (MQ=255)                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                       >  1:744773/1‑139 (MQ=255)                                                                                                              gCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtata                      <  1:344160/139‑1 (MQ=255)                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTgcg                                                   >  2:781408/1‑106 (MQ=255)                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTgcg                                                   <  1:781408/106‑1 (MQ=255)                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATg                  >  2:623399/1‑139 (MQ=255)                                                                                                                                  gCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGc  <  2:744773/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       CGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGC  >  NC_000913/1411867‑1412134

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