Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,113,616	C→A	V108V (GTC→GTA)	uspD →	universal stress protein D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,616	0	C	A	100.0%	189.0 / NA	54	V108V (GTC→GTA)	uspD	universal stress protein D
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/31);  total (23/31)

ATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCT  >  NC_000913/4113483‑4113744                                                                                                                                      |                                                                                                                                 aTAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGt                                                                                                                            >  2:987936/1‑140 (MQ=255)    aGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGGCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGtcat                                                                                                                         >  2:85573/1‑140 (MQ=255)    aGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGtcat                                                                                                                         >  1:627756/1‑140 (MQ=255)        ggATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCATGTATGTGGtcatcatc                                                                                                                     <  1:322935/140‑1 (MQ=255)            aaCAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTc                                                                                                                 >  1:349551/1‑140 (MQ=255)                aGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGTGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCAttt                                                                                                             >  2:30321/1‑140 (MQ=255)                          aCTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGtt                                                                                                   <  1:27748/140‑1 (MQ=255)                             gACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGa                                                                                                >  1:847112/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  1:167720/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  2:743909/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  1:822878/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  1:417489/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  1:152839/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  1:303314/1‑140 (MQ=255)                               cGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATg                                                                                              >  1:264957/1‑140 (MQ=255)                                 aaaaaTATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGcc                                                                                            >  1:934439/1‑140 (MQ=255)                                    aaTATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGc                                                                                         <  2:935237/140‑1 (MQ=255)                                    aaTATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGc                                                                                         <  2:735991/140‑1 (MQ=255)                                          cAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATcg                                                                                   <  2:864204/140‑1 (MQ=255)                                           aaTGGGCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATcgc                                                                                  <  1:85573/140‑1 (MQ=255)                                             tGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCgg                                                                                <  2:779868/140‑1 (MQ=255)                                               gCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGga                                                                              <  2:303314/140‑1 (MQ=255)                                                ccGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGgat                                                                             >  1:395797/1‑140 (MQ=255)                                                ccGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTAT‑‑‑GTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCggg                                                                               >  2:791137/1‑135 (MQ=255)                                                ccGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTAT‑‑‑GTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCggg                                                                               <  1:791137/135‑1 (MQ=255)                                                 cGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGgatg                                                                            <  2:512158/140‑1 (MQ=255)                                                  gAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGgatga                                                                           <  2:375965/140‑1 (MQ=255)                                                  gAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGgatga                                                                           <  2:417489/140‑1 (MQ=255)                                                          aaCTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCaa                                                                       <  1:575648/136‑1 (MQ=255)                                                          aaCTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCaa                                                                       >  2:575648/1‑136 (MQ=255)                                                          aaCTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAg                                                                   <  2:264957/140‑1 (MQ=255)                                                          aaCTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAg                                                                   <  1:700900/140‑1 (MQ=255)                                                          aaCTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAg                                                                   >  1:226595/1‑140 (MQ=255)                                                              gCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGa                                                               >  2:18000/1‑140 (MQ=255)                                                               cGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAg                                                              >  1:284741/1‑140 (MQ=255)                                                                gTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGt                                                             >  2:95718/1‑140 (MQ=255)                                                                 tATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTg                                                            <  1:743909/140‑1 (MQ=255)                                                                         gcgGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTg                                                    >  2:187015/1‑140 (MQ=255)                                                                             agaAATGCCGGAAACACTGCTGGAAATTATGCAAAACGAACAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCa                                                <  1:996014/140‑1 (MQ=255)                                                                             agaAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCa                                                <  2:391173/140‑1 (MQ=255)                                                                             agaAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCa                                                <  2:442640/140‑1 (MQ=255)                                                                                 aTGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGt                                            <  2:934439/140‑1 (MQ=255)                                                                                  tGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTg                                           <  2:152839/140‑1 (MQ=255)                                                                                   gCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGc                                          <  2:533856/140‑1 (MQ=255)                                                                                    ccGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGcc                                         <  2:686744/140‑1 (MQ=255)                                                                                       gAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGtt                                      <  2:627756/140‑1 (MQ=255)                                                                                             ctgctgGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGGGTGCGGATTTGCTCATCGTGCCGTTTATCGa                                <  2:167720/140‑1 (MQ=255)                                                                                                       ttATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCg                      <  1:187015/140‑1 (MQ=255)                                                                                                       ttATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCg                      <  1:95718/140‑1 (MQ=255)                                                                                                       ttATGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCg                      <  2:395797/140‑1 (MQ=255)                                                                                                          tGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCg                      >  2:529498/1‑137 (MQ=255)                                                                                                          tGCAAAAAGAGCAGTGCGACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCg                      <  1:529498/137‑1 (MQ=255)                                                                                                                           gACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCt  <  1:18000/140‑1 (MQ=255)                                                                                                                           gACCTCCTTGTATGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCt  <  2:284741/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 ATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCT  >  NC_000913/4113483‑4113744

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑