Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,203,224	G→A	*67* (TAG→TAA)	croE →	e14 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,203,224	0	G	A	100.0%	185.8 / NA	57	*67* (TAG→TAA)	croE	e14 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/34);  total (23/34)

ACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAGCAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATCGCGGCGTTGCCGAG  >  NC_000913/1203100‑1203350                                                                                                                             |                                                                                                                               aCTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCaa                                                                                                                  <  2:288096/139‑1 (MQ=255)      cGGCTGGTTTTCAATGGCTATAAAATAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCAGATATCTATCTGAAAGATTAACAGAACACTTTCAAttttt                                                                                                             <  2:1433788/139‑1 (MQ=255)      cGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAAttttt                                                                                                             <  1:329986/139‑1 (MQ=255)      cGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAAttttt                                                                                                             <  2:1873288/139‑1 (MQ=255)      cGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAAttttt                                                                                                             <  1:1498862/139‑1 (MQ=255)         cTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAAc                                                                                                          >  1:1793332/1‑139 (MQ=255)          tGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAAcc                                                                                                         <  2:590772/139‑1 (MQ=255)              tttCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAg                                                                                                     <  2:1089028/139‑1 (MQ=255)                    tGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               <  2:1552681/139‑1 (MQ=255)                    tGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               <  2:506107/139‑1 (MQ=255)                       cTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAg                                                                                            >  1:1608901/1‑139 (MQ=255)                         ataAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGc                                                                                          <  2:1839656/139‑1 (MQ=255)                          taAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCt                                                                                         >  1:434132/1‑139 (MQ=255)                             aaaaGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAAt                                                                                      <  2:49192/139‑1 (MQ=255)                              aaaGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATc                                                                                     <  1:1877844/139‑1 (MQ=255)                              aaaGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATc                                                                                     <  2:1613007/139‑1 (MQ=255)                                agacAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACAATTTCAATTTTTAACAACAGAACGATGAgg                                                                                           >  1:327009/4‑131 (MQ=255)                                aGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAgg                                                                                           <  2:327009/131‑1 (MQ=255)                                 gCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTg                                                                                  <  2:1793332/139‑1 (MQ=255)                                     gTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTCTCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTa                                                                              >  1:974543/1‑139 (MQ=255)                                       tttGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAg                                                                            <  1:1768841/139‑1 (MQ=255)                                                          ttGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               <  2:462361/101‑1 (MQ=255)                                                          ttGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               >  1:462361/1‑101 (MQ=255)                                                             aGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               <  1:6069/98‑1 (MQ=255)                                                             aGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               >  2:6069/1‑98 (MQ=255)                                                                aaTACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGc                                                   <  2:434132/139‑1 (MQ=255)                                                                       aTCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAAc                                                      >  1:109596/1‑129 (MQ=255)                                                                       aTCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAAc                                                      <  2:109596/129‑1 (MQ=255)                                                                       aTCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTg                                            <  2:1429504/139‑1 (MQ=255)                                                                        tCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTgg                                           <  1:500848/139‑1 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAGGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGc                                                   >  1:1514779/1‑127 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGaa                                                                 <  1:878083/113‑1 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGaa                                                                 >  2:878083/1‑113 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGc                                                   <  2:1514779/127‑1 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCt                                                 >  2:1035183/1‑129 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCt                                                 <  1:1035183/129‑1 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTgg                                           >  2:1587220/1‑135 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTgg                                           <  1:1587220/135‑1 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACg                                       >  2:1191561/1‑139 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACg                                       >  2:1630870/1‑139 (MQ=255)                                                                            gTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACg                                       >  1:555629/1‑139 (MQ=255)                                                                                    aCGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCt                               <  2:555629/139‑1 (MQ=255)                                                                                              aCTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACca                                                                                                        >  2:1589802/1‑56 (MQ=255)                                                                                              aCTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACca                                                                                                        <  1:1589802/56‑1 (MQ=255)                                                                                              aCTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAg                                                                            <  1:421043/84‑1 (MQ=255)                                                                                              aCTTAAGACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAg                                                                            >  2:421043/1‑84 (MQ=255)                                                                                                    gACCGGATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATcg               <  1:1191561/139‑1 (MQ=255)                                                                                                        ggATATCTATCCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATCgcggc           >  2:907817/1‑139 (MQ=255)                                                                                                               tatcCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATCGCGGCGTTGCCg    <  1:1630870/139‑1 (MQ=255)                                                                                                                atcCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATCGCGGCGTTGCCGa   >  2:697323/1‑139 (MQ=255)                                                                                                                 tcCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAAc                                                      <  2:1669031/87‑1 (MQ=255)                                                                                                                 tcCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAAc                                                      >  1:1669031/1‑87 (MQ=255)                                                                                                                 tcCGAAAGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATCGCGGCGTTGCCGAg  >  1:1150768/1‑139 (MQ=255)                                                                                                                      aaGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTaa                                                                                       <  2:1071288/49‑1 (MQ=255)                                                                                                                      aaGATTAACAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTaa                                                                                       >  1:1071288/1‑49 (MQ=255)                                                                                                                         aTTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               >  1:842126/1‑38 (MQ=255)                                                                                                                         aTTAACAGAACACTTTCAATTTTTAACCACAGAACGAt                                                                                               <  2:842126/38‑1 (MQ=255)                                                                                                                             |                                                                                                                               ACTTACGGCTGGTTTTCAATGGCTATAAAAAAGCCAGTTTTGTGCTGGCTAAAAAACTTGAGCAATACACATCAGGTGCAATTACGAAATCTGACTTAAGACCGGATATCTATCCGAAAGATTAGCAGAACACTTTCAATTTTTAACCACAGAACGATGAGGCTAATCGTGGGTAAGCATCACTGGAAAATAGAAAAACAGCCTGAGTGGTACGTGAAAGCTGTCAGAAAAACTATCGCGGCGTTGCCGAG  >  NC_000913/1203100‑1203350

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: ‑