Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	571,444	G→A	*184* (TAG→TAA)	ybcL →	inactive polymorphonuclear leukocyte migration suppressor, DLP12 prophage, UPF0098 family secreted protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	571,444	0	G	A	100.0%	203.3 / NA	63	*184* (TAG→TAA)	ybcL	inactive polymorphonuclear leukocyte migration suppressor, DLP12 prophage, UPF0098 family secreted protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/32);  total (31/32)

CAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTT  >  NC_000913/571310‑571576                                                                                                                                       |                                                                                                                                     cAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTg                                                                                                                                  <  2:777251/139‑1 (MQ=255)        aaGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTAt                                                                                                                           >  1:458341/1‑139 (MQ=255)              gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgag                                                                                                                                <  2:963485/128‑1 (MQ=255)              gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgag                                                                                                                                >  1:963485/1‑128 (MQ=255)                  tctAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAaga                                                                                                                 >  1:840320/1‑139 (MQ=255)                         aCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                          <  2:991818/139‑1 (MQ=255)                           tGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        <  1:1258194/139‑1 (MQ=255)                                         tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          <  2:693181/139‑1 (MQ=255)                                         tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATACTGGCAAGAGGt                                                                                                              <  2:782743/119‑1 (MQ=255)                                         tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATACTGGCAAGAGGt                                                                                                              >  1:782743/1‑119 (MQ=255)                                            aTTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCaa                                                                                                                   <  1:37572/111‑1 (MQ=255)                                            aTTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCaa                                                                                                                   >  2:37572/1‑111 (MQ=255)                                                tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                                                                                   <  2:1170999/139‑1 (MQ=255)                                                  aCTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  1:639364/139‑1 (MQ=255)                                                        gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               >  1:405297/1‑135 (MQ=255)                                                        gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               <  2:405297/135‑1 (MQ=255)                                                              cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        >  2:1094618/1‑104 (MQ=255)                                                              cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        <  1:1094618/104‑1 (MQ=255)                                                                      ggTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgaga                                                             >  2:708442/1‑139 (MQ=255)                                                                      ggTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgaga                                                             >  2:1002886/1‑139 (MQ=255)                                                                          atatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCa                                                         <  2:458341/139‑1 (MQ=255)                                                                           tatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCAt                                                        >  1:85142/1‑139 (MQ=255)                                                                           tatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCACAAACATCATGAGATTCa                                                         >  2:588174/1‑138 (MQ=255)                                                                           tatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCACAAACATCATGAGATTCa                                                         <  1:588174/138‑1 (MQ=255)                                                                                 tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                             >  1:373903/1‑64 (MQ=255)                                                                                 tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                             <  2:373903/64‑1 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgag                                                                                                                                <  1:404751/60‑1 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgag                                                                                                                                >  2:404751/1‑60 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               <  2:843357/109‑1 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               >  1:381333/1‑109 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               >  1:843357/1‑109 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               <  2:381333/109‑1 (MQ=255)                                                                                  aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTa                                                 >  1:1021472/1‑139 (MQ=255)                                                                                    tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               <  1:1002886/139‑1 (MQ=255)                                                                                    tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               <  1:708442/139‑1 (MQ=255)                                                                                         ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                                                                                   <  1:463764/98‑1 (MQ=255)                                                                                         ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                                                                                   >  2:463764/1‑98 (MQ=255)                                                                                         ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               >  2:1056941/1‑134 (MQ=255)                                                                                         ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               <  1:1056941/134‑1 (MQ=255)                                                                                          taaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCtt                                         >  2:465732/1‑139 (MQ=255)                                                                                              aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    <  1:601983/60‑1 (MQ=255)                                                                                              aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    >  2:601983/1‑60 (MQ=255)                                                                                              aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 >  1:581069/1‑95 (MQ=255)                                                                                              aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  2:581069/95‑1 (MQ=255)                                                                                              aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGc                                     <  2:1021472/139‑1 (MQ=255)                                                                                              aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGCTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGACGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGc                                     <  2:875984/139‑1 (MQ=255)                                                                                                        cTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTtgtg                           >  1:106614/1‑139 (MQ=255)                                                                                                         tGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTgtgt                          >  1:633631/1‑139 (MQ=255)                                                                                                            gaTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGt                       <  1:465732/139‑1 (MQ=255)                                                                                                             aTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa                                                                           >  2:651013/1‑86 (MQ=255)                                                                                                             aTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa                                                                           <  1:651013/86‑1 (MQ=255)                                                                                                                  acCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGt                 <  2:106614/139‑1 (MQ=255)                                                                                                                      gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                                <  1:1236934/72‑1 (MQ=255)                                                                                                                      gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                                >  2:1236934/1‑72 (MQ=255)                                                                                                                       tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    >  2:47931/1‑35 (MQ=255)                                                                                                                       tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    <  1:47931/35‑1 (MQ=255)                                                                                                                       tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatc                                                                   <  2:58233/84‑1 (MQ=255)                                                                                                                       tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatc                                                                   >  1:58233/1‑84 (MQ=255)                                                                                                                         tATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGgaa          >  1:889681/1‑139 (MQ=255)                                                                                                                              gaTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGgaaagaaa     >  2:1129451/1‑139 (MQ=255)                                                                                                                                 aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCAt                                                        >  2:687667/1‑85 (MQ=255)                                                                                                                                 aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCAt                                                        <  1:687667/85‑1 (MQ=255)                                                                                                                                 aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGtt  >  2:699938/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                     CAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTT  >  NC_000913/571310‑571576

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