Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,011,270	C→T	T436M (ACG→ATG)	uup →	ATP‑binding protein with possible role in replication

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,011,270	0	C	T	100.0%	93.4 / NA	28	T436M (ACG→ATG)	uup	ATP‑binding protein with possible role in replication
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/13);  total (15/13)

CGATAAAACGGTGATGGATAACCTTGCCGAAGGTAAGCAAGAGGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGACGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGAAACGCTGGA  >  NC_000913/1011151‑1011393                                                                                                                        |                                                                                                                            cGATAAAACGGTGATGGATAACCTTGCCGAAGGTAAGCAAGAGGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTgcg                                                                                                                >  2:156022/1‑133 (MQ=255) cGATAAAACGGTGATGGATAACCTTGCCGAAGGTAAGCAAGAGGTGATGGTTAACAGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTgcg                                                                                                                >  2:150286/1‑133 (MQ=255)           gTGATGGATAACCTTGCCGAAGGTAAGCAAGAGGTGATGGTTAACGGCAAGTCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCg                                                                                                      <  2:23575/133‑1 (MQ=255)                      ccTTGCCGAAGGTAAGCAAGAGGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAAc                                                                                           >  1:24253/1‑133 (MQ=255)                           ccGAAGGTAAGCAAGAGGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCtt                                                                                      >  1:451879/1‑133 (MQ=255)                                   aaGCAAGAGGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTgg                                                                              <  1:156022/133‑1 (MQ=255)                                          agGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGttt                                                                       >  2:547558/1‑133 (MQ=255)                                           ggTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTg                                                                      <  1:478585/133‑1 (MQ=255)                                              gATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTc                                                                   >  1:210036/1‑133 (MQ=255)                                               aTGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTcc                                                                  <  1:20808/133‑1 (MQ=255)                                               aTGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTcc                                                                  <  1:601180/133‑1 (MQ=255)                                                tGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCt                                                                 >  1:542447/1‑133 (MQ=255)                                                     aaCGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAAc                                                            >  1:54019/1‑133 (MQ=255)                                                      aCGGCAAGCCACGCCGCGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAAcc                                                           >  1:573645/1‑133 (MQ=255)                                                          cAAGCCACGCCACGTATTGGGCTACTTGAAGGACTTTCTGTCCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGc                                                       <  2:232313/133‑1 (MQ=255)                                                           aaGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCa                                                      <  2:359626/133‑1 (MQ=255)                                                            aGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCaa                                                     >  2:359708/1‑133 (MQ=255)                                                            aGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCaa                                                     <  1:547558/133‑1 (MQ=255)                                                                     acgTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGAt                                            <  1:359708/133‑1 (MQ=255)                                                                               cTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGaa                                  <  1:150286/133‑1 (MQ=255)                                                                                  tttGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGTGCGTTTGTTCCTCAAACCAAGCAGCTTCTTAATTCTTCACTAAcct                               >  1:222461/1‑132 (MQ=255)                                                                                  tttGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAGCCg                               >  2:38682/1‑133 (MQ=255)                                                                                  tttGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCg                               <  2:542447/133‑1 (MQ=255)                                                                                            tttCTGTTCCATCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATc                     >  2:570789/1‑133 (MQ=255)                                                                                                      aTCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGa           <  2:24253/133‑1 (MQ=255)                                                                                                       tCCGAAACGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGaa          >  1:373930/1‑133 (MQ=255)                                                                                                            aaCGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGAAACGCt     >  1:569946/1‑133 (MQ=255)                                                                                                              cGGGCGATGATGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGAAACGCTgg   <  2:373930/133‑1 (MQ=255)                                                                                                               gggCGATGACGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGAAACGCTGGa  <  1:569091/133‑1 (MQ=255)                                                                                                                        |                                                                                                                            CGATAAAACGGTGATGGATAACCTTGCCGAAGGTAAGCAAGAGGTGATGGTTAACGGCAAGCCACGCCACGTATTGGGCTATTTGCAGGACTTTCTGTTCCATCCGAAACGGGCGATGACGCCGGTACGTGCGCTTTCTGGCGGTGAGCGGAACCGCTTGCTGCTGGCGCGTTTGTTCCTCAAACCAAGCAACTTATTGATTCTTGACGAACCGACCAACGATCTTGATGTCGAAACGCTGGA  >  NC_000913/1011151‑1011393

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