Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,783,249	T→C	pseudogene (56/240 nt)	pinH ←	pseudogene, invertase resolvase family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,783,249	0	T	C	100.0%	302.6 / NA	85	pseudogene (56/240 nt)	pinH	pseudogene, invertase resolvase family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (45/40);  total (45/40)

CCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGTTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTGTTATTGCTGG  >  NC_000913/2783118‑2783378                                                                                                                                    |                                                                                                                                  cccATTGCTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTcc                                                                                                                                >  1:758228/1‑135 (MQ=255) cccATTGCTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTcc                                                                                                                                <  2:758228/135‑1 (MQ=255)      tgttgGGCAAATATAGATTGGGCCAGCGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCgatg                                                                                                                       >  1:1295890/1‑139 (MQ=255)            gCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTcc                                                                                                                 >  1:1089265/1‑139 (MQ=255)                atatAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGa                                                                                                              >  2:1367323/1‑138 (MQ=255)                atatAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGa                                                                                                              <  1:1367323/138‑1 (MQ=255)                         gggCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACGAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGtttttt                                                                                                    <  2:1295890/139‑1 (MQ=255)                          ggCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTcc                                                                                                                                <  1:19998/110‑1 (MQ=255)                          ggCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTcc                                                                                                                                >  2:19998/1‑110 (MQ=255)                          ggCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTc                                                                                                   >  2:440875/1‑139 (MQ=255)                          ggCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTc                                                                                                   >  1:245837/1‑139 (MQ=255)                           gCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGtttttt                                                                                                    >  1:416586/1‑137 (MQ=255)                           gCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGAACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGtttttt                                                                                                    <  2:416586/137‑1 (MQ=255)                                agCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCgatgat                                                                                                                     <  2:177000/115‑1 (MQ=255)                                agCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCgatgat                                                                                                                     >  1:177000/1‑115 (MQ=255)                                    cGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGAcc                                                                                         >  1:118990/1‑139 (MQ=255)                                    cGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGAcc                                                                                         >  1:969762/1‑139 (MQ=255)                                    cGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTACCGATGATGTCCTGACCAGTTTTTTCAGTCCGGAcc                                                                                         >  1:1396353/1‑139 (MQ=255)                                          ttAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGc                                                                                   >  1:773742/1‑139 (MQ=255)                                            aaTACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCa                                                                                 <  2:969762/139‑1 (MQ=255)                                            aaTACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCa                                                                                 <  2:118990/139‑1 (MQ=255)                                            aaTACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCa                                                                                 <  2:1396353/139‑1 (MQ=255)                                            aaTACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCa                                                                                 >  2:845078/1‑139 (MQ=255)                                                 cACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGtt                                                                            <  2:773742/139‑1 (MQ=255)                                                 cACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGtt                                                                            >  2:958749/1‑139 (MQ=255)                                                   cGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGtttt                                                                          >  1:755445/1‑139 (MQ=255)                                                    gTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTa                                                                         >  1:1484858/1‑139 (MQ=255)                                                    gTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTa                                                                         <  2:245837/139‑1 (MQ=255)                                                       cGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTccc                                                                                                                               <  2:1035826/82‑1 (MQ=255)                                                       cGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTccc                                                                                                                               >  1:1035826/1‑82 (MQ=255)                                                       cGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGt                                                                      >  2:81202/1‑139 (MQ=255)                                                        gCACAGCTCCTCCAACAGCACGACAAGATGTCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGtt                                                                     >  1:167869/1‑139 (MQ=255)                                                        gCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGtt                                                                     <  2:1484858/139‑1 (MQ=255)                                                          acaGCTCCTCCAACAGCACGACAAGATGTCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGa                                                                   <  2:167869/139‑1 (MQ=255)                                                              ctcctcCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCg                                                                                             <  1:624974/109‑1 (MQ=255)                                                              ctcctcCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCg                                                                                             >  2:624974/1‑109 (MQ=255)                                                              ctcctcCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTg                                                               >  1:1452693/1‑139 (MQ=255)                                                                          cACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGc                                                   >  2:502072/1‑139 (MQ=255)                                                                           aCGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCa                                                  >  1:837480/1‑139 (MQ=255)                                                                            cGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCAc                                                 >  1:1124079/1‑139 (MQ=255)                                                                             gACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTc                                                                                               <  2:189280/92‑1 (MQ=255)                                                                             gACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTc                                                                                               >  1:189280/1‑92 (MQ=255)                                                                             gACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                >  1:485874/1‑139 (MQ=255)                                                                             gACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                >  2:1177638/1‑139 (MQ=255)                                                                             gACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                >  1:257027/1‑139 (MQ=255)                                                                               cAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCa                                              <  2:1452693/139‑1 (MQ=255)                                                                                   aTGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGttc                                          <  1:1177638/139‑1 (MQ=255)                                                                                            aCTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGAcctcc                                                                                      >  1:1480841/1‑86 (MQ=255)                                                                                            aCTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGAcctcc                                                                                      <  2:1480841/86‑1 (MQ=255)                                                                                            aCTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTtcagtcag                                 <  2:755445/139‑1 (MQ=255)                                                                                            aCTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTtcagtcag                                 <  2:1124079/139‑1 (MQ=255)                                                                                              tGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGcc                                                                                  <  2:978504/88‑1 (MQ=255)                                                                                              tGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGcc                                                                                  >  1:978504/1‑88 (MQ=255)                                                                                               gcgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCg                              <  2:485874/139‑1 (MQ=255)                                                                                               gcgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCg                              <  1:958749/139‑1 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTcc                                                                                                                                <  2:735939/40‑1 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTcc                                                                                                                                >  1:735939/1‑40 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGa                                                                   >  2:841796/1‑101 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGa                                                                   <  1:841796/101‑1 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                <  2:1284595/120‑1 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                >  1:1284595/1‑120 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCgt                             >  2:862394/1‑139 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCgt                             >  2:351096/1‑139 (MQ=255)                                                                                                cgcCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCgt                             >  1:991963/1‑139 (MQ=255)                                                                                                   cccAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACg                                                                       >  2:959285/1‑94 (MQ=255)                                                                                                   cccAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACg                                                                       <  1:959285/94‑1 (MQ=255)                                                                                                   cccAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                >  1:1127046/1‑117 (MQ=255)                                                                                                   cccAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACg                                                <  2:1127046/117‑1 (MQ=255)                                                                                                    ccAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCgtggtg                         <  2:837480/139‑1 (MQ=255)                                                                                                     cAGTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTAACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGc                        >  2:1346923/1‑139 (MQ=255)                                                                                                       gTCGGTTCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAATTGGTAACCCTGAGCACACAGTTCTTCAGTCAGGCGTGGTGCAc                      >  1:1086949/1‑139 (MQ=255)                                                                                                            ttCAGTTTACACACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTaa                 <  2:1089265/139‑1 (MQ=255)                                                                                                             tCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAAc                <  1:845078/139‑1 (MQ=255)                                                                                                             tCAGTTTACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTAACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAAc                <  1:1346923/139‑1 (MQ=255)                                                                                                                 tttACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCAGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTg            <  1:502072/139‑1 (MQ=255)                                                                                                                  ttACAGACTAGTAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACACAGTTCTTCAGTCAGGCGTGGTGCACCGCAACGCTGt           <  2:1086949/139‑1 (MQ=255)                                                                                                                  ttACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTGt           <  2:257027/139‑1 (MQ=255)                                                                                                                  ttACAGACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTGt           <  1:351096/139‑1 (MQ=255)                                                                                                                       gACTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTGTTATTg      <  2:991963/139‑1 (MQ=255)                                                                                                                        aCTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTGAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTGTTATTGc     >  1:1178690/1‑139 (MQ=255)                                                                                                                         cTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTc                                                                                               >  1:1009192/1‑48 (MQ=255)                                                                                                                         cTAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTc                                                                                               <  2:1009192/48‑1 (MQ=255)                                                                                                                          tAGCAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGT‑ACGCTGTTATTGCTgg  >  2:106686/1‑139 (MQ=255)                                                                                                                             cAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCa                                                                                 <  2:713916/58‑1 (MQ=255)                                                                                                                             cAGTGCGCTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCa                                                                                 >  1:713916/1‑58 (MQ=255)                                                                                                                                    |                                                                                                                                  CCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGCCACATACTGCGCCCCAGTCGGTTCAGTTTACAGACTAGCAGTGCGTTCCCTGCCGATGATGTCCTGACCAGTTTTTTCAGTCCGGACCTCCTGCCACGGTTTTACGTTGAAGTGGTAACCCTGAGCACGCAGTTCTTCAGTCAGGCGTGGTGCACCGTAACGCTGTTATTGCTGG  >  NC_000913/2783118‑2783378

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