Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,863	C→T	pseudogene (19/48 nt)	rph ←	ribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,863	0	C	T	100.0%	249.6 / NA	77	pseudogene (19/48 nt)	rph	ribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/39);  total (38/39)

AAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATG  >  NC_000913/3815745‑3815992                                                                                                                       |                                                                                                                                  aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTcccc                                                                                                               >  1:682364/1‑139 (MQ=255)     tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCgg                                                                                                           >  1:227082/1‑139 (MQ=255)        cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                                        >  2:1119629/1‑139 (MQ=255)          cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                                   >  2:1679799/1‑126 (MQ=255)          cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                                   <  1:1679799/126‑1 (MQ=255)          cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCag                                                                                                      <  1:1031653/139‑1 (MQ=255)          cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCag                                                                                                      >  1:1826893/1‑139 (MQ=255)          cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCag                                                                                                      <  2:682364/139‑1 (MQ=255)            gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCagag                                                                                                    <  2:227082/139‑1 (MQ=255)            gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCagag                                                                                                    <  2:1826893/139‑1 (MQ=255)                 ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                                                                               <  1:1729461/139‑1 (MQ=255)                       aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgatt                                                                                                                        <  1:1340831/108‑1 (MQ=255)                       aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgatt                                                                                                                        >  2:1340831/1‑108 (MQ=255)                       aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTAGCTAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCGGAGCCAACAAGATg                                                                                         >  2:1123713/1‑139 (MQ=255)                        tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                                                        >  1:965257/1‑139 (MQ=255)                        tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                                                        <  1:1119629/139‑1 (MQ=255)                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     <  1:1492192/139‑3 (MQ=255)                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     <  2:519888/139‑3 (MQ=255)                             tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACGATGGATTCGATTCCCCTCGGGCCa                                                                                                       <  1:1243989/119‑1 (MQ=255)                             tCGCTCCTCATCTTACTTTTATACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCa                                                                                                       >  2:1243989/1‑119 (MQ=255)                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTc                                                                                                             >  1:722204/1‑107 (MQ=255)                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTc                                                                                                             <  2:722204/107‑1 (MQ=255)                                        cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTg                                                                       >  1:1182674/1‑139 (MQ=255)                                           aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCagag                                                                                                    >  2:1217799/1‑108 (MQ=255)                                           aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCagag                                                                                                    <  1:1217799/108‑1 (MQ=255)                                                tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACa                                                                                                                               >  2:618741/1‑76 (MQ=255)                                                tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACa                                                                                                                               <  1:618741/76‑1 (MQ=255)                                                tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTc                                                               >  2:368431/1‑139 (MQ=255)                                                tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTc                                                               >  2:595132/1‑139 (MQ=255)                                                 cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCg                                                              >  2:31574/1‑139 (MQ=255)                                                  tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGc                                                             >  2:1594460/1‑139 (MQ=255)                                                     aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCtt                                                          <  2:1182674/139‑1 (MQ=255)                                                      gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                         <  1:441686/92‑1 (MQ=255)                                                      gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                         >  2:441686/1‑92 (MQ=255)                                                       aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     <  1:529109/111‑3 (MQ=255)                                                       aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     >  2:529109/1‑109 (MQ=255)                                                        cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                       <  1:1594460/139‑1 (MQ=255)                                                             aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                  <  2:300516/139‑1 (MQ=255)                                                              aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcc                                                 >  2:1633505/1‑139 (MQ=255)                                                                 ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCtt                                                                               <  1:1602167/106‑1 (MQ=255)                                                                 ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCtt                                                                               >  2:1602167/1‑106 (MQ=255)                                                                  gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                             <  1:368431/139‑1 (MQ=255)                                                                  gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                             <  1:595132/139‑1 (MQ=255)                                                                     aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTg                                                                       >  1:780084/1‑110 (MQ=255)                                                                     aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTg                                                                       <  2:780084/110‑1 (MQ=255)                                                                                gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCagag                                                                                                    >  2:35189/1‑71 (MQ=255)                                                                                gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCagag                                                                                                    <  1:35189/71‑1 (MQ=255)                                                                                 ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                                                        <  2:1046394/82‑1 (MQ=255)                                                                                 ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                                                        >  1:1046394/1‑82 (MQ=255)                                                                                 ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCt                                                                >  2:1002727/1‑105 (MQ=255)                                                                                 ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCt                                                                <  1:1002727/105‑1 (MQ=255)                                                                                   ttAAAAATCAGTTTGTCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTCCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGcc                            >  1:1224607/1‑139 (MQ=255)                                                                                   ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                                        <  2:1535318/64‑1 (MQ=255)                                                                                   ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                                        >  1:1535318/1‑64 (MQ=255)                                                                                   ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                                                                               >  2:1521819/1‑73 (MQ=255)                                                                                   ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                                                                               <  1:1521819/73‑1 (MQ=255)                                                                                    tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                           <  2:1224607/139‑1 (MQ=255)                                                                                     aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCg                          >  1:1566829/1‑139 (MQ=255)                                                                                           cAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                    <  1:1633505/139‑1 (MQ=255)                                                                                            aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGc                                                                 >  2:141211/1‑93 (MQ=255)                                                                                            aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGc                                                                 <  1:141211/93‑1 (MQ=255)                                                                                            aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                       <  2:1573697/103‑1 (MQ=255)                                                                                            aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                       >  1:1573697/1‑103 (MQ=255)                                                                                              tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     >  1:934353/1‑70 (MQ=255)                                                                                              tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     <  2:934353/72‑3 (MQ=255)                                                                                               ttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                <  2:1566829/139‑1 (MQ=255)                                                                                                tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat               >  1:1294610/1‑139 (MQ=255)                                                                                                  ccAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca             <  1:31574/139‑1 (MQ=255)                                                                                                  ccAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca             <  2:965257/139‑1 (MQ=255)                                                                                                        ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     >  2:332827/1‑60 (MQ=255)                                                                                                        ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAgag                                                                                     <  1:332827/62‑3 (MQ=255)                                                                                                        ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt       >  2:316265/1‑139 (MQ=255)                                                                                                        ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt       >  2:931719/1‑139 (MQ=255)                                                                                                           ccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa    >  2:318686/1‑139 (MQ=255)                                                                                                             ttCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg  <  2:1294610/139‑1 (MQ=255)                                                                                                                 gCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt           <  1:23566/126‑1 (MQ=255)                                                                                                                 gCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAG‑AGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt           >  2:23566/1‑126 (MQ=255)                                                                                                                       |                                                                                                                                  AAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATG  >  NC_000913/3815745‑3815992

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