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%	179.0 / NA	55	*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 (26/29);  total (26/29)

GTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGA  >  NC_000913/571312‑571571                                                                                                                                     |                                                                                                                                gTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGAAAGGTgag                                                                                                                           >  1:389799/1‑139 (MQ=255)    tAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                        >  2:267690/1‑139 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATATCACCAGTTTATGAGATAAAGTAAGGt                                                                                                                              >  2:949538/1‑125 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                              >  1:1145420/1‑125 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                              <  2:1145420/125‑1 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                              <  1:949538/125‑1 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGc                                                                                                                >  2:1248529/1‑139 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGc                                                                                                                >  2:1321423/1‑139 (MQ=255)             ggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                               <  1:267690/139‑1 (MQ=255)             ggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                               <  1:135201/139‑1 (MQ=255)                   aaaaaCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                         >  2:1209430/1‑139 (MQ=255)                        cTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGAc                                                                                                    <  1:1248529/139‑1 (MQ=255)                                aTTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTaa                                                                                                       <  2:586254/128‑1 (MQ=255)                                aTTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTaa                                                                                                       >  1:586254/1‑128 (MQ=255)                                  tCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAg                                                                                          >  2:383646/1‑139 (MQ=255)                                       tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                     <  2:1053222/139‑1 (MQ=255)                                            tCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATGACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGaaa                                                                                >  1:660037/1‑139 (MQ=255)                                              tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATGACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                                                                              <  2:660037/139‑1 (MQ=255)                                              tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                              >  2:817012/1‑91 (MQ=255)                                              tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGt                                                                                                                              <  1:817012/91‑1 (MQ=255)                                                 cTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                           <  1:40275/139‑1 (MQ=255)                                                      gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa                                                                      >  2:936463/1‑139 (MQ=255)                                                      gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa                                                                      <  1:383646/139‑1 (MQ=255)                                                        ggAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATa                                                                    <  2:389799/139‑1 (MQ=255)                                                        ggAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATa                                                                    <  1:1321423/139‑1 (MQ=255)                                                               tAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                     >  1:337653/1‑99 (MQ=255)                                                               tAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                     <  2:337653/99‑1 (MQ=255)                                                                       tataTGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTc                                                     >  1:1436187/1‑139 (MQ=255)                                                                        atatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                               >  2:529501/1‑80 (MQ=255)                                                                        atatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                               <  1:529501/80‑1 (MQ=255)                                                                         tatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCAt                                                   >  1:27916/1‑139 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                           >  1:706377/1‑113 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                           <  2:706377/113‑1 (MQ=255)                                                                              ttAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGt                                              <  2:27916/139‑1 (MQ=255)                                                                               tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgag                                                         >  1:194087/1‑127 (MQ=255)                                                                               tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgag                                                         <  2:194087/127‑1 (MQ=255)                                                                                       ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                     >  2:582763/1‑139 (MQ=255)                                                                                            aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                         >  2:595440/1‑66 (MQ=255)                                                                                            aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                         <  1:595440/66‑1 (MQ=255)                                                                                              cGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                     <  1:879899/68‑1 (MQ=255)                                                                                              cGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                     >  2:879899/1‑68 (MQ=255)                                                                                                 aaCCGGTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACGTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCt                           <  1:936463/139‑1 (MQ=255)                                                                                                 aaCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCt                           <  1:582763/139‑1 (MQ=255)                                                                                                             aaCACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCt                                                                                                                   >  1:287998/1‑39 (MQ=255)                                                                                                             aaCACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCt                                                                                                                   <  2:287998/39‑1 (MQ=255)                                                                                                                  cAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCa          <  1:298715/139‑1 (MQ=255)                                                                                                                    gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                           <  1:690266/72‑1 (MQ=255)                                                                                                                    gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACgg                                                                           >  2:690266/1‑72 (MQ=255)                                                                                                                    gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                            >  1:191176/1‑71 (MQ=255)                                                                                                                    gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                            <  2:191176/71‑1 (MQ=255)                                                                                                                     tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCt                           <  2:586709/119‑1 (MQ=255)                                                                                                                     tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCt                           >  1:586709/1‑119 (MQ=255)                                                                                                                       tATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGTCAGCTTTGTGTCGTATCCGTCAGGgaa     <  1:1209430/139‑1 (MQ=255)                                                                                                                         tGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAAGTAAGACTAACTTAAAGATCGAAGAAATACTGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGgaaag   <  2:1436187/139‑1 (MQ=255)                                                                                                                          gagaTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGgaaaga  >  1:1027222/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                GTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGA  >  NC_000913/571312‑571571

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