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

ACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCA  >  NC_000913/1803434‑1803698                                                                                                                                      |                                                                                                                                    aCCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCAt                                                                                                                                >  1:519964/1‑139 (MQ=255)  ccATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATg                                                                                                                               >  1:1173637/1‑139 (MQ=255)      ttGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTCATGCATGtttt                                                                                                                           >  2:1048953/1‑139 (MQ=255)         cGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTg                                                                                                                        <  1:180857/139‑1 (MQ=255)               aTACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAtct                                                                                                                  <  1:919193/139‑1 (MQ=255)                         atatTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGAATGAAACTGTCTAAAAAAATCCATAGAGAAAGAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                        <  1:1048953/139‑1 (MQ=255)                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCa                                                                                                                      <  2:740343/119‑1 (MQ=255)                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCa                                                                                                                      >  1:740343/1‑119 (MQ=255)                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCT‑AAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAAAGGTAATGCATGTTTTTTGCAAtct                                                                                                                  >  1:966499/1‑122 (MQ=255)                               aacaggaacagCTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCT‑AAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAtct                                                                                                                  <  2:966499/122‑1 (MQ=255)                                    gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             >  1:1531192/1‑139 (MQ=255)                                    gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             >  1:813701/1‑139 (MQ=255)                                    gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             >  1:190367/1‑139 (MQ=255)                                         gTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                        >  1:339331/1‑139 (MQ=255)                                           tAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATc                                                                                      <  2:339331/139‑1 (MQ=255)                                                 ctgAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                                <  2:1173637/139‑1 (MQ=255)                                                    aaGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                        >  1:280811/1‑112 (MQ=255)                                                    aaGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACg                                                                                                        <  2:280811/112‑1 (MQ=255)                                                    aaGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGc                                                                             >  1:995423/1‑139 (MQ=255)                                                                  tGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGAt                                                               >  1:741233/1‑139 (MQ=255)                                                                            tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                                     <  1:464321/139‑1 (MQ=255)                                                                              tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTTCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                                   <  2:1519123/107‑1 (MQ=255)                                                                              tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTTCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                                   >  1:1519123/1‑107 (MQ=255)                                                                              tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                                   >  2:1043565/1‑139 (MQ=255)                                                                                  aaaaaaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGaa                                               <  2:519964/139‑1 (MQ=255)                                                                                      aaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCa                                                                                 >  2:410910/1‑101 (MQ=255)                                                                                      aaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATCTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCa                                                                                 <  1:410910/101‑1 (MQ=255)                                                                                        tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAAt                                                                          >  1:871818/1‑106 (MQ=255)                                                                                        tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAAt                                                                          <  2:871818/106‑1 (MQ=255)                                                                                        tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCa                                         >  1:1356817/1‑139 (MQ=255)                                                                                             agagAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTc                                    <  2:1245736/139‑1 (MQ=255)                                                                                                gaAAAATACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCacc                                 <  1:122525/139‑3 (MQ=255)                                                                                                gaAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAc                                 <  2:813701/139‑1 (MQ=255)                                                                                                gaAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAc                                 <  2:1531192/139‑1 (MQ=255)                                                                                                    aaaaCTATACACCCAACAGTGAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGGAAGGTCAACCTGa                             <  1:1080897/139‑1 (MQ=255)                                                                                                         tataCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTaaa                        >  1:1144096/1‑139 (MQ=255)                                                                                                           taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                        >  1:1595216/1‑73 (MQ=255)                                                                                                           taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                        <  2:1595216/73‑1 (MQ=255)                                                                                                           taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACa                      >  1:543120/1‑139 (MQ=255)                                                                                                              acCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTa                   <  2:741233/139‑1 (MQ=255)                                                                                                                  aaCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcc               <  2:543120/139‑1 (MQ=255)                                                                                                                   aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcct              >  2:1065036/1‑139 (MQ=255)                                                                                                                       tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt          >  2:1385272/1‑139 (MQ=255)                                                                                                                       tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt          <  2:1144096/139‑1 (MQ=255)                                                                                                                          aTCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                                >  2:453208/1‑66 (MQ=255)                                                                                                                          aTCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                                <  1:453208/66‑1 (MQ=255)                                                                                                                            ctcAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGAtt                           >  2:235485/1‑117 (MQ=255)                                                                                                                            ctcAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGAtt                           <  1:235485/117‑1 (MQ=255)                                                                                                                              cAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAAt                                                                          <  2:1408109/68‑1 (MQ=255)                                                                                                                              cAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAAt                                                                          >  1:1408109/1‑68 (MQ=255)                                                                                                                               aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCGTTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCa  <  2:190367/139‑1 (MQ=255)                                                                                                                               aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCa  >  1:1138039/1‑139 (MQ=255)                                                                                                                                 tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGa       >  2:1390355/1‑132 (MQ=255)                                                                                                                                 tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGa       <  1:1390355/132‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    ACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCA  >  NC_000913/1803434‑1803698

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