BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,617	C→G	8.5%	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	8.5%	340.5 / 19.7	106	noncoding (35/76 nt)	glyY	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (52/45); new base G (5/4); total (57/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 108 total aligned reads displayed.

AAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTA  >  NC_000913/4392470‑4392765

aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGcc                                                                                                                                                     >  1:745825/1‑149 (MQ=255)

  gCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAATCAGCTGAAAGGAACAGACCCCCCAAGCGGAAATAGGTGAGTTGCTAGAGCACGACCTTGCCa                                                                                                                                                    >  4:458368/1‑149 (MQ=1)

     ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                                                                >  1:432000/1‑149 (MQ=39)

     ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                                                                <  3:828886/149‑1 (MQ=38)

      gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTc                                                                                                                                               >  1:581430/1‑149 (MQ=39)

        aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                             <  1:164293/149‑2 (MQ=255)

        aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                             <  1:313545/149‑1 (MQ=38)

             tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                                                                                                                                        >  3:685076/1‑149 (MQ=38)

              cAGTTGGTAGAGCACGACCTTGCGAAGGTCGGGGCCGCG‑AGTTCGAGTCCCGTTTCCCGCTCGAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                       <  2:376361/149‑1 (MQ=11)

                gtgggTAGGGCACGACCGGGCCAGGGTCGGGGTCGCG‑AGTCCGAGTCTCGTCTCCGGCGACAGACTTGAAAA‑GTGCGGGAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                     <  2:124765/146‑1 (MQ=11)

                gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                     <  3:21970/149‑1 (MQ=38)

                   tgTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                                  <  2:14507/148‑1 (MQ=38)

                       gagCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGTTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                              >  4:447580/1‑149 (MQ=38)

                       gagCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                              <  3:310387/149‑1 (MQ=38)

                        agCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                             >  3:321840/1‑149 (MQ=39)

                         gCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCATGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                                                                                            >  4:559307/1‑149 (MQ=38)

                          cccGACCTTGGCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                                                                                           >  1:282955/3‑149 (MQ=255)

                           actgctgtGCCAAGGTC‑GTGTCGGGTATTTCGAGTCCCGTTACCTGCCCCAAAATTTGAAA‑GTGCTGCAAGGCAAAGACCCCCTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                          <  3:232754/142‑1 (MQ=11)

                           aCGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                          <  1:521555/149‑1 (MQ=38)

                                    ccAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGc                                                                                                                 >  1:866258/1‑149 (MQ=40)

                                     cAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                                >  1:191002/1‑149 (MQ=39)

                                      aaGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                               >  1:123350/1‑149 (MQ=38)

                                          tCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                                                                                           <  3:396052/149‑1 (MQ=38)

                                            ggggTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAtt                                                                                                         >  2:555869/1‑149 (MQ=40)

                                                  gcg‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctc                                                                                                   <  1:304387/149‑1 (MQ=255)

                                                    g‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctctc                                                                                                 >  3:635565/1‑149 (MQ=255)

                                                       gtgcGAGTCTCGTTTCCCGGTCGAAAATTGAAAA‑GTGCGGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCaa                                                                                               <  2:395914/146‑1 (MQ=14)

                                                       gTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCaa                                                                                               <  3:778618/149‑1 (MQ=255)

                                                       gTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCaa                                                                                               <  3:458368/149‑1 (MQ=255)

                                                           gAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGGCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaa                                                                                           >  2:555890/1‑149 (MQ=255)

                                                            aGTCTCGTTTCCCGCTCCAAATTTTGAAA‑GTGCTGTAAGGCACAGACCAGCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaa                                                                                          <  4:513794/149‑1 (MQ=255)

                                                            aGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaa                                                                                          <  4:310691/149‑1 (MQ=255)

                                                              tctcGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                                                        >  4:358345/1‑149 (MQ=255)

                                                              tctcGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                                                        >  1:575662/1‑149 (MQ=255)

                                                              tcgcGTATCCCGCCCGAAAATTTGAAA‑GTGCTGTAGGGCACAGACGACCTAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                                                        <  2:75609/146‑1 (MQ=17)

                                                                tcGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                     >  3:767938/1‑134 (MQ=255)

                                                                tcGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctc                                                                                                   <  2:843734/136‑1 (MQ=255)

                                                                tcGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctc                                                                                                   >  1:843734/1‑136 (MQ=255)

                                                                tcGTTTCCCGCTCCAAAATTCGAAA‑TTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                     <  4:767938/134‑1 (MQ=255)

                                                                    ttCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                                                                  <  3:878906/149‑1 (MQ=255)

                                                                     tCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAg                                                                                 <  4:24811/149‑1 (MQ=255)

                                                                        cGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGa                                                                              >  1:24490/1‑149 (MQ=255)

                                                                             cAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                                                         >  1:336557/1‑149 (MQ=255)

                                                                             cAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                                                         >  2:5021/1‑149 (MQ=255)

                                                                             cAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                                                         >  1:585024/1‑149 (MQ=255)

                                                                               aaaTTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATg                                                                       >  1:419568/1‑149 (MQ=255)

                                                                                     gAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAt                                                                 >  2:255066/1‑149 (MQ=255)

                                                                                       aa‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                                                               >  1:233241/1‑149 (MQ=255)

                                                                                            gCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                                                           <  1:150128/149‑1 (MQ=255)

                                                                                              tGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                                                         >  3:691884/1‑149 (MQ=255)

                                                                                              tGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGATCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                                                         >  3:768643/1‑149 (MQ=255)

                                                                                               gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttt                                                        <  1:521290/149‑1 (MQ=255)

                                                                                                 aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCTAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                      <  2:222363/149‑1 (MQ=255)

                                                                                                 aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                      >  4:66799/1‑149 (MQ=255)

                                                                                                  aGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                              >  4:374122/1‑76 (MQ=32)

                                                                                                  aGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                              <  3:374122/76‑1 (MQ=32)

                                                                                                   ggCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttgt                                                    <  2:585024/149‑1 (MQ=255)

                                                                                                      acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt                                                 >  4:323773/1‑149 (MQ=255)

                                                                                                      acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt                                                 >  3:421868/1‑149 (MQ=255)

                                                                                                      acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATGCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt                                                 <  1:351970/149‑1 (MQ=255)

                                                                                                       caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAACGACGCGATGCGTTATTGCTGGTTTTTGTTGtc                                                <  3:559307/149‑1 (MQ=255)

                                                                                                         taccacgCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc                                              <  2:756084/148‑1 (MQ=255)

                                                                                                         gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc                                              <  3:5501/149‑1 (MQ=255)

                                                                                                          accaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct                                               <  3:125186/147‑1 (MQ=255)

                                                                                                          accaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct                                               >  4:125186/1‑147 (MQ=255)

                                                                                                              cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACa                                         >  3:882310/1‑149 (MQ=255)

                                                                                                              cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACa                                         >  4:332114/1‑149 (MQ=255)

                                                                                                               ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa                                        >  4:50846/1‑149 (MQ=255)

                                                                                                               ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa                                        >  3:65593/1‑149 (MQ=255)

                                                                                                                cacgCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaaa                                       <  4:754219/146‑1 (MQ=255)

                                                                                                                   gCGGGAATAGCTCAGTTGGTAGAGCACGATCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc                                              <  1:655006/139‑1 (MQ=255)

                                                                                                                   gCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc                                              >  2:655006/1‑139 (MQ=255)

                                                                                                                   gCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAActc                                    >  3:668172/1‑149 (MQ=255)

                                                                                                                   gCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAActc                                    >  1:374861/1‑149 (MQ=255)

                                                                                                                     gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCtt                                  <  4:635565/149‑1 (MQ=255)

                                                                                                                             cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACag                          <  2:822109/149‑1 (MQ=255)

                                                                                                                             cTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGCGACGGAAGTTCGATTCTCGTGTCTCGCTCCAAATTCTACTGTCAATAAAAGAGCCACAGCGACGCGAGGCGTTATTGCGGGGTTTTGTTTTCTCTGACAAACTCTTCTAAACag                          >  2:16316/1‑149 (MQ=21)

                                                                                                                              tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                         >  1:701948/1‑149 (MQ=255)

                                                                                                                              tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGACTCTCGTTTCCCGCTCCAAATTCTTCTCGCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                         <  3:746718/149‑1 (MQ=255)

                                                                                                                               cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagag                        <  2:233241/149‑1 (MQ=255)

                                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTCGTAAACAGAGtt                      <  3:66799/149‑1 (MQ=255)

                                                                                                                                   tGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTAt                    <  4:668172/149‑1 (MQ=255)

                                                                                                                                    ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATc                   <  3:299505/149‑1 (MQ=255)

                                                                                                                                     gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATcc                  >  2:212950/1‑149 (MQ=255)

                                                                                                                                       agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagag                        >  2:49782/1‑141 (MQ=255)

                                                                                                                                       agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagag                        <  1:49782/141‑1 (MQ=255)

                                                                                                                                        gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCcaca               <  4:882310/149‑1 (MQ=255)

                                                                                                                                          gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCtt                                  >  3:435427/1‑128 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                                                         >  1:337222/3‑104 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagag                        >  2:601848/3‑137 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagag                        <  1:601848/135‑1 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                         >  4:423027/3‑136 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                         <  3:423027/134‑1 (MQ=255)

                                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGcc            >  3:474108/3‑149 (MQ=255)

                                                                                                                                           cccGACATTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                                                         <  2:337222/102‑1 (MQ=255)

                                                                                                                                              gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCt           <  4:854866/147‑1 (MQ=255)

                                                                                                                                              gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCt           >  3:854866/1‑147 (MQ=255)

                                                                                                                                               aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATCTGCTCTCACTGTAATATCCACAGCGCCCCGACGCGTTTTTGCTGGTTTTTGTTGCCTCTGACACACTCTTGTAAACAGAGTTATTCACAGCCTCAg        >  3:7498/1‑149 (MQ=25)

                                                                                                                                                  ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCt     >  3:248756/1‑149 (MQ=255)

                                                                                                                                                     ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTa  <  2:370156/149‑1 (MQ=255)

AAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCG‑AGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAA‑GTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGTA  >  NC_000913/4392470‑4392765

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