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%	212.9 / NA	66	pseudogene (50/1476 nt)	arpB	pseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/31);  total (35/31)

TAACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTT  >  NC_000913/1803432‑1803691                                                                                                                                        |                                                                                                                             tAACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATAAACCCAACAGTAAATCTCAATGGTAATGc                                                                                                                           >  2:1459474/1‑139 (MQ=255)   aCCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCAt                                                                                                                         >  2:257539/1‑139 (MQ=255)    ccATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATg                                                                                                                        >  1:865278/1‑139 (MQ=255)           cGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTg                                                                                                                 <  2:618754/139‑1 (MQ=255)              aGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCaa                                                                                                              >  2:1211971/1‑139 (MQ=255)               gCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAt                                                                                                             >  2:1225391/1‑139 (MQ=255)                   aCTGGATTATATTGAACAGGAACATTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCTACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTcc                                                                                                         <  2:705174/139‑1 (MQ=255)                        aTTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTg                                                                                                    >  2:938000/1‑139 (MQ=255)                          tataTTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGAc                                                                                                  <  1:1005560/139‑1 (MQ=255)                                 aacaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTAttt                                                                                           >  2:1408593/1‑139 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                      >  1:1406650/1‑139 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                      >  2:274721/1‑139 (MQ=255)                                      gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                      <  1:274721/139‑1 (MQ=255)                                           gTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                 <  1:1408593/139‑1 (MQ=255)                                           gTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                 <  1:257539/139‑1 (MQ=255)                                             tAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATc                                                                               <  1:1211971/139‑1 (MQ=255)                                                 tgctgAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTc                                                                           >  1:410757/1‑139 (MQ=255)                                                          tCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATa                                                                  <  1:1459474/139‑1 (MQ=255)                                                           cGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                 >  2:554480/1‑139 (MQ=255)                                                           cGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                 <  2:1266621/139‑1 (MQ=255)                                                             tAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                            <  1:1311110/126‑1 (MQ=255)                                                             tAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                            >  2:1311110/1‑126 (MQ=255)                                                             tAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGc                                                               >  1:214732/1‑139 (MQ=255)                                                                  aTTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTg                                                          <  2:214732/139‑1 (MQ=255)                                                                     gTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGAtt                                                       >  1:863280/1‑139 (MQ=255)                                                                      tATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTc                                                      >  2:353649/1‑139 (MQ=255)                                                                        tGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                 <  1:202997/94‑1 (MQ=255)                                                                        tGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                 >  2:202997/1‑94 (MQ=255)                                                                         gAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCAtt                                                   <  2:865278/139‑1 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                 >  2:167118/1‑88 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                 <  1:167118/88‑1 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                         <  1:455398/112‑1 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                         >  2:455398/1‑112 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                              >  1:1399955/1‑139 (MQ=255)                                                                              tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                              <  2:1406650/139‑1 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                 >  1:93637/1‑118 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                 <  2:93637/118‑1 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                            >  1:1062700/1‑139 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                            >  1:469067/1‑139 (MQ=255)                                                                                tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                            >  1:618775/1‑139 (MQ=255)                                                                                   aaaaaaaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAGGGa                                         <  1:938000/139‑1 (MQ=255)                                                                                          tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCa                                  <  1:554480/139‑1 (MQ=255)                                                                                          tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCa                                  >  2:109375/1‑139 (MQ=255)                                                                                                 agaAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCaa                           <  2:1399955/139‑1 (MQ=255)                                                                                                   aaaaaaaCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAcc                         <  2:1592905/139‑1 (MQ=255)                                                                                                        aaCTATACACCCAACTGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGACGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGAtt                    <  2:426545/139‑1 (MQ=255)                                                                                                        aaCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGAtt                    <  2:618775/139‑1 (MQ=255)                                                                                                           acagcacCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCGTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAATGTCAACCTGATTaaa                 <  2:1104807/135‑1 (MQ=255)                                                                                                             taCACCCATCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACa               >  2:444446/1‑139 (MQ=255)                                                                                                             taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCg                                                                                    >  1:1498507/1‑70 (MQ=255)                                                                                                             taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCg                                                                                    <  2:1498507/70‑1 (MQ=255)                                                                                                             taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                            >  1:663668/1‑78 (MQ=255)                                                                                                             taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                            <  2:663668/78‑1 (MQ=255)                                                                                                                    aaCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTCTTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcc        >  1:892376/1‑139 (MQ=255)                                                                                                                     aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGa                                                         >  2:338187/1‑89 (MQ=255)                                                                                                                     aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGa                                                         <  1:338187/89‑1 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATc                                                                               >  1:898947/1‑63 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATc                                                                               <  2:898947/63‑1 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcct       <  2:1117736/135‑1 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt   <  1:353649/139‑1 (MQ=255)                                                                                                                         tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTGTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcct       >  1:1117736/1‑135 (MQ=255)                                                                                                                          aaaTCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGtt  >  1:399659/1‑139 (MQ=255)                                                                                                                          aaaTCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGtt  >  2:19124/1‑139 (MQ=255)                                                                                                                          aaaTCTCAATGGTAATGCAGGTTTTTTGCTATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACtcctcctgt  >  1:1104807/1‑136 (MQ=255)                                                                                                                                   tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                              <  2:659687/86‑1 (MQ=255)                                                                                                                                   tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                              >  1:659687/1‑86 (MQ=255)                                                                                                                                        |                                                                                                                             TAACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTT  >  NC_000913/1803432‑1803691

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