Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,803,567	G→A	pseudogene (50/1476 nt)	arpB →	pseudogene, ankyrin repeats

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,803,567	0	G	A	100.0%	153.5 / NA	49	pseudogene (50/1476 nt)	arpB	pseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/25);  total (24/25)

TAACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATA  >  NC_000913/1803432‑1803700                                                                                                                                        |                                                                                                                                      tAACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGc                                                                                                                                    <  1:1275647/139‑1 (MQ=255)   aCCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCAt                                                                                                                                  >  2:98789/1‑139 (MQ=255)     cATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAATAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGt                                                                                                                                <  1:807099/139‑1 (MQ=255)                          tataTTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGAc                                                                                                           >  1:828807/1‑139 (MQ=255)                             atTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCt                                                                                                        >  2:494420/1‑139 (MQ=255)                                  acaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTAtttt                                                                                                   <  1:706213/139‑1 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                               >  2:182167/1‑139 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                               >  1:500270/1‑139 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                               >  2:623336/1‑139 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                               <  2:828807/139‑1 (MQ=255)                                        acagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCATCAGTAAATATCAATTGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCg                                                                                             <  2:760571/139‑1 (MQ=255)                                          agTTAACTGCTGAAGATCGTAACCATTGTATGAAACTATCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                           <  1:623336/139‑1 (MQ=255)                                           gTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                          >  1:761180/1‑139 (MQ=255)                                            ttAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCAt                                                                                         <  2:1176425/139‑1 (MQ=255)                                                    tgAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAtt                                                                                 >  2:203005/1‑139 (MQ=255)                                                      aaGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                          >  1:34860/1‑112 (MQ=255)                                                      aaGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                          <  2:34860/112‑1 (MQ=255)                                                             tAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGc                                                                        <  1:203005/139‑1 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                                     >  2:167485/1‑118 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                                     <  1:167485/118‑1 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATa                                                                           >  2:544531/1‑128 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATa                                                                           <  1:544531/128‑1 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGAt                                                                 >  2:1112650/1‑138 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGAt                                                                 <  1:1112650/138‑1 (MQ=255)                                                                         gAAACTTTCTAAAAAAATCCATAGGGAAAAATACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCAtt                                                            <  1:638112/139‑1 (MQ=255)                                                                            aCTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                           <  1:1173002/105‑1 (MQ=255)                                                                            aCTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                           >  2:1173002/1‑105 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGc                                                                                                                                    >  2:355369/1‑62 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGc                                                                                                                                    <  1:355369/62‑1 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                                       <  2:500270/139‑1 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                                     >  2:738825/1‑139 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                                     >  2:386728/1‑139 (MQ=255)                                                                                     aaaaaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGaaa                                                >  2:398719/1‑139 (MQ=255)                                                                                        aaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAgg                                             <  1:738825/139‑1 (MQ=255)                                                                                        aaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAgg                                             <  1:494420/139‑1 (MQ=255)                                                                                                         aCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTa                            <  1:182167/139‑1 (MQ=255)                                                                                                                acCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTa                     <  2:1092027/139‑1 (MQ=255)                                                                                                                    aaCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCaa                                    >  2:1294023/1‑120 (MQ=255)                                                                                                                    aaCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCaa                                    <  1:1294023/120‑1 (MQ=255)                                                                                                                     aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTc                                                                                              <  2:455387/61‑1 (MQ=255)                                                                                                                     aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTc                                                                                              >  1:455387/1‑61 (MQ=255)                                                                                                                     aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcct                >  2:136228/1‑139 (MQ=255)                                                                                                                     aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTCGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAGACATTActtct                <  2:761180/139‑1 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt            >  2:636296/1‑139 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt            <  1:98789/139‑1 (MQ=255)                                                                                                                                 aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTc                                                                                    <  1:490267/59‑1 (MQ=255)                                                                                                                                 aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTc                                                                                    >  2:490267/1‑59 (MQ=255)                                                                                                                                 aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCa    >  1:1258417/1‑139 (MQ=255)                                                                                                                                   tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATa  >  1:916094/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                      TAACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATA  >  NC_000913/1803432‑1803700

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