Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,477,822	G→A	A61A (GCG→GCA)	slyX →	protein SlyX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,477,822	0	G	A	100.0%	65.5 / NA	20	A61A (GCG→GCA)	slyX	protein SlyX
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/2);  total (18/2)

CGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCC  >  NC_000913/3477685‑3477957                                                                                                                                          |                                                                                                                                        cGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATc                                                                                                                                       >  1:1560894/1‑140 (MQ=255)  gACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAAcatcgcatcg                                                                                                                                      <  2:470605/140‑1 (MQ=255)  gACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAAcatcgcatcg                                                                                                                                      <  2:258096/140‑1 (MQ=255)  gACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAAcatcgcatcg                                                                                                                                      <  1:347384/140‑1 (MQ=255)   aCTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACatcgcatcgc                                                                                                                                     >  1:785756/1‑140 (MQ=255)    cTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACAtcgcatcgca                                                                                                                                    <  2:576692/140‑1 (MQ=255)        cTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCt                                                                                                                                >  1:117798/1‑140 (MQ=255)                gAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAAc                                                                                                                        <  1:1033291/140‑1 (MQ=255)                     cACCATTGAAGAACTGAACGTCAGGGGGACCGCTAAGAAAATGGAGATGGCAAAACAGCGCAATCATCTGCTTCTGCGGACCGAAAAATAAAAACCCGCCCAGCCGTCAAACATCGCATCGCAGGCGAAAAAAAcgcaac                                                                                                                   <  2:785756/140‑1 (MQ=255)                          ttGAAGAGCGGAACGTCACGTTGCCCGCTCATGAAAGGGAAAGGGCGAAACGGCGCGATCATCGGCGTTTGCGGACGGAAAAGTAAAAAGCCAGCCAGCCGTAGAACATCGCACGGCGGGCTGAAGAAAcgccaccgcca                                                                                                              <  2:951807/140‑1 (MQ=255)                          ttGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAAcgccaccgcca                                                                                                              >  2:1149456/1‑140 (MQ=255)                           tGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAAcgccaccgccac                                                                                                             >  1:1449948/1‑140 (MQ=255)                             aagaaCTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACat                                                                                                           >  1:563597/1‑140 (MQ=255)                                 aCTGAACGTCACGGTGCCCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACattatt                                                                                                       <  2:880137/140‑1 (MQ=255)                                     aaCGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAgg                                                                                                   >  1:174592/1‑140 (MQ=255)                                     aaCGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAgg                                                                                                   >  1:1131691/1‑140 (MQ=255)                                          cACGGTGACCGCTCATAAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTaa                                                                                              >  1:23724/1‑140 (MQ=255)                                               tGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTaaaaaaa                                                                                         >  1:1620718/1‑140 (MQ=255)                                                        aTGAAATGGAGAGGGCGAAACTGCGCGATCATCTGCGTCTGCGGACGGAGAAGTAAAAACCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAACCCCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAtt                                                                                <  2:910616/140‑1 (MQ=255)                                                         tGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTc                                                                               >  1:548380/1‑140 (MQ=255)                                                           aaaTGGAGAGGGCAAAACTGCGCGATCATCTGCGTCTGTGGACCGAGAAGTAAAAAGCCAGCCAGCCGTCAAACACCGCACCGAAGGCTGAAAAACCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTccc                                                                             <  2:120441/140‑1 (MQ=255)                                                                      gCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGt                                                                  >  2:20867/1‑140 (MQ=255)                                                                          aaCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACt                                                              >  1:876720/1‑140 (MQ=255)                                                                            cTGCGCGATCCTCTGCGTCTGCTGACCGAAAAGTAAAACCCCGGCCAGCCGTCGAACATCGCACCGCAGGCTGAAGAAACCCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACttt                                                            <  2:563597/140‑1 (MQ=255)                                                                                 cgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCgg                                                       >  1:196005/1‑140 (MQ=255)                                                                                         tGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTg                                               >  1:1575307/1‑140 (MQ=255)                                                                                              ctgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaac                                          >  1:1043461/1‑140 (MQ=255)                                                                                                        agaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgcc                                >  1:1231291/1‑140 (MQ=255)                                                                                                                               tCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc         >  1:1081308/1‑140 (MQ=255)                                                                                                                                 gAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCt       >  1:1268413/1‑140 (MQ=255)                                                                                                                                      tcgcatcgcaGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGcc  >  1:388049/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        CGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCC  >  NC_000913/3477685‑3477957

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