Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,617	C→G	5.2%	noncoding (35/76 nt)	glyY →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,617	0	C	G	5.2%	259.8 / 6.4	77	noncoding (35/76 nt)	glyY	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base G (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01

AAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCT  >  NC_000913/4392470‑4392762                                                                                                                                                     |                                                                                                                                                  aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGT‑AGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCa                                                                                                                                                 >  1:579362/1‑149 (MQ=255)         aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                          <  3:664467/149‑1 (MQ=40)         aTAGCTCAGTTGGTAGACCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                          <  2:183248/149‑1 (MQ=38)          tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                         <  1:851234/149‑1 (MQ=39)          tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCAAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                         <  4:674267/149‑1 (MQ=38)              tCAGTTGGTAGAGTACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                         <  3:19742/145‑1 (MQ=38)              tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTAATTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                         >  4:19742/1‑145 (MQ=255)                aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                          >  3:721258/1‑142 (MQ=39)                aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                          <  4:721258/142‑1 (MQ=38)                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTACTGTAAGGCACAGACCACCCAAGTGGGAATAGCTCAGTTGGTAGAGTACGACATTGACAAGGGCGAGGTtgctt                                                                                                                                  >  4:829742/1‑147 (MQ=11)                     gTAGAGCACGACCTAGCCAAGGACGGGGTCGCGAGTTCGAGTCCCGTTTCC‑CGCTCCAAATTTTAAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                                                              <  4:392941/149‑1 (MQ=11)                      tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                             >  3:102451/1‑149 (MQ=39)                           cccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                                                                                        >  3:642968/3‑149 (MQ=255)                            aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                       <  1:763518/149‑1 (MQ=38)                               aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                    >  4:197617/1‑149 (MQ=38)                               aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                    >  3:209585/1‑149 (MQ=39)                                ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                                                                                   >  3:318220/1‑149 (MQ=38)                                 cTTGCCAAGGTCGGGGTGGAGAGTGCGAGTCCCGTTTCCACGC‑CCAAATTTTAAAAGTGCTGTAACGAACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                  <  2:700116/149‑1 (MQ=11)                                 cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                  >  4:801467/1‑149 (MQ=38)                                      cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                             >  4:217557/1‑149 (MQ=38)                                           tCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGTCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                                                                                        >  4:95359/1‑149 (MQ=38)                                                gTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAACTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttctt                                                                                                   <  4:642968/149‑1 (MQ=255)                                                 tggagAGTTCGAGTCCCGTTTCC‑CGCTACAAAAGTGGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                  <  2:112602/145‑1 (MQ=14)                                                  cgcgAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATtcttct                                                                                                 <  1:524392/149‑1 (MQ=255)                                                   gcgAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctc                                                                                                >  3:94650/1‑149 (MQ=255)                                                      aGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                 >  3:202084/1‑129 (MQ=255)                                                      aGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                 <  4:202084/129‑1 (MQ=40)                                                           gAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaa                                                                                        >  1:141237/1‑149 (MQ=255)                                                            aGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaa                                                                                       <  4:94650/149‑1 (MQ=255)                                                             gTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAat                                                                                      >  2:61387/1‑149 (MQ=255)                                                              tctcGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                                                     >  4:503864/1‑149 (MQ=255)                                                              tctcGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                                                     >  3:218328/1‑149 (MQ=255)                                                                tcGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATc                                                                                   <  1:850278/149‑1 (MQ=255)                                                                    ttCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                                                               <  1:431339/149‑1 (MQ=255)                                                                      cc‑cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGc                                                                             >  2:840436/1‑149 (MQ=255)                                                                      cc‑cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGc                                                                             >  1:618311/1‑149 (MQ=255)                                                                       c‑cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCg                                                                            >  4:168925/1‑149 (MQ=255)                                                                            tCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcgc                                                                        >  2:530041/1‑149 (MQ=255)                                                                                      gAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAt                                                              >  1:775790/1‑149 (MQ=255)                                                                                              tGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                                                      >  1:201576/1‑149 (MQ=255)                                                                                               gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttt                                                     >  2:310814/1‑149 (MQ=255)                                                                                                 aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                   <  1:115538/149‑1 (MQ=255)                                                                                                 aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                   <  3:339344/149‑1 (MQ=255)                                                                                                 aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                   <  3:822446/149‑1 (MQ=255)                                                                                                      acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt                                              >  2:441418/1‑149 (MQ=255)                                                                                                       caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtc                                             <  3:113470/149‑1 (MQ=255)                                                                                                       caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTCCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtc                                             <  1:544006/149‑1 (MQ=255)                                                                                                            caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCCATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa                                        <  2:618311/149‑1 (MQ=255)                                                                                                            caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa                                        <  2:141237/149‑1 (MQ=255)                                                                                                            caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa                                        >  4:320776/1‑149 (MQ=255)                                                                                                               ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATcc                                                                                  <  3:126566/104‑1 (MQ=255)                                                                                                               ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATcc                                                                                  >  4:126566/1‑104 (MQ=255)                                                                                                               ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa                                     >  1:329015/1‑149 (MQ=255)                                                                                                                cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGCTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCCCCAAATACTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaaa                                    <  4:561939/149‑1 (MQ=255)                                                                                                                 aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc                                   <  2:775790/149‑1 (MQ=255)                                                                                                                  aGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt                                              >  3:366981/1‑137 (MQ=255)                                                                                                                  aGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt                                              <  4:366981/137‑1 (MQ=255)                                                                                                                     gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCtt                               <  4:102451/149‑1 (MQ=255)                                                                                                                      ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTg                              <  1:840436/149‑1 (MQ=255)                                                                                                                       gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGt                             <  1:208510/149‑1 (MQ=255)                                                                                                                       gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGt                             <  3:197617/149‑1 (MQ=255)                                                                                                                              tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagg                      >  4:60727/1‑148 (MQ=255)                                                                                                                              tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaaa                      <  1:310814/149‑3 (MQ=255)                                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAActct                                <  3:442834/136‑1 (MQ=255)                                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAActct                                >  4:442834/1‑136 (MQ=255)                                                                                                                                       agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCcac             >  2:645238/1‑149 (MQ=255)                                                                                                                                       agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAAGATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCcac             >  4:566886/1‑149 (MQ=255)                                                                                                                                         agCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAg           >  1:23112/1‑149 (MQ=255)                                                                                                                                          gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                      >  1:333089/1‑137 (MQ=255)                                                                                                                                          gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                      <  2:333089/137‑1 (MQ=255)                                                                                                                                          gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGc          >  2:814822/1‑149 (MQ=255)                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTATTGTCTCTGACAAACTCtt                               >  1:409042/3‑127 (MQ=255)                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTATTGTCTCTGACAAACTCtt                               <  2:409042/125‑1 (MQ=255)                                                                                                                                           cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGcc         >  2:738949/1‑149 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                                                        >  1:176024/2‑101 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                                                        <  2:176024/100‑1 (MQ=255)                                                                                                                                                  ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCt  >  2:199256/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC‑CGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCT  >  NC_000913/4392470‑4392762

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