Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,026,152	A→G	E561G (GAG→GGG)	helD →	DNA helicase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,026,152	0	A	G	100.0%	115.9 / NA	34	E561G (GAG→GGG)	helD	DNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/16);  total (18/16)

TTCAGCAGAACCCAGGCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATT  >  NC_000913/1026017‑1026280                                                                                                                                        |                                                                                                                                 ttCAGCAGAACCCAGGCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGgcgc                                                                                                                              <  2:75924/140‑1 (MQ=255)       agAACCCAGGCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTg                                                                                                                        >  1:347714/1‑140 (MQ=255)              aGGCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTgg                                                                                                                 <  1:502500/140‑1 (MQ=255)              aGGCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTgg                                                                                                                 <  1:98334/140‑1 (MQ=255)                gCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGgcg                                                                                                               >  2:616368/1‑140 (MQ=255)                     cTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTa                                                                                                          >  1:313951/1‑140 (MQ=255)                     cTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTa                                                                                                          >  2:322620/1‑140 (MQ=255)                     cTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTa                                                                                                          >  2:279958/1‑140 (MQ=255)                       gAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTAcc                                                                                                        <  2:305268/140‑1 (MQ=255)                        aaaaaGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCa                                                                                                       <  2:347714/140‑1 (MQ=255)                         aaaaGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCAt                                                                                                      <  2:632216/140‑1 (MQ=255)                               cGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATg                                                                                                >  2:257258/1‑140 (MQ=255)                                      cAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTg                                                                                         <  1:616368/140‑1 (MQ=255)                                       aGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGc                                                                                        >  1:413619/1‑140 (MQ=255)                                                                   tCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTgg                                                            >  2:210344/1‑140 (MQ=255)                                                                     aCGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGTGGCaa                                                                    >  1:358346/1‑130 (MQ=255)                                                                     aCGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGTGGCaa                                                                    <  2:358346/130‑1 (MQ=255)                                                                               aTGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCaa                                                >  1:48839/1‑140 (MQ=255)                                                                                 gagaGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAAt                                              >  2:390554/1‑140 (MQ=255)                                                                                   gagTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                            <  1:48447/140‑1 (MQ=255)                                                                                    agTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATc                                             <  2:603754/138‑1 (MQ=255)                                                                                    agTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATc                                             >  1:603754/1‑138 (MQ=255)                                                                                      tCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACt                                         >  1:276942/1‑140 (MQ=255)                                                                                        aaCTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                       >  1:270762/1‑140 (MQ=255)                                                                                         aCTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTa                                      >  2:151973/1‑140 (MQ=255)                                                                                                             aGCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaaa                  <  1:257258/140‑1 (MQ=255)                                                                                                              gCTCTCTGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAg                 >  2:716103/1‑140 (MQ=255)                                                                                                                tctctGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAggg               >  1:645515/1‑140 (MQ=255)                                                                                                                 ctctGGTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGc              <  2:48839/140‑1 (MQ=255)                                                                                                                     ggTTATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACa          <  2:621397/140‑1 (MQ=255)                                                                                                                       ttATGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAgg        <  2:313951/140‑1 (MQ=255)                                                                                                                         aTGCCAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCg      >  1:528732/1‑140 (MQ=255)                                                                                                                            ccAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAt   <  1:716103/140‑1 (MQ=255)                                                                                                                             cAAACCGGAAGGGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt  <  2:528732/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 TTCAGCAGAACCCAGGCCAGCTGAAAAAGCCGCTAAACAGCTTAACCAATGGAGACAAAAAAGCCGTCACGTTATTGGATGAGAGTCAACTTGACGCTTTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATT  >  NC_000913/1026017‑1026280

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