Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,090,672	C→T	G85D (GGC→GAC)	frvX ←	peptidase M42 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,090,672	0	C	T	100.0%	171.8 / NA	49	G85D (GGC→GAC)	frvX	peptidase M42 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/26);  total (23/26)

ATCAAATGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGCCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAGCTGCCCA  >  NC_000913/4090536‑4090801                                                                                                                                           |                                                                                                                                  atcaAATGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCTTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAgtcgc                                                                                                                                >  1:653261/1‑140 (MQ=255)   caAATGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAgtcgc                                                                                                                                <  1:946012/138‑1 (MQ=255)   caAATGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAgtcgc                                                                                                                                >  2:946012/1‑138 (MQ=255)    aaaTGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCaa                                                                                                                             <  1:803625/140‑1 (MQ=255)       tGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAAtgg                                                                                                                          <  2:589729/140‑1 (MQ=255)            gCGGTTGTTGCTTTTGCTTTTCCGCTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTa                                                                                                                     <  1:345413/140‑1 (MQ=255)                   ttgcttttgcttTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCa                                                                                                              <  2:559492/140‑1 (MQ=255)                          tgcttttCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAAcc                                                                                                       >  2:1215544/1‑140 (MQ=255)                                  cGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCtcgt                                                                                               <  1:1215544/140‑1 (MQ=255)                                     tAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGa                                                                                            >  1:789327/1‑140 (MQ=255)                                        cGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGTGCGCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGAt                                                                                           <  1:18636/140‑1 (MQ=255)                                          cATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCACTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTgg                                                                                       <  2:1124662/140‑1 (MQ=255)                                                      cGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAg                                                                           <  1:118544/140‑1 (MQ=255)                                                       ggAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGc                                                                          <  1:240524/140‑1 (MQ=255)                                                          aCCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGa                                                                       >  1:1141926/1‑140 (MQ=255)                                                          aCCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGa                                                                       <  2:1141926/140‑1 (MQ=255)                                                                  cacCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATcc                                                               <  2:15496/140‑1 (MQ=255)                                                                     cAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCata                                                            >  1:669730/1‑140 (MQ=255)                                                                     cAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCata                                                            >  2:963879/1‑140 (MQ=255)                                                                       gggATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATg                                                          >  2:1123809/1‑140 (MQ=255)                                                                        ggATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGt                                                         >  2:1147593/1‑140 (MQ=255)                                                                            tttCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGa                                                     >  2:1110061/1‑140 (MQ=255)                                                                             ttCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGAc                                                    <  1:920740/140‑1 (MQ=255)                                                                             ttCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGAc                                                    >  1:921339/1‑140 (MQ=255)                                                                             ttCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGAc                                                    >  2:894576/1‑140 (MQ=255)                                                                                   cccTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGc                                              <  1:260001/140‑1 (MQ=255)                                                                                   cccTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGc                                              <  1:963879/140‑1 (MQ=255)                                                                                    ccTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCa                                             <  1:925934/140‑1 (MQ=255)                                                                                       tgtgtgt‑‑gCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACACAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACt                                          >  1:664162/1‑140 (MQ=255)                                                                                           tgt‑‑gCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTg                                      <  1:894576/140‑1 (MQ=255)                                                                                           tgt‑‑gCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTg                                      <  2:706270/140‑1 (MQ=255)                                                                                            gt‑‑gCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTgg                                     >  2:1166137/1‑140 (MQ=255)                                                                                                           cccGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTg                                      >  1:317202/1‑126 (MQ=255)                                                                                                           cccGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTg                                      <  2:317202/126‑1 (MQ=255)                                                                                                           cccGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTAt                             >  1:17547/1‑135 (MQ=255)                                                                                                           cccGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTAt                             <  2:17547/135‑1 (MQ=255)                                                                                                           cccGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTcccc                        <  1:1110061/140‑1 (MQ=255)                                                                                                           cccGGTGGTTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTcccc                        >  1:12077/1‑140 (MQ=255)                                                                                                                tggtTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACg                   >  1:300285/1‑140 (MQ=255)                                                                                                                tggtTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACACAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACg                   <  2:664162/140‑1 (MQ=255)                                                                                                                  gtTGAGCATCGACTGATTCCACCAGTCGCCAATGGTGATAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGGCCGACAAGGGCAACTTTTGGACCTTTATTCCCCTGACggg                 >  2:1074286/1‑140 (MQ=255)                                                                                                                    tGAGCATCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGca               >  1:552661/1‑140 (MQ=255)                                                                                                                          tCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAg         >  1:979463/1‑140 (MQ=255)                                                                                                                          tCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAg         <  2:552661/140‑1 (MQ=255)                                                                                                                          tCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAg         <  2:361655/140‑1 (MQ=255)                                                                                                                          tCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAg         >  2:1175324/1‑140 (MQ=255)                                                                                                                          tCGACTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAg         <  1:855113/140‑1 (MQ=255)                                                                                                                             aCTGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAGCTg      >  1:239123/1‑140 (MQ=255)                                                                                                                               tGATTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAGCTGcc    <  2:979463/140‑1 (MQ=255)                                                                                                                                 aTTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAGCTGCCCa  >  2:315625/1‑140 (MQ=255)                                                                                                                                 aTTCCACCAGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAGCTGCCCa  >  2:150217/1‑140 (MQ=255)                                                                                                                                        caGTCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACaaa          <  2:1014552/125‑1 (MQ=255)                                                                                                                                           |                                                                                                                                  ATCAAATGACAGCGGTTGTTGCTTTTGCTTTTCCGTTAACGCATGAGGCGCGACGGAACCAATCACACCAGGGATTTTCACTCCCTTGTGTGT‑‑GCGTATGGTTACCCGGTGGTTGAGCATCGACTGATTCCACCAGCCGCCAATGGTGGTAAAACGCAGAAAACCGCTCTCGTCGATGTGGGTGACCATAAAGCCGACTTCATCCATATGTCCGACAACGGCAACTTTTGGACCTTTATTCCCCTTACGGGCAACAAAGCTGCCCA  >  NC_000913/4090536‑4090801

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