Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,803,711	T→C	pseudogene (194/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,711	0	T	C	100.0%	239.0 / NA	67	pseudogene (194/1476 nt)	arpB	pseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/34);  total (33/34)

TTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTC  >  NC_000913/1803575‑1803845                                                                                                                                         |                                                                                                                                       tttttGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGccc                                                                                                                                      >  2:790353/1‑139 (MQ=255)   tttGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCcccc                                                                                                                                    <  1:790353/139‑1 (MQ=255)   tttGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCcccc                                                                                                                                    <  2:497967/139‑1 (MQ=255)    ttGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCt                                                                                                                                   <  1:1235875/139‑1 (MQ=255)     tGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCtt                                                                                                                                  >  2:1108944/1‑139 (MQ=255)         aTCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTc                                                                                                                              >  2:975137/1‑139 (MQ=255)                  tGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     <  1:570097/139‑1 (MQ=255)                  tGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     <  2:1534473/139‑1 (MQ=255)                       cTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCACCCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGtt                                                                                                                <  1:639486/139‑1 (MQ=255)                        tATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     <  2:512313/133‑1 (MQ=255)                        tATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     >  1:512313/1‑133 (MQ=255)                                         tGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCa                                                                                              <  2:298246/139‑1 (MQ=255)                                                 cAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTg                                                                                      >  1:452783/1‑139 (MQ=255)                                                   aTACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGtt                                                                                                                <  2:1398954/111‑1 (MQ=255)                                                   aTACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGtt                                                                                                                >  1:1398954/1‑111 (MQ=255)                                                     aCGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAgg                                                                                  >  2:500618/1‑139 (MQ=255)                                                         ccTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGt                                                                              >  1:551291/1‑139 (MQ=255)                                                           ttGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTcg                                                                            >  2:826938/1‑139 (MQ=255)                                                            tGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     >  1:294481/1‑97 (MQ=255)                                                            tGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     <  2:294481/97‑1 (MQ=255)                                                                 cATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTAt                                                                      <  1:826938/139‑1 (MQ=255)                                                                   ttGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCg                                                                                                                            <  2:706575/83‑1 (MQ=255)                                                                   ttGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCg                                                                                                                            >  1:706575/1‑83 (MQ=255)                                                                   ttGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATcg                                                                    <  2:602773/139‑1 (MQ=255)                                                                             gAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGa                                                          >  1:1255760/1‑139 (MQ=255)                                                                                aaGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAg                                                                               >  1:1161751/1‑115 (MQ=255)                                                                                aaGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAg                                                                               <  2:1161751/115‑1 (MQ=255)                                                                                     aaaGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCt                                                  <  1:500618/139‑1 (MQ=255)                                                                                     aaaGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCt                                                  <  2:1255760/139‑1 (MQ=255)                                                                                      aaGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGTAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                 <  1:975137/139‑1 (MQ=255)                                                                                             ccTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTcc                                                                                                                             >  1:1535289/1‑56 (MQ=255)                                                                                             ccTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTcc                                                                                                                             <  2:1535289/56‑1 (MQ=255)                                                                                             ccTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTa                                          >  2:855668/1‑139 (MQ=255)                                                                                                   tAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCATTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCTCTTTAGATAtt                                    >  2:1093603/1‑139 (MQ=255)                                                                                                      aCATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTc                                 >  2:689793/1‑139 (MQ=255)                                                                                                       cATTACTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCt                                <  2:551291/139‑1 (MQ=255)                                                                                                           aCTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAc                                                                                                                        >  2:1263021/1‑47 (MQ=255)                                                                                                           aCTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAc                                                                                                                        <  1:1263021/47‑1 (MQ=255)                                                                                                           aCTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                 <  2:921414/118‑1 (MQ=255)                                                                                                           aCTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTg                            >  2:1306645/1‑139 (MQ=255)                                                                                                           aCTCCTCCGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTAAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                 >  1:921414/1‑118 (MQ=255)                                                                                                               ctccGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGt                                                               >  1:439640/1‑100 (MQ=255)                                                                                                               ctccGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGt                                                               <  2:439640/100‑1 (MQ=255)                                                                                                               ctccGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGata                                      >  2:727395/1‑125 (MQ=255)                                                                                                               ctccGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGata                                      <  1:727395/125‑1 (MQ=255)                                                                                                                  cGTTGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAatg                     >  2:67625/1‑139 (MQ=255)                                                                                                                    ttGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     >  1:732765/1‑41 (MQ=255)                                                                                                                    ttGAATCCATACAGCAGCATGCCCCCTTAGTCCGGGAcgcg                                                                                                                     <  2:732765/41‑1 (MQ=255)                                                                                                                      gAATCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATgg                 <  1:689793/139‑1 (MQ=255)                                                                                                                       aaTCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCAtcctc                                                                                         <  2:149625/66‑1 (MQ=255)                                                                                                                       aaTCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCAtcctc                                                                                         >  1:149625/1‑66 (MQ=255)                                                                                                                         tCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCtt                               <  2:1126101/122‑1 (MQ=255)                                                                                                                         tCCATACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCtt                               >  1:1126101/1‑122 (MQ=255)                                                                                                                            aTACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTc                                 <  2:168747/117‑1 (MQ=255)                                                                                                                            aTACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTc                                 >  1:168747/1‑117 (MQ=255)                                                                                                                             tACAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACaa          >  1:123359/1‑139 (MQ=255)                                                                                                                              aCAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                 >  1:876697/1‑99 (MQ=255)                                                                                                                              aCAGCAGCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCtt                                                 <  2:876697/99‑1 (MQ=255)                                                                                                                               cagcagCATGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATg        <  2:123359/139‑1 (MQ=255)                                                                                                                                 gcagcaTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCa               >  2:268600/1‑130 (MQ=255)                                                                                                                                 gcagcaTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCa               <  1:268600/130‑1 (MQ=255)                                                                                                                                  cagcaTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTa                                                              >  1:987025/1‑82 (MQ=255)                                                                                                                                  cagcaTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTa                                                              <  2:987025/82‑1 (MQ=255)                                                                                                                                     caTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATcgc                                                                   <  2:1288599/74‑1 (MQ=255)                                                                                                                                     caTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATcgc                                                                   >  1:1288599/1‑74 (MQ=255)                                                                                                                                     caTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaaa             <  1:475674/128‑1 (MQ=255)                                                                                                                                     caTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATgtcagtc  >  1:623844/1‑139 (MQ=255)                                                                                                                                     caTGCCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTCCGCTTTAGATATTTTCTTTTGTCGAATGATGGCaaa             >  2:475674/1‑128 (MQ=255)                                                                                                                                         |                                                                                                                                       TTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCATACAGCAGCATGTCCCCTTAGTCCGGGACGCGGAGTTCAGAGCATTACTTCGCCATCCTCCTGCAGGGAGTCGCGTTATCGCGAGTAAGGATTTTGGCTTCGCTTTAGATATTTTCTTTTGTCGAATGATGGCAAACAATGTCAGTC  >  NC_000913/1803575‑1803845

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: ‑