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

ACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCA  >  NC_000913/1803434‑1803698                                                                                                                                       |                                                                                                                                    aCCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCAt                                                                                                                                >  1:1400465/1‑139 (MQ=255)      ttGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGtttt                                                                                                                           >  2:846198/1‑139 (MQ=255)                tACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAtctc                                                                                                                 >  1:1874381/1‑139 (MQ=255)                   tGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCtt                                                                                                              >  2:1076395/1‑139 (MQ=255)                    ggATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTc                                                                                                             <  2:1191918/139‑1 (MQ=255)                     gATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCt                                                                                                            <  2:1621433/139‑1 (MQ=255)                        tataTTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGAc                                                                                                         >  1:202607/1‑139 (MQ=255)                        tataTTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGAc                                                                                                         <  2:244278/139‑1 (MQ=255)                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTAttt                                                                                                  >  2:1563219/1‑139 (MQ=255)                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTAttt                                                                                                  >  1:552261/1‑139 (MQ=255)                                    gaacagTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGt                                                                                                                              <  2:358793/106‑1 (MQ=255)                                    gaacagTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGt                                                                                                                              >  1:358793/1‑106 (MQ=255)                                    gaacagTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTTACGCTATTTTTTGt                                                                                             <  1:846198/139‑1 (MQ=255)                                         gTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                        >  1:89567/1‑139 (MQ=255)                                                  tgAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                         >  2:693674/1‑129 (MQ=255)                                                  tgAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                         <  1:693674/129‑1 (MQ=255)                                                    aaGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGc                                                                             >  2:400746/1‑139 (MQ=255)                                                     aaaTCGTAACC‑TTGTAATAAAATGGTCTAAAAAAATCAATGGGCAAAAAATCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCa                                                                            <  1:489080/137‑1 (MQ=255)                                                         cGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             <  2:770570/118‑1 (MQ=255)                                                         cGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             >  1:770570/1‑118 (MQ=255)                                                         cGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                        <  1:1563219/139‑1 (MQ=255)                                                             aCCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGccc                                                                    <  1:1579839/139‑1 (MQ=255)                                                                 ttGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAAt                                                                          <  2:80594/129‑1 (MQ=255)                                                                 ttGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAAt                                                                          >  1:80594/1‑129 (MQ=255)                                                                  tGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                                <  2:138973/122‑1 (MQ=255)                                                                  tGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCAt                                                                                >  1:138973/1‑122 (MQ=255)                                                                  tGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGAt                                                               <  2:1400465/139‑1 (MQ=255)                                                                   gT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                        >  2:854378/1‑129 (MQ=255)                                                                   gT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                        <  1:854378/129‑1 (MQ=255)                                                                   gT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGAtt                                                              >  1:1585148/1‑139 (MQ=255)                                                                             tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                                     >  1:407684/1‑139 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                                   <  1:561660/107‑1 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGt                                                                                   >  2:561660/1‑107 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGa                                                                >  2:1079610/1‑126 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGa                                                                <  1:1079610/126‑1 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                                   >  2:214616/1‑139 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                                   >  1:151314/1‑139 (MQ=255)                                                                               tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAAt                                                   >  1:1005970/1‑139 (MQ=255)                                                                                    aaaaaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGaaa                                              >  2:680617/1‑139 (MQ=255)                                                                                    aaaaaTCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGaaa                                              >  1:264461/1‑139 (MQ=255)                                                                                         tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             <  1:1551022/87‑1 (MQ=255)                                                                                         tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             >  2:1551022/1‑87 (MQ=255)                                                                                         tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCa                                         >  1:1198717/1‑139 (MQ=255)                                                                                         tCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCa                                         >  2:646172/1‑139 (MQ=255)                                                                                             tAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGt                                     >  2:332386/1‑139 (MQ=255)                                                                                                 gaAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAc                                 <  2:1005970/139‑1 (MQ=255)                                                                                                     aaaaCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGa                             <  2:1874381/139‑1 (MQ=255)                                                                                                       aaCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGAtt                           <  1:1076395/139‑1 (MQ=255)                                                                                                        aCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTa                          <  1:646172/139‑1 (MQ=255)                                                                                                          tataCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAt                                                                                                                    <  2:361897/47‑1 (MQ=255)                                                                                                          tataCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAt                                                                                                                    >  1:361897/1‑47 (MQ=255)                                                                                                          tataCACCCAACAGTAAATCTCAACGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             <  1:742272/70‑1 (MQ=255)                                                                                                          tataCACCCAACAGTAAATCTCAACGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                             >  2:742272/1‑70 (MQ=255)                                                                                                            taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCTATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACa                      >  2:1227806/1‑139 (MQ=255)                                                                                                            taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                         <  2:547520/72‑1 (MQ=255)                                                                                                            taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                         >  1:547520/1‑72 (MQ=255)                                                                                                            taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACa                      >  2:401398/1‑139 (MQ=255)                                                                                                            taCACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTAAACCTGATTAAACa                      >  1:1690444/1‑139 (MQ=255)                                                                                                                 ccAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                         >  2:927393/1‑67 (MQ=255)                                                                                                                 ccAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                         <  1:927393/67‑1 (MQ=255)                                                                                                                  cAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCa                                                                            <  2:684195/79‑1 (MQ=255)                                                                                                                  cAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCa                                                                            >  1:684195/1‑79 (MQ=255)                                                                                                                   aaCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcc               >  1:1407802/1‑139 (MQ=255)                                                                                                                    aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                        >  2:1769871/1‑65 (MQ=255)                                                                                                                    aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                        <  1:1769871/65‑1 (MQ=255)                                                                                                                    aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcct              >  2:1803368/1‑139 (MQ=255)                                                                                                                      aGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAc                                 >  2:358408/1‑118 (MQ=255)                                                                                                                      aGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAc                                 <  1:358408/118‑1 (MQ=255)                                                                                                                      aGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcc               >  2:1794584/1‑136 (MQ=255)                                                                                                                      aGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcc               <  1:1794584/136‑1 (MQ=255)                                                                                                                        tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt          >  1:652734/1‑139 (MQ=255)                                                                                                                        tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt          <  2:407684/139‑1 (MQ=255)                                                                                                                        tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt          <  1:332386/139‑1 (MQ=255)                                                                                                                        tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt          <  2:1585148/139‑1 (MQ=255)                                                                                                                           aTCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAcc                                >  1:164576/1‑114 (MQ=255)                                                                                                                           aTCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAAcc                                <  2:164576/114‑1 (MQ=255)                                                                                                                               cAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGa                                                >  2:1587749/1‑94 (MQ=255)                                                                                                                               cAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGa                                                <  1:1587749/94‑1 (MQ=255)                                                                                                                                aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGa                             <  2:556527/112‑1 (MQ=255)                                                                                                                                aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGa                             >  1:556527/1‑112 (MQ=255)                                                                                                                                aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCa  <  2:151314/139‑1 (MQ=255)                                                                                                                                aaTGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCa  >  1:1721544/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                    ACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGT‑ATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAATCCA  >  NC_000913/1803434‑1803698

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