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%	203.9 / NA	54	pseudogene (1472/1476 nt)	arpB	pseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/32);  total (22/32)

TGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTAT  >  NC_000913/1804845‑1805135                                                                                                                                                 |                                                                                                                                                   tGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAAc                                                                                                                                               <  1:111212/150‑1 (MQ=255)  ggCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCa                                                                                                                                             <  2:589560/151‑1 (MQ=255)         aaaaaTCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTAttt                                                                                                                                      <  1:621600/151‑1 (MQ=255)         aaaaaTCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTAt                                                                                                                                        <  2:266602/149‑1 (MQ=255)             aTCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGa                                                                                                                                  <  2:242848/151‑1 (MQ=255)                   acaacaTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGatatata                                                                                                                            <  2:520934/151‑1 (MQ=255)                        aTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGa                                                                                                                        <  1:140559/150‑1 (MQ=255)                        aTGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGAt                                                                                                                       <  2:450656/151‑1 (MQ=255)                             aGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACAt                                                                                                                   >  2:223126/1‑150 (MQ=255)                              gCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATaa                                                                                                                 <  2:55562/151‑1 (MQ=255)                               cTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAat                                                                                                                >  1:443749/1‑151 (MQ=255)                                tGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAatt                                                                                                               >  1:754521/1‑151 (MQ=255)                                 gTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAattt                                                                                                              <  2:776667/151‑1 (MQ=255)                                 gTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAattt                                                                                                              <  1:543446/151‑1 (MQ=255)                                 gTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAatt                                                                                                               <  2:477993/150‑1 (MQ=255)                                    cATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAattttat                                                                                                           <  1:637269/151‑1 (MQ=255)                                       ataGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAattttattt                                                                                                         <  1:365915/150‑1 (MQ=255)                                        taGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTAttttt                                                                                                       >  1:232100/1‑151 (MQ=255)                                        taGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTAttttt                                                                                                       >  1:254012/1‑151 (MQ=255)                                             aTTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCGAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCAt                                                                                                  <  2:740956/151‑1 (MQ=255)                                             aTTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCa                                                                                                   >  1:648800/1‑150 (MQ=255)                                                catcatAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATaaa                                                                                               <  1:261145/151‑1 (MQ=255)                                                    atAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGc                                                                                           <  1:477199/151‑1 (MQ=255)                                                       aGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGa                                                                                        <  1:679288/151‑1 (MQ=255)                                                        gCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGaa                                                                                       <  1:21283/151‑1 (MQ=255)                                                          aTTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGaaa                                                                                      <  1:730523/150‑1 (MQ=255)                                                             ataAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAg                                                                                   <  2:588443/150‑1 (MQ=255)                                                                 aaCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGAcc                                                                              >  2:508564/1‑151 (MQ=255)                                                                  aCTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCa                                                                             <  1:303743/151‑1 (MQ=255)                                                                            cATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAAc                                                                   >  2:415350/1‑151 (MQ=255)                                                                             atatTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACg                                                                  >  1:410472/1‑151 (MQ=255)                                                                                  atGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCAc                                                              <  2:99687/150‑1 (MQ=255)                                                                                    gctgctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCAccc                                                            >  2:328955/1‑150 (MQ=255)                                                                                     ctgctATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCg                                                           >  2:404012/1‑150 (MQ=255)                                                                                          aTTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGTACAACGCCACCCGACGGCt                                                     <  1:797057/151‑1 (MQ=255)                                                                                              aTGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGta                                                  >  2:623416/1‑150 (MQ=255)                                                                                              aTGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGta                                                  >  2:493009/1‑150 (MQ=255)                                                                                               tGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGtaat                                                >  1:678975/1‑151 (MQ=255)                                                                                                 cAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGtaata                                               >  2:277234/1‑150 (MQ=255)                                                                                                  aaTCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAAc                                             <  1:456546/151‑1 (MQ=255)                                                                                                     cAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACaa                                           <  2:554409/150‑1 (MQ=255)                                                                                                       gCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACaaaaa                                        <  1:439403/151‑1 (MQ=255)                                                                                                       gCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACaaaa                                         >  2:321282/1‑150 (MQ=255)                                                                                                             gCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCt                                  >  2:56812/1‑151 (MQ=255)                                                                                                              ccTGAGTTTTAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCtt                                 >  2:318550/1‑151 (MQ=255)                                                                                                                      ttAGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCt                          <  2:10248/150‑1 (MQ=255)                                                                                                                        aGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTc                       >  2:389743/1‑151 (MQ=255)                                                                                                                        aGTGCTGATGAAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTc                       <  2:65814/151‑1 (MQ=255)                                                                                                                              gatgaAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCtttt                 <  2:701538/151‑1 (MQ=255)                                                                                                                                tgaAATAAAGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGt               >  1:782860/1‑151 (MQ=255)                                                                                                                                       aaGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGACACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTc        >  1:19128/1‑151 (MQ=255)                                                                                                                                        aGACGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTcc       >  2:269943/1‑151 (MQ=255)                                                                                                                                          aCGTATTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGt     <  2:168295/151‑1 (MQ=255)                                                                                                                                              aTTGATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTAt  <  1:361633/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TGGCCGCTAAAAATCATGACAACATGTCAGCTGTTCATATAGCCATTCATCATAAGCATTATAAAACTGTAGAAACATATTATGCTGCTATTAATGCAATCAGCCAAAGCCTGAGTTTTAGTGCTGATGAAATAAAGACGTATTTATAACCAGCTATTTACGATATATATCGATACATAATTTTATTTTTTTCATAAAACGCAGAAACAGTGACCATGTGCACAACGCCACCCGACGGCTGTAATAACAAAAAATCCCTTGATGCCTGTCCCTTTTGTTACACTCCGTTAT  >  NC_000913/1804845‑1805135

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