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%	191.1 / NA	59	*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 (29/30);  total (29/30)

GTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTAT  >  NC_000913/571312‑571578                                                                                                                                     |                                                                                                                                       gTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgag                                                                                                                                  >  2:225832/1‑139 (MQ=255)    tAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTGAGGTGAGAGt                                                                                                                               <  2:964625/139‑1 (MQ=255)    tAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                               <  1:431711/139‑1 (MQ=255)           tGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTgg                                                                                                                        >  1:8719/1‑139 (MQ=255)            gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGc                                                                                                                       >  2:203968/1‑139 (MQ=255)                 ctAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAagag                                                                                                                  >  1:160743/1‑139 (MQ=255)                  tAAAAACTGAAAAGATCCCTGAAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAgg                                                                                                                 >  2:978583/1‑139 (MQ=255)                     aaaCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTaa                                                                                                              <  1:225832/139‑1 (MQ=255)                       aCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  1:338397/139‑1 (MQ=255)                                                aCTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   <  1:978583/139‑1 (MQ=255)                                                aCTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   <  2:219317/139‑1 (MQ=255)                                                  tCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 <  2:8719/139‑1 (MQ=255)                                                      gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      <  1:662816/98‑1 (MQ=255)                                                      gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      >  2:662816/1‑98 (MQ=255)                                                      gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa                                                                             <  2:370139/139‑1 (MQ=255)                                                            cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAca                                                                       >  1:273806/1‑139 (MQ=255)                                                                      ttATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGAtt                                                             >  2:1108956/1‑139 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            >  2:954254/1‑87 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            >  2:647831/1‑87 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  1:647831/87‑1 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  1:954254/87‑1 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatca                                                                    <  2:1058339/127‑1 (MQ=255)                                                                           tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatca                                                                    >  1:1058339/1‑127 (MQ=255)                                                                            gCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 <  1:933348/113‑1 (MQ=255)                                                                            gCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 >  2:933348/1‑113 (MQ=255)                                                                               tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTga                                                                                                                                   <  1:1026288/60‑1 (MQ=255)                                                                               tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTga                                                                                                                                   >  2:1026288/1‑60 (MQ=255)                                                                                aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTa                                                   <  1:621391/139‑1 (MQ=255)                                                                                  tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                               >  1:226023/1‑61 (MQ=255)                                                                                  tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                               <  2:226023/61‑1 (MQ=255)                                                                                  tGCTAATAAAATCGCAACCGCTGAGATAACACCAGGTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGCTCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAg                                                 <  2:490158/139‑1 (MQ=255)                                                                                     taataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  1:465124/77‑1 (MQ=255)                                                                                     taataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            >  2:465124/1‑77 (MQ=255)                                                                                       ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                            >  1:693247/1‑139 (MQ=255)                                                                                       ataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCctct                                            >  2:725791/1‑139 (MQ=255)                                                                                            aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGTTTAAGCCTCTTATGc                                       >  2:444705/1‑139 (MQ=255)                                                                                                  aCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCtt                                 <  1:1108956/139‑1 (MQ=255)                                                                                                   ccGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCttt                                <  2:273806/139‑1 (MQ=255)                                                                                                      cTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAgg                                                                                                                 <  2:662446/55‑1 (MQ=255)                                                                                                      cTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAgg                                                                                                                 >  1:662446/1‑55 (MQ=255)                                                                                                      cTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGTTTAAGCCTCTTATGCCAGCTTtgtg                             <  1:444705/139‑1 (MQ=255)                                                                                                       tGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATg                                                                                                                            <  2:973876/43‑1 (MQ=255)                                                                                                       tGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATg                                                                                                                            >  1:973876/1‑43 (MQ=255)                                                                                                               cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  1:251684/51‑1 (MQ=255)                                                                                                               cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            >  2:251684/1‑51 (MQ=255)                                                                                                                   aGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCct                                              >  1:200096/1‑109 (MQ=255)                                                                                                                   aGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCct                                              <  2:200096/109‑1 (MQ=255)                                                                                                                     tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatc                                                                     >  2:952067/1‑84 (MQ=255)                                                                                                                     tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatc                                                                     <  1:952067/84‑1 (MQ=255)                                                                                                                     tttATGAGATAAAGTAAGGTGAGAGTAGGCTGACAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAggg              >  2:675656/1‑139 (MQ=255)                                                                                                                      ttATGAGATGAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGga             <  1:725791/139‑1 (MQ=255)                                                                                                                       tATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGgaa            >  1:704916/1‑139 (MQ=255)                                                                                                                          gagaTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   >  2:458234/1‑65 (MQ=255)                                                                                                                          gagaTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   <  1:458234/65‑1 (MQ=255)                                                                                                                             aTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCAt                                                          <  1:1081467/87‑1 (MQ=255)                                                                                                                             aTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCAt                                                          >  2:1081467/1‑87 (MQ=255)                                                                                                                               aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAc                                                                        <  2:798963/71‑1 (MQ=255)                                                                                                                               aaaGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAc                                                                        >  1:798963/1‑71 (MQ=255)                                                                                                                                 aGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTAt  >  2:18323/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                       GTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTAT  >  NC_000913/571312‑571578

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