Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,290,753	G→A	C171Y (TGC→TAC)	rssB →	PcnB‑degradosome interaction factor; response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,290,753	0	G	A	100.0%	140.9 / NA	37	C171Y (TGC→TAC)	rssB	PcnB‑degradosome interaction factor; response regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/21);  total (16/21)

CTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGTTATTAC  >  NC_000913/1290613‑1290899                                                                                                                                                |                                                                                                                                                   ctATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAAt                                                                                                                                          >  1:554541/1‑151 (MQ=255)    ccccggattttttATTCCGCGGGTGGAGGAGGAGAAAGGGTTTTTG‑CGCGAGGGGGATGAAAGGT‑TTAATACCCC‑GGCGGCGGCGGCAAAATTATTACAGAAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTa                                                                                                                                            <  1:538346/136‑1 (MQ=255)      ccagattttttATTTCGCGTGTGGAGAAAGA‑GAAAGGTTTTTTGCGCGATGGGGATGCAATTT‑TTAATCCCCC‑TCCCGAAGGGGCAAATTTATTACGGGAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCg                                                                                                                                     <  2:502788/143‑1 (MQ=255)         gCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTc                                                                                                                                   >  1:78917/1‑150 (MQ=255)           aTGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAAt                                                                                                                                >  2:279724/1‑151 (MQ=255)                  aTTCGCGCGTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGc                                                                                                                         >  1:108646/1‑151 (MQ=255)                      ggggttggagaaggagaaaGGGTTTTTT‑CGCGA‑TGGGAATACATGG‑TTGAAAACCCGGGCCGCGGCGGCAAAATTATTAAGGAAATTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCcgcg                                                                                                                     <  2:232357/138‑1 (MQ=255)                      gggggtggaggagga‑gaaaGGTTTTTTGCGCGAGGGGAAGGAATTGG‑TTGACACCCC‑TGCCGAGGCGGATAATTTTTTACGGAAACAAAAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCcgc                                                                                                                      <  2:524260/143‑1 (MQ=255)                        gcttggaagaagagaaaATTTTTTTT‑CGCGGTTGGGTGGAATTGT‑TTAATCCCCC‑TGCGGCGGCGGAAAATTTTTAAAAGGAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCcgcg                                                                                                                     <  2:325206/141‑1 (MQ=255)                        gcgTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCgg                                                                                                                    >  2:245649/1‑150 (MQ=255)                                       aaggttttttg‑CGCGACTGGGATG‑AATGGTTTAATACCCC‑GCCCGCGGGGGCAAATTTATTACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGt                                                                                                    <  1:195610/145‑1 (MQ=255)                                             ttttg‑CGCGACGGGGATGCATTGG‑TTAAAACCCC‑TGCCGCGGCGGAAAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCtt                                                                                                <  1:264368/149‑1 (MQ=255)                                               ttt‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGat                                                                                             >  1:123793/1‑150 (MQ=255)                                                   cgcgACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTg                                                                                         >  2:109957/1‑151 (MQ=255)                                                    gcgatggggATGAATTGT‑TTGAAACCCC‑TGCCGCGGCGGAAAAATTATTAAAGAAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATAtt                                                                                          <  2:535382/143‑1 (MQ=255)                                                     agACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGcc                                                                                       >  2:287466/2‑151 (MQ=255)                                                      gACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCg                                                                                      >  1:652653/1‑151 (MQ=255)                                                       aCTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGc                                                                                     >  1:491832/1‑151 (MQ=255)                                                        cTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGc                                                                                     >  1:157350/1‑150 (MQ=255)                                                        cTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGc                                                                                     >  2:273397/1‑150 (MQ=255)                                                          gggATGAAATGG‑TTGATACCCC‑GGCCGCGGCGGCAAAATTATACCGGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACt                                                                                  <  1:262092/151‑1 (MQ=255)                                                               tgaattgtttGATACCCC‑TGCCGCAGGGGAAAAATTATAAAAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACttt                                                                                <  2:23984/141‑1 (MQ=255)                                                                  agtgtgtaaTACCCC‑GGCCGAAGCGGAGAATTTATGAAAGAATATAAAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACttt                                                                                <  1:201580/138‑1 (MQ=255)                                                                       ttaaTACCCC‑GGCCGCAGGGGGAAATTTATTAAAGAAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGaaaa                                                                         <  1:610715/145‑1 (MQ=255)                                                                            acccc‑cGCCGCGGCGGAAAAATTATTACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGaaa                                                                          <  2:203573/140‑1 (MQ=255)                                                                              ccccgcccGCGGCGGAAAAATTATAACAGAACCTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTg                                                                  <  1:321119/143‑1 (MQ=255)                                                                                             cgAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTg                                                  >  1:618835/1‑150 (MQ=255)                                                                                               aaaTTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATg                                               >  2:427166/1‑151 (MQ=255)                                                                                               aaaTTATTACAGGAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTg                                                  <  1:421269/148‑1 (MQ=255)                                                                                                  ttattaCAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTc                                             >  2:281261/1‑150 (MQ=255)                                                                                                  ttattaAGGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGt                                              <  2:652482/149‑1 (MQ=255)                                                                                                           aaatcccAACCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTg                                   <  2:8122/145‑1 (MQ=255)                                                                                                               acaacccgccgGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAGCCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGGTGTCACCCGAGCTGGACa                               <  2:67306/146‑1 (MQ=255)                                                                                                                  aaCCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATa                             <  1:115739/150‑1 (MQ=255)                                                                                                                  aaCCGCCGGTTCAGCAGGTGATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATa                             <  1:109959/150‑1 (MQ=255)                                                                                                                                     gATTTCCCATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTg         <  1:216242/151‑1 (MQ=255)                                                                                                                                            cATTACCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGTTATTAc  >  1:589775/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                   CTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTT‑CGCGACTGGGATGCAATGG‑TTGATAACCC‑TGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGTTATTAC  >  NC_000913/1290613‑1290899

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