Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,801,794	T→C	E259E (GAA→GAG)	thrS ←	threonine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,801,794	0	T	C	100.0%	161.4 / NA	45	E259E (GAA→GAG)	thrS	threonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/19);  total (26/19)

ACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAA  >  NC_000913/1801666‑1801928                                                                                                                                 |                                                                                                                                       acacGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCat                                                                                                                             <  1:650830/140‑1 (MQ=255) acacGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCat                                                                                                                             <  1:323983/140‑1 (MQ=255)       tCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTa                                                                                                                       >  2:557152/1‑140 (MQ=255)          atcatGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCt                                                                                                                                 <  2:442824/127‑1 (MQ=255)          atcatGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCt                                                                                                                                 >  1:442824/1‑127 (MQ=255)          atcatGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAg                                                                                                                    >  1:463377/1‑140 (MQ=255)                 aCGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGc                                                                                                             >  1:8930/1‑140 (MQ=255)                  cGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCt                                                                                                            >  2:771891/1‑140 (MQ=255)                        tttAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCATCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTAc                                                                                                      <  2:986049/140‑1 (MQ=255)                          tAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCg                                                                                                    >  2:160130/1‑140 (MQ=255)                               tCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGAtttt                                                                                               >  2:591992/1‑140 (MQ=255)                                   gATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACgg                                                                                           >  2:767390/1‑140 (MQ=255)                                         ggTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTcg                                                                                     <  1:160130/140‑1 (MQ=255)                                          gTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTcgc                                                                                    >  1:1078429/1‑140 (MQ=255)                                                tttCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTcg                                                                              >  1:966416/1‑140 (MQ=255)                                                  tCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTcgcg                                                                            >  2:1027592/1‑140 (MQ=255)                                                         aGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGcttctt                                                                     >  2:807241/1‑140 (MQ=255)                                                          gaacgaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTcgcg                                                                                   >  2:97235/1‑125 (MQ=255)                                                          gaacgaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTcgcg                                                                                   <  1:97235/125‑1 (MQ=255)                                                              gaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAgg                                                                <  2:440371/140‑1 (MQ=255)                                                              gaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAgg                                                                <  2:1078429/140‑1 (MQ=255)                                                              gaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAgg                                                                <  1:1027592/140‑1 (MQ=255)                                                              gaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAgg                                                                <  1:771891/140‑1 (MQ=255)                                                              gaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAgg                                                                <  1:767390/140‑1 (MQ=255)                                                                       cTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAgg                                                       <  2:841756/140‑1 (MQ=255)                                                                         tCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTa                                                     <  1:58842/140‑1 (MQ=255)                                                                          ccAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTaa                                                    >  1:311347/1‑140 (MQ=255)                                                                                  cGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAg                                            >  2:994878/1‑140 (MQ=255)                                                                                        aTGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCttt                                      >  2:400864/1‑140 (MQ=255)                                                                                        aTGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCttt                                      >  2:1019739/1‑140 (MQ=255)                                                                                        aTGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCttt                                      >  2:415392/1‑140 (MQ=255)                                                                                            tCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTg                                  >  1:243513/1‑140 (MQ=255)                                                                                                 ccGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGc                             >  2:786663/1‑140 (MQ=255)                                                                                                 ccGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGc                             >  2:829806/1‑140 (MQ=255)                                                                                                 ccGTCGTTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGc                             >  2:937083/1‑140 (MQ=255)                                                                                                      gtTGTGCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACg                        >  2:143655/1‑140 (MQ=255)                                                                                                         gtgCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGC‑TTTTTGTCTGCCCACGccgt                    <  1:1019739/140‑1 (MQ=255)                                                                                                          tgCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgt                    <  2:966416/140‑1 (MQ=255)                                                                                                          tgCCAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgt                    <  1:1027882/140‑1 (MQ=255)                                                                                                            ccAGAATACCATACCCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtac                  >  1:444089/1‑140 (MQ=255)                                                                                                                    ccataccCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATa          <  1:143655/140‑1 (MQ=255)                                                                                                                    ccataccCGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATa          <  2:902575/140‑1 (MQ=255)                                                                                                                          ccGGCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGc    <  1:197447/140‑1 (MQ=255)                                                                                                                            ggCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCaa  >  1:1056278/1‑140 (MQ=255)                                                                                                                            ggCGCCTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGATTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCaa  >  2:196114/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                       ACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAA  >  NC_000913/1801666‑1801928

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: ‑