Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,296,984	G→A	*192* (TAG→TAA)	yraP →	outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,296,984	0	G	A	100.0%	129.3 / NA	42	*192* (TAG→TAA)	yraP	outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/23);  total (19/23)

AAAGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGAACAA  >  NC_000913/3296850‑3297106                                                                                                                                       |                                                                                                                           aaaGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaacagc                                                                                                                        <  2:161364/139‑1 (MQ=255)  aaGTGACCACCGAATACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcc                                                                                                                       <  1:833465/139‑1 (MQ=255)  aaGTGACCACCGAAAACGGTGCAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcc                                                                                                                       <  1:759817/139‑1 (MQ=255)    gTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcccc                                                                                                                     <  2:390439/139‑1 (MQ=255)       accaccGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTg                                                                                                                  >  1:1019276/1‑139 (MQ=255)                     tGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCg                                                                                                    >  2:591829/1‑139 (MQ=255)                     tGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCg                                                                                                    <  1:1119856/139‑1 (MQ=255)                             tCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            <  1:591829/139‑1 (MQ=255)                                    ggggCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                     >  2:457747/1‑139 (MQ=255)                                      ggCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                          >  2:631231/1‑100 (MQ=255)                                       gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                  <  2:499884/139‑1 (MQ=255)                                        cTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                 >  2:1066523/1‑139 (MQ=255)                                             gACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                            <  1:457747/139‑1 (MQ=255)                                               cTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                          <  2:57986/139‑1 (MQ=255)                                                  aaCGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaa                                                                       <  1:87224/139‑1 (MQ=255)                                                               aGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGt                                                          <  1:1066523/139‑1 (MQ=255)                                                                 cggcggCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                               <  2:1009441/116‑1 (MQ=255)                                                                 cggcggCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                               >  1:1009441/1‑116 (MQ=255)                                                                        aGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGAt                                                 <  2:202915/139‑1 (MQ=255)                                                                           tatTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAg                                              <  1:462582/139‑1 (MQ=255)                                                                                   gCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAc                                      >  2:734074/1‑139 (MQ=255)                                                                                    ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            >  1:528833/1‑84 (MQ=255)                                                                                    ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            <  2:528833/84‑1 (MQ=255)                                                                                    ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGt                                                                                       <  1:354180/89‑1 (MQ=255)                                                                                    ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGt                                                                                       >  2:354180/1‑89 (MQ=255)                                                                                           aGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                          <  2:1130453/47‑1 (MQ=255)                                                                                           aGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                          >  1:1130453/1‑47 (MQ=255)                                                                                                  tggaGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAAACGTAGGGCAGATAAGGCGTTTACGCCGCATCGGACAtc                       >  2:970737/4‑138 (MQ=37)                                                                                                  tGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGGCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGtt                                         <  1:807676/121‑1 (MQ=255)                                                                                                  tGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGtt                                         >  2:807676/1‑121 (MQ=255)                                                                                                              cTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            >  2:124004/1‑58 (MQ=255)                                                                                                              cTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            <  1:124004/58‑1 (MQ=255)                                                                                                                 cGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAAt          <  2:881111/137‑1 (MQ=255)                                                                                                                 cGGAGTTTACGTTTATTAAAAAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAAt          >  1:881111/1‑137 (MQ=255)                                                                                                                   gCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                            <  1:205224/69‑1 (MQ=255)                                                                                                                   gCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                            >  2:205224/1‑69 (MQ=255)                                                                                                                    cgtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            <  1:408845/52‑1 (MQ=255)                                                                                                                    cgtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                            >  2:408845/1‑52 (MQ=255)                                                                                                                    cgtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                  >  2:473283/1‑62 (MQ=255)                                                                                                                    cgtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                  <  1:473283/62‑1 (MQ=255)                                                                                                                    cgtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGaa     >  2:744568/1‑139 (MQ=255)                                                                                                                     gtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACa                         >  1:648998/1‑118 (MQ=255)                                                                                                                     gtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACa                         <  2:648998/118‑1 (MQ=255)                                                                                                                       ttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGaacaa  >  2:743445/1‑139 (MQ=255)                                                                                                                        tacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                  <  1:143028/58‑1 (MQ=255)                                                                                                                        tacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                  >  2:143028/1‑58 (MQ=255)                                                                                                                                       |                                                                                                                           AAAGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGAACAA  >  NC_000913/3296850‑3297106

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