Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,428,079	C→T	pseudogene (386/402 nt) *213* (TAG→TAA)	ytfA → ytfB ←	pseudogene, related to transcriptional regulators OapA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,428,079	0	C	T	100.0%	249.7 / NA	77	pseudogene (386/402 nt) *213* (TAG→TAA)	ytfA ytfB	pseudogene, related to transcriptional regulators OapA family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/39);  total (38/39)

TATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTT  >  NC_000913/4427949‑4428209                                                                                                                                   |                                                                                                                                   tattatTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  2:999280/139‑1 (MQ=255) tattatTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  2:913555/139‑1 (MQ=255) tattatTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  1:982004/139‑1 (MQ=255) tattatTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  1:514039/139‑1 (MQ=255)  attattTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCa                                                                                                                           <  2:541067/139‑1 (MQ=255)   ttattTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAACCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            >  2:429297/1‑137 (MQ=255)   ttattTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAACCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  1:429297/137‑1 (MQ=255)          aTCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            >  2:855871/1‑130 (MQ=255)          aTCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  1:855871/130‑1 (MQ=255)                     tACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCa                                                                                                                           >  1:1436181/1‑120 (MQ=255)                     tACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCa                                                                                                                           <  2:1436181/120‑1 (MQ=255)                        gCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGTCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                                     <  1:1797822/139‑1 (MQ=255)                         cTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                                    >  2:1061033/1‑139 (MQ=255)                         cTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                                    >  2:165321/1‑139 (MQ=255)                          ttGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                                   >  2:1295686/1‑139 (MQ=255)                          ttGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                                   >  2:305737/1‑139 (MQ=255)                            ggCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCg                                                                                                 >  2:335138/1‑139 (MQ=255)                             gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACcgcg                                                                                                                             >  1:1877938/1‑110 (MQ=255)                             gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACcgcg                                                                                                                             <  2:1877938/110‑1 (MQ=255)                             gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                                          <  1:1335339/129‑1 (MQ=255)                             gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                >  2:2014289/1‑139 (MQ=255)                             gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                <  1:305737/139‑1 (MQ=255)                             gCCTGGTTAATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCagc                                                                                                          >  2:1335339/1‑127 (MQ=255)                              ccTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCggg                                                                                               <  1:2014289/139‑1 (MQ=255)                                                tATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGat                                                                             >  2:1836380/1‑139 (MQ=255)                                                     tttCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTTACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGatta                                                                           <  2:1280850/136‑1 (MQ=255)                                                     tttCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                                   >  2:2014114/1‑112 (MQ=255)                                                     tttCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                                   <  1:2014114/112‑1 (MQ=255)                                                     tttCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGatta                                                                           >  1:1280850/1‑136 (MQ=255)                                                       tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            <  1:243027/85‑1 (MQ=255)                                                       tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                            >  2:243027/1‑85 (MQ=255)                                                        cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCt                                                                                   <  2:553417/125‑1 (MQ=255)                                                        cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCt                                                                                   >  1:553417/1‑125 (MQ=255)                                                        cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTcccc                                                                     >  2:856468/1‑139 (MQ=255)                                                        cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTccc                                                                      >  2:961014/1‑138 (MQ=255)                                                        cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTccc                                                                      <  1:961014/138‑1 (MQ=255)                                                         aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCg                                                                    >  1:1930408/1‑139 (MQ=255)                                                          ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGt                                                                   >  1:1675163/1‑139 (MQ=255)                                                          ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGt                                                                   >  1:233326/1‑139 (MQ=255)                                                              aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAAt                                                                                                                      >  1:2000229/1‑84 (MQ=255)                                                              aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAAt                                                                                                                      <  2:2000229/84‑1 (MQ=255)                                                              aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                               <  2:233326/139‑1 (MQ=255)                                                              aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                               <  1:1061033/139‑1 (MQ=255)                                                              aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                               <  2:1930408/139‑1 (MQ=255)                                                               ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                                    <  1:1564579/101‑1 (MQ=255)                                                               ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                                    >  2:1564579/1‑101 (MQ=255)                                                                       aaGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAg                                                      <  2:1675163/139‑1 (MQ=255)                                                                            aGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGAGTATTCCCCGTATCAATGGTTAAGCCCGt                                                 <  1:165321/139‑1 (MQ=255)                                                                              aCGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTca                                               >  2:1286837/1‑139 (MQ=255)                                                                               cGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTcac                                              <  1:856468/139‑1 (MQ=255)                                                                                   aGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                >  2:1719094/1‑85 (MQ=255)                                                                                   aGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                <  1:1719094/85‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                >  1:2040713/1‑80 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                <  2:2040713/80‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCt                                     >  2:1994527/1‑139 (MQ=255)                                                                                                 aCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTcc                                                                       >  2:751486/1‑96 (MQ=255)                                                                                                 aCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTcc                                                                       <  1:751486/96‑1 (MQ=255)                                                                                                  cAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGtt                                                         >  1:534291/1‑109 (MQ=255)                                                                                                  cAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGtt                                                         <  2:534291/109‑1 (MQ=255)                                                                                                          aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                               <  2:822582/95‑1 (MQ=255)                                                                                                          aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                               >  1:822582/1‑95 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                        >  1:1428342/1‑65 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                        <  2:1428342/65‑1 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGattatt                                                                         <  1:984297/80‑1 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGattatt                                                                         >  2:984297/1‑80 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCa              >  1:120703/1‑139 (MQ=255)                                                                                                                 ttttGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                        <  1:1667620/63‑1 (MQ=255)                                                                                                                 ttttGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                        >  2:1667620/1‑63 (MQ=255)                                                                                                                     ggTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGatta                                                                           >  2:1598350/1‑72 (MQ=255)                                                                                                                     ggTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGatta                                                                           <  1:1598350/72‑1 (MQ=255)                                                                                                                      gTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                        >  1:1569542/1‑58 (MQ=255)                                                                                                                      gTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                        <  2:1569542/58‑1 (MQ=255)                                                                                                                         gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGtt    <  1:1994527/139‑1 (MQ=255)                                                                                                                          aCCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGttt   >  1:359312/1‑139 (MQ=255)                                                                                                                           cccTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGtttt  <  1:1836380/139‑1 (MQ=255)                                                                                                                               tCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                                     <  1:1458916/52‑1 (MQ=255)                                                                                                                               tCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                                     >  2:1458916/1‑52 (MQ=255)                                                                                                                                   |                                                                                                                                   TATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTT  >  NC_000913/4427949‑4428209

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