Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,804,371	T→C	pseudogene (1278/1900 nt)	arpB →	putative ankyrin repeat protein B, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,804,371	0	T	C	100.0%	116.1 / NA	33	pseudogene (1278/1900 nt)	arpB	putative ankyrin repeat protein B, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/18);  total (15/18)

CTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGTAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGAAATAATCGATATCTTTCTGA  >  NC_000913/1804245‑1804496                                                                                                                               |                                                                                                                              cTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCAcatc                                                                                                                  >  1:346114/1‑140 (MQ=255)  tGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCAcatca                                                                                                                 <  2:346114/140‑1 (MQ=255)      aaaGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACAtcattca                                                                                                             <  2:368304/140‑1 (MQ=255)       aaGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACAtcattcat                                                                                                            <  2:404903/140‑1 (MQ=255)                  gAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACTGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCt                                                                                                 <  2:486914/140‑1 (MQ=255)                     ataataACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGt                                                                                              <  1:363901/140‑1 (MQ=255)                               cgcAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACta                                                                                    <  1:17408/140‑2 (MQ=255)                               cgcAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTg                                                                                    <  1:2349/140‑1 (MQ=255)                               cgcAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTg                                                                                    >  2:2349/1‑140 (MQ=255)                               cgcAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTg                                                                                    >  2:230289/1‑140 (MQ=255)                                     aTAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACt                                                                              <  1:234013/140‑1 (MQ=255)                                              gTACACCTGCTTTATTTTTGGCTATGATGAATGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCt                                                                     >  2:42473/1‑140 (MQ=255)                                              gTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTaa                                                                            >  2:319485/1‑133 (MQ=255)                                              gTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTaa                                                                            <  1:319485/133‑1 (MQ=255)                                              gTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCt                                                                     >  2:201310/1‑140 (MQ=255)                                                   ccTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGa                                                                <  1:327794/140‑1 (MQ=255)                                                     cctgctATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGaaa                                                              >  1:494263/6‑140 (MQ=255)                                                       cTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAaca                                                            >  1:517170/1‑140 (MQ=255)                                                             ttttGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGa                                                      <  1:313713/140‑1 (MQ=255)                                                             ttttGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGa                                                      >  2:259956/1‑140 (MQ=255)                                                                    tATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTtata                                               <  2:517170/140‑1 (MQ=255)                                                                         tgaACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTcc                                          <  1:245686/140‑1 (MQ=255)                                                                            aCGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCAtg                                       <  1:126090/140‑1 (MQ=255)                                                                             cGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCAtgt                                      >  2:193116/1‑140 (MQ=255)                                                                                 cATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTa                                  <  1:193116/140‑1 (MQ=255)                                                                                          aaCGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGa                         <  2:235988/140‑1 (MQ=255)                                                                                          aaCGTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGa                         >  1:61674/1‑140 (MQ=255)                                                                                             gTGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGa                      >  1:1987/1‑140 (MQ=255)                                                                                              tGAAAATATTTATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGaa                     >  2:168039/1‑140 (MQ=255)                                                                                                        tattcAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGAAATAATCGata           >  2:127574/5‑140 (MQ=255)                                                                                                        tATGCAAGAAATTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGAAATAATCGata           >  2:1859/1‑140 (MQ=255)                                                                                                                 aaTTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGAAATAATCGATATCTTTCTGa  <  1:91229/140‑1 (MQ=255)                                                                                                                 aaTTCAGTCACTGGCAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGAAATAATCGATATCTTTCTGa  <  2:76070/140‑1 (MQ=255)                                                                                                                               |                                                                                                                              CTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATACTGACAACGTGAAAATATTTATGCAAGAAATTCAGTCACTGGTAGATAATCACATCATTCATGAAGATAATCTGGTTAAATTACTGCAAACTAAAAGTGCTAACGAAACACCTGGACTTTATATCTCCATGTTGTATGGATTCGATGAAATAATCGATATCTTTCTGA  >  NC_000913/1804245‑1804496

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