Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,804,989	T→G	pseudogene (1472/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,804,989	0	T	G	100.0%	216.5 / NA	56	pseudogene (1472/1476 nt)	arpB	pseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/29);  total (27/29)

TTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATC  >  NC_000913/1804844‑1805136                                                                                                                                                  |                                                                                                                                                    ttGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAAc                                                                                                                                                >  1:397985/1‑151 (MQ=255)     ccGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGc                                                                                                                                            <  1:30793/151‑1 (MQ=255)     ccGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGc                                                                                                                                            >  2:446213/1‑151 (MQ=255)       gCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCt                                                                                                                                           >  1:501030/1‑150 (MQ=255)            aaaTCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTAc                                                                                                                                     <  1:24242/151‑1 (MQ=255)             aaTCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACg                                                                                                                                    >  2:141880/1‑151 (MQ=255)               tCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGat                                                                                                                                  >  2:124605/1‑151 (MQ=255)                cATTACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTAGAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGat                                                                                                                                  <  2:334928/150‑1 (MQ=255)                cATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGat                                                                                                                                  >  2:390783/1‑150 (MQ=255)                   gACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGatatat                                                                                                                              >  1:736815/1‑151 (MQ=255)                      aacaTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGatatatat                                                                                                                            <  2:254336/150‑1 (MQ=255)                      aacaTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATc                                                                                                                           <  1:803349/151‑1 (MQ=255)                        caTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGa                                                                                                                         <  1:411184/151‑1 (MQ=255)                         aTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGAt                                                                                                                        >  2:867479/1‑151 (MQ=255)                              aGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATa                                                                                                                   <  1:66169/151‑1 (MQ=255)                              aGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATa                                                                                                                   <  2:308656/151‑1 (MQ=255)                               gCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATaa                                                                                                                  <  2:108533/151‑1 (MQ=255)                                cTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAat                                                                                                                 >  2:704604/1‑151 (MQ=255)                                 tGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAatt                                                                                                                >  2:191750/1‑151 (MQ=255)                                   ttCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAatttt                                                                                                              >  2:179704/1‑151 (MQ=255)                                   ttCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAat                                                                                                                 <  2:96480/148‑1 (MQ=255)                                      atatAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAattttatt                                                                                                           <  2:232976/151‑1 (MQ=255)                                      atatAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAattttatt                                                                                                           >  2:590804/1‑151 (MQ=255)                                           gCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTAttttttt                                                                                                      <  1:84180/151‑1 (MQ=255)                                            ccATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTc                                                                                                     <  1:387233/151‑1 (MQ=255)                                            ccATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTc                                                                                                     <  1:649345/151‑1 (MQ=255)                                                tcatcaTAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATa                                                                                                  >  2:898677/1‑150 (MQ=255)                                                tcatcaTAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATa                                                                                                  <  2:216969/150‑1 (MQ=255)                                                    catAAGCATTATAAAACTGTAGAAACATATTATGCTTCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAAc                                                                                              <  1:506939/150‑1 (MQ=255)                                                        aGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGa                                                                                         <  1:550042/151‑1 (MQ=255)                                                             tataAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAAc                                                                                      >  1:827993/1‑149 (MQ=255)                                                              ataAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAg                                                                                    >  1:370007/1‑150 (MQ=255)                                                              ataAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGt                                                                                   <  2:399069/151‑1 (MQ=255)                                                                aaaaCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTg                                                                                  <  1:264214/150‑1 (MQ=255)                                                                 aaaCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGAc                                                                                >  1:500360/1‑151 (MQ=255)                                                                  aaCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGAcc                                                                               <  2:876304/151‑1 (MQ=255)                                                                  aaCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGAcc                                                                               >  1:260659/1‑151 (MQ=255)                                                                     tGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCAtg                                                                            >  1:32326/1‑151 (MQ=255)                                                                             cATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAAc                                                                    <  1:918163/151‑1 (MQ=255)                                                                                attatGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGc                                                                  <  1:344574/150‑1 (MQ=255)                                                                                          tATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGc                                                       >  1:329635/1‑151 (MQ=255)                                                                                              aaTGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGt                                                    <  2:211423/150‑1 (MQ=255)                                                                                               aTGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGtaa                                                  <  1:141587/151‑1 (MQ=255)                                                                                                  cAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGtaata                                                >  1:406241/1‑150 (MQ=255)                                                                                                    aTCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAAc                                              >  1:271212/1‑150 (MQ=255)                                                                                                     tCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACa                                             >  2:259248/1‑150 (MQ=255)                                                                                                          cAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAAt                                       >  2:829021/1‑151 (MQ=255)                                                                                                               ccTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCt                                   <  2:290107/150‑1 (MQ=255)                                                                                                                  gAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTg                                 <  1:186911/149‑1 (MQ=255)                                                                                                                  gAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGAt                               >  2:798383/1‑151 (MQ=255)                                                                                                                        tAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGt                         >  2:87559/1‑151 (MQ=255)                                                                                                                                 tgaAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGt                >  1:209471/1‑151 (MQ=255)                                                                                                                                           aCGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGt      <  2:277272/151‑1 (MQ=255)                                                                                                                                              tATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTa    <  1:446223/150‑1 (MQ=255)                                                                                                                                              tATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTAt   >  1:729486/1‑151 (MQ=255)                                                                                                                                               aTTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATc  <  2:225452/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TTGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTATC  >  NC_000913/1804844‑1805136

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