Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,789,783	G→A	E267K (GAA→AAA)	csiD →	PF08943 family protein CsiD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,789,783	0	G	A	100.0%	152.8 / NA	45	E267K (GAA→AAA)	csiD	PF08943 family protein CsiD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/16);  total (29/16)

AGCAAAGATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTGAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACTGATGCGCCAGCG  >  NC_000913/2789651‑2789916                                                                                                                                     |                                                                                                                                      aGCAAAGATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAg                                                                                                                                >  1:407037/1‑140 (MQ=255)     aaGATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCaaa                                                                                                                            <  1:590157/140‑1 (MQ=255)     aaGATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCaaa                                                                                                                            <  1:778642/140‑1 (MQ=255)       gATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAgg                                                                                                                          <  2:523815/140‑1 (MQ=255)       gATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAgg                                                                                                                          <  2:104342/140‑1 (MQ=255)        aTGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGc                                                                                                                         >  2:211878/1‑140 (MQ=255)         tGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCa                                                                                                                        <  1:56031/140‑1 (MQ=255)           tttttCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCAtt                                                                                                                      <  2:9816/140‑1 (MQ=255)             tttCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCAttct                                                                                                                    >  2:741848/1‑140 (MQ=255)                cATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCAttctttc                                                                                                                 >  2:576485/1‑140 (MQ=255)                   ccGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCtgt                                                                                                              >  2:522186/1‑140 (MQ=255)                   ccGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCtgg                                                                                                              >  2:449205/1‑139 (MQ=255)                         ttCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGcccgt                                                                                                        >  2:831476/1‑140 (MQ=255)                           cGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGcccgttc                                                                                                      >  2:409018/1‑140 (MQ=255)                               gTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGcccgttcccgt                                                                                                  <  1:222667/140‑1 (MQ=255)                                 cGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTg                                                                                                <  2:991626/140‑1 (MQ=255)                                     cAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCaa                                                                                            <  1:170209/140‑1 (MQ=255)                                              cGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGtt                                                                                   >  2:503775/1‑140 (MQ=255)                                                cccGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAGGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGa                                                                                 >  1:71114/1‑140 (MQ=255)                                                cccGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGa                                                                                 >  1:550722/1‑140 (MQ=255)                                                cccGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGa                                                                                 >  2:144744/1‑140 (MQ=255)                                                               tatCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCt                                                                  >  2:63132/1‑140 (MQ=255)                                                               tatCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCt                                                                  >  1:526927/1‑140 (MQ=255)                                                               tatCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGGGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCt                                                                  >  1:59926/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGCCGCCATTAAAACCAGCAAAGGCATTCTTTCTTTGCCCGTTCCCGTTGGCAAAATCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  1:687658/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  1:300760/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  1:380307/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  1:27716/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  1:80430/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  1:157686/1‑140 (MQ=255)                                                                    aCCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTg                                                             >  2:887579/1‑140 (MQ=255)                                                                         ttCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACgg                                                        >  1:836871/1‑140 (MQ=255)                                                                         ttCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACgg                                                        >  2:478214/1‑140 (MQ=255)                                                                         ttCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACgg                                                        >  1:670505/1‑140 (MQ=255)                                                                           cGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTc                                                      >  1:1023453/1‑140 (MQ=255)                                                                                          cTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTc                                       <  1:741848/140‑1 (MQ=255)                                                                                                agaagGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCAcc                                 <  1:503775/140‑1 (MQ=255)                                                                                                agaagGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCAcc                                 <  1:144744/140‑1 (MQ=255)                                                                                                   agGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCgg                              >  2:722612/1‑140 (MQ=255)                                                                                                          tgGTTGAGCGAGCTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTgcg                       <  1:63132/140‑1 (MQ=255)                                                                                                                   gagcTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACt              >  1:346125/1‑140 (MQ=255)                                                                                                                    agcTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACTg             <  1:887579/140‑1 (MQ=255)                                                                                                                    agcTTTCCGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACTg             >  1:671161/1‑140 (MQ=255)                                                                                                                           cGACGCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACTGATGCGcc      <  2:380307/140‑1 (MQ=255)                                                                                                                               gCCATTAAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACTGATGCGCCAGCg  <  2:670505/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      AGCAAAGATGTTTTTCATCCGGTGTTCGACGTCGATCAACAGGGTCGCCCGGTGATGCGCTATATCGACCAGTTCGTCCAGCCGAAAGACTTCGAAGAAGGCGTGTGGTTGAGCGAGCTTTCCGACGCCATTGAAACCAGCAAAGGCATTCTTTCTGTGCCCGTTCCCGTTGGCAAATTCCTGTTGATTAACAACCTGTTCTGGCTGCACGGTCGCGACCGCTTTACTCCGCACCCGGATCTGCGCCGTGAACTGATGCGCCAGCG  >  NC_000913/2789651‑2789916

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