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%	244.0 / NA	75	*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 (36/39);  total (36/39)

TTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTT  >  NC_000913/571314‑571576                                                                                                                                   |                                                                                                                                     ttAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGagag                                                                                                                              >  1:6606/1‑139 (MQ=255)  tAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                             <  2:396532/139‑1 (MQ=255)          gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgaga                                                                                                                               <  1:649214/129‑1 (MQ=255)          gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgaga                                                                                                                               >  2:649214/1‑129 (MQ=255)          gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGc                                                                                                                     >  1:611946/1‑139 (MQ=255)           ggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    <  2:1641293/139‑1 (MQ=255)                tAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAagg                                                                                                                                    >  2:1501162/1‑118 (MQ=255)                 aaaaaCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                              <  2:6606/139‑1 (MQ=255)                 aaaaaCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                              >  2:318048/1‑139 (MQ=255)                       tGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        <  1:318048/139‑1 (MQ=255)                          aaaGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                          >  1:480866/1‑134 (MQ=255)                          aaaGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                          <  2:480866/134‑1 (MQ=255)                                     tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        >  1:73532/1‑125 (MQ=255)                                     tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        <  2:73532/125‑1 (MQ=255)                                          tCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCTGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCg                                                                                           <  1:1078013/133‑1 (MQ=255)                                          tCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGAAAGACTAACTTAAAGATCg                                                                                           >  2:1078013/1‑133 (MQ=255)                                            tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                                   >  2:1157705/1‑91 (MQ=255)                                            tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                                   <  1:1157705/91‑1 (MQ=255)                                            tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          >  1:205009/1‑132 (MQ=255)                                            tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          <  2:205009/132‑1 (MQ=255)                                                tCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               <  2:54465/139‑1 (MQ=255)                                                 ccAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCaa                                                                                                                   <  2:1432509/102‑1 (MQ=255)                                                 ccAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCaa                                                                                                                   >  1:1432509/1‑102 (MQ=255)                                                    gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATg                                                                                                                          <  2:822469/92‑1 (MQ=255)                                                    gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATg                                                                                                                          >  1:822469/1‑92 (MQ=255)                                                      ggAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATa                                                                         <  1:1606485/139‑1 (MQ=255)                                                              aGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatcat                                                                 <  2:269273/139‑1 (MQ=255)                                                                  ggTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               <  2:1302299/121‑1 (MQ=255)                                                                  ggTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               >  1:1302299/1‑121 (MQ=255)                                                                         tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAg                                                      >  1:961219/1‑139 (MQ=255)                                                                             tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    >  2:520147/1‑73 (MQ=255)                                                                             tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                    <  1:520147/73‑1 (MQ=255)                                                                                tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               >  1:618850/1‑107 (MQ=255)                                                                                tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                               <  2:618850/107‑1 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCt                                                                                                                        >  1:752151/1‑61 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCt                                                                                                                        <  2:752151/61‑1 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          <  2:302200/91‑1 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          >  1:302200/1‑91 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                          >  1:525034/1‑139 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                          >  1:470800/1‑139 (MQ=255)                                                                                     ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                          <  2:961219/139‑1 (MQ=255)                                                                                      taaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatcat                                                                 <  2:454307/115‑1 (MQ=255)                                                                                      taaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatcat                                                                 >  1:454307/1‑115 (MQ=255)                                                                                          aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 >  1:865825/1‑95 (MQ=255)                                                                                          aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  2:865825/95‑1 (MQ=255)                                                                                          aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgag                                                              <  2:662369/114‑1 (MQ=255)                                                                                          aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgag                                                              >  1:662369/1‑114 (MQ=255)                                                                                           tCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                          <  2:1380140/69‑1 (MQ=255)                                                                                           tCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                          >  1:1380140/1‑69 (MQ=255)                                                                                                         aTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 >  2:604879/1‑80 (MQ=255)                                                                                                         aTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  1:604879/80‑1 (MQ=255)                                                                                                             cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        <  1:437501/53‑1 (MQ=255)                                                                                                             cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                        >  2:437501/1‑53 (MQ=255)                                                                                                             cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                          <  2:1484280/115‑1 (MQ=255)                                                                                                             cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                          >  1:1484280/1‑115 (MQ=255)                                                                                                                 aGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCttt                              <  2:258879/123‑1 (MQ=255)                                                                                                                 aGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCttt                              >  1:258879/1‑123 (MQ=255)                                                                                                                  gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               >  1:695845/1‑105 (MQ=255)                                                                                                                  gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               <  2:695845/105‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          >  1:282131/1‑61 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                          <  2:282131/61‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 >  2:937622/1‑70 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 >  2:1664688/1‑70 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  1:1664688/70‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  1:937622/70‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATaaa                                                                       >  2:325315/1‑80 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATaaa                                                                       <  1:325315/80‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               >  2:854222/1‑104 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               <  1:854222/104‑1 (MQ=255)                                                                                                                        gagaTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTTAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               <  2:993964/99‑1 (MQ=255)                                                                                                                        gagaTAAAGTAAGGTGAGAGTATGCTGCCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                               >  1:993964/1‑99 (MQ=255)                                                                                                                             aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAACGATCGAAGAAATACGGATGCATATACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGtt  <  1:229059/139‑1 (MQ=255)                                                                                                                             aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 >  1:139901/1‑60 (MQ=255)                                                                                                                             aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                 <  2:139901/60‑1 (MQ=255)                                                                                                                               aGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCa                                                         <  2:525217/82‑1 (MQ=255)                                                                                                                               aGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCa                                                         >  1:525217/1‑82 (MQ=255)                                                                                                                                   |                                                                                                                                     TTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTT  >  NC_000913/571314‑571576

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