BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,254	A→T	9.4%	intergenic (+32/‑115)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,254	0	A	T	9.4%	332.8 / 18.9	118	intergenic (+32/‑115)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (59/46); new base T (7/5); total (67/51)
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.00e+00

Only 100 of 119 total aligned reads displayed.

GAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781107‑781403

gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCACATGTTAAGCGTGAAg                                                                                                                                                        >  1:115831/1‑150 (MQ=255)

     ggTCGGGGGTGCGATCCCGTCATCACCCACCACTCGGGTCTTTAGCTCAGTTGGTAGAGGAGTTGACTTTGTATCAATTGGTCGAAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGATGGATaa                                                                                                                                                   <  1:1997874/150‑1 (MQ=14)

      gTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTTGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                                                  <  2:94754/150‑1 (MQ=255)

         ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTTAACATAAATATAAGTTGTTAAGCGTGAAGGATCACggt                                                                                                                                               >  1:823816/1‑148 (MQ=14)

         ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCAGTTTTAACATCATACACAGATGTTAAGCGGGAAGTAGAACGtt                                                                                                                                               >  1:2360664/1‑150 (MQ=14)

         ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGTGAAGGATAACggt                                                                                                                                               >  1:1302642/1‑148 (MQ=255)

           cggTGCTATCCCGTCATCACCCACCACTCGGGTGGTTTTCTCCGTTGGTAGATCAGTTGACTTTTAATCAATTGGTCGCTGGTTCGAATCCTTCACGACCAACCAGTTTTAACAACAAACTCAGATGTTAAGCGTGAAGGATAACGTTgc                                                                                                                                             <  1:1514892/150‑1 (MQ=14)

            ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcg                                                                                                                                            >  2:1310716/1‑150 (MQ=255)

              ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTAGAATCCTTCACGACCCACCAGTTTTAAAATAAAAATCAGATGTTAATAGTGAATGATAAAGTTGCGcc                                                                                                                                          >  1:817927/1‑150 (MQ=14)

              ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCCGTTTTAACGTCAAACGCGGATGTTAACCGTGAAGGATAACGTTGCGcc                                                                                                                                          >  2:2435956/1‑150 (MQ=14)

              ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGcc                                                                                                                                          >  2:1688601/1‑150 (MQ=255)

                cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACCCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAg                                                                                                                                        <  1:150220/150‑1 (MQ=255)

                                caccacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTCACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgc                                                                                                                        <  2:1970812/150‑1 (MQ=255)

                                 accacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGTTGTTAAGCGTGAAGGATAACTATGCGCAAGAAACGGCCCGTAGGccg                                                                                                                       >  2:2347321/1‑150 (MQ=14)

                                        gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGag‑cgaag                                                                                                               <  2:2565141/150‑1 (MQ=255)

                                                  cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATcc                                                                                                    >  1:1368789/1‑150 (MQ=255)

                                                   tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCt                                                                                                   >  1:1754834/1‑150 (MQ=255)

                                                     aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGAAGTTAACCGTGACGGATAACAATGCCCACGCAACGGCCCGTAGGGAGAG‑CGAAG‑CTAGTCATcaagc                                                                                                 >  2:1223681/1‑146 (MQ=11)

                                                     aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAAATAAGATGATAAGCGTGAAGGATAACTATTCGCAAGAAACGGCACGTAGGGAGAG‑CGAAG‑CTAGTCATcaagc                                                                                                 >  2:2500462/1‑146 (MQ=14)

                                                          gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGAc                                                                                            >  1:2189168/1‑150 (MQ=255)

                                                          gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAg                                                                                                                           >  1:998819/1‑121 (MQ=38)

                                                                   ttGACTTTTAATCAATTGGTCGCTGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCCGAAGTCGAGTCATACTGCACGACccaccag                                                                                     >  2:2185106/1‑149 (MQ=14)

                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATc                                                                                   >  1:520637/1‑150 (MQ=255)

                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATc                                                                                   >  2:1882846/1‑150 (MQ=255)

                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATc                                                                                   >  2:1455864/1‑150 (MQ=255)

                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAACCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCaaac                                                                                   >  2:2573022/1‑148 (MQ=255)

                                                                   ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAACACGGAAGATAAGCGTTAAGGAAAAAGTAGCGACAGCAACGGAACGAATGGCGCG‑CGAAG‑AGAGTCATCCTGCACGACCCACCaaac                                                                                   >  2:888237/1‑148 (MQ=11)

                                                                     gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCtt                                                                                 <  2:2189168/150‑1 (MQ=255)

                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCTAACACAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGag‑cgaag‑c                                                                                                             >  1:547121/1‑119 (MQ=255)

                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGag‑cgaag‑c                                                                                                             <  2:547121/119‑1 (MQ=255)

                                                                        ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                              >  1:297037/1‑150 (MQ=255)

                                                                         tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCGGCACGACCCACCAATCTTAAAg                                                                             >  2:303302/1‑150 (MQ=255)

                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCGGATGTAAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGACG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                        >  1:1257427/1‑150 (MQ=255)

                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAtttg                                                                        >  2:387597/1‑148 (MQ=255)

                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                        >  1:1721419/1‑150 (MQ=255)

                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                        >  1:2197927/1‑150 (MQ=255)

                                                                              tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                        >  2:1794511/1‑150 (MQ=255)

                                                                                aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAAACTAAAAGATTGGcc                                                                      >  2:190011/1‑150 (MQ=255)

                                                                                aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATACTGCACCACCCACCAATCATAAAGGTTGTcc                                                                      >  2:500544/1‑150 (MQ=14)

                                                                                    ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                  <  1:1729878/150‑1 (MQ=255)

                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTGAGTTGGTAAGCGTGAAGGTTAACGATGAGCAAGCAACGGCCCGGAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGGCCCACCAATCTTAAAGAGTGGCCCCGAGt                                                                >  1:1673722/1‑150 (MQ=17)

                                                                                      tCGCAGGTTCGAATCCTGCACGACCCACCAGTTATAACATCAAACTCAGATGATAATCGAGAAGGATAACGGTGCGCCAGCAACGGCCCGGAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATGGGCCCCGAgg                                                                >  1:170037/1‑149 (MQ=17)

                                                                                       cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTa                                                               <  1:769988/150‑1 (MQ=255)

                                                                                       cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTa                                                               <  2:2092390/150‑1 (MQ=255)

                                                                                         cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAGAGATTGGCCCCGAGTaaa                                                             >  2:1132412/1‑150 (MQ=255)

                                                                                         cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATACAGCACGACACCCCAAACGTAGAGATTTGAACCCAGTaaa                                                             >  2:1650726/1‑150 (MQ=14)

                                                                                          aGGTTCGAATCCTGCACGACCCACCTTTTTTAACATCAAACTAAGATGTTAAGCTTGATGTTTATGGTTTCGACAGCTTCGGCCCGTAGGTCGTG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                                            <  1:2345832/150‑1 (MQ=18)

                                                                                                taaTCGTGCACGACCCACCTGTTTTAACATCAAACTCAGTTGTTATGCGTGTAGGATTACGTTGCGACATCTACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCTCCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                      <  2:981671/149‑1 (MQ=18)

                                                                                                gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                      >  2:1183479/1‑150 (MQ=255)

                                                                                                gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttg                                                      >  2:139059/1‑149 (MQ=255)

                                                                                                 aaTCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggtgg                                                                     >  1:560548/1‑131 (MQ=11)

                                                                                                 aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                     <  1:1688601/150‑1 (MQ=255)

                                                                                                 aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                     <  2:817927/150‑1 (MQ=255)

                                                                                                    ccTGCACGACCCACCAGTTTTAACATCCAAGTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                                  <  1:2434607/150‑1 (MQ=255)

                                                                                                     cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                                 >  1:1051702/1‑150 (MQ=255)

                                                                                                                     ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTc                                >  1:626653/1‑150 (MQ=255)

                                                                                                                     ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTc                                >  2:740961/1‑150 (MQ=255)

                                                                                                                     ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACACGTATGGGTc                                >  2:473599/1‑150 (MQ=255)

                                                                                                                       ttACCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTACAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGt                              <  2:1001821/150‑1 (MQ=255)

                                                                                                                       ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGt                              <  1:1310716/150‑1 (MQ=255)

                                                                                                                        tAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGtt                             <  2:1212497/150‑1 (MQ=255)

                                                                                                                          aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACGACCCATACGGGtgttac                            <  2:1368789/150‑6 (MQ=21)

                                                                                                                           cATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CTAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCACGAGTAAAAATCTTTCAGGTAAC‑ACACGTATGGGTCGTTAGc                          >  1:1100980/1‑150 (MQ=255)

                                                                                                                             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTc                        >  1:1619506/1‑150 (MQ=255)

                                                                                                                             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTc                        >  1:1619507/1‑150 (MQ=255)

                                                                                                                             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTc                        >  1:1619510/1‑150 (MQ=255)

                                                                                                                                aagacAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGt                     >  1:1356001/5‑150 (MQ=255)

                                                                                                                                aagacAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGAATGGGTCGTTAGCTCAGt                     >  1:379211/5‑150 (MQ=255)

                                                                                                                                aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGt                     <  1:1011/150‑1 (MQ=255)

                                                                                                                                    cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGt                 <  1:1455864/150‑1 (MQ=255)

                                                                                                                                    cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGt                 <  1:1882846/150‑1 (MQ=255)

                                                                                                                                     aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                  <  1:1372557/101‑1 (MQ=255)

                                                                                                                                     aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                  >  2:1372557/1‑101 (MQ=255)

                                                                                                                                     aGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTa                >  1:276401/1‑150 (MQ=255)

                                                                                                                                      gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTag               >  2:1783357/1‑150 (MQ=255)

                                                                                                                                       aTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTaga              <  1:1000327/150‑1 (MQ=255)

                                                                                                                                       aTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTaga              >  1:1951831/1‑150 (MQ=255)

                                                                                                                                          ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  1:866461/150‑1 (MQ=255)

                                                                                                                                          ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  1:2295927/150‑1 (MQ=255)

                                                                                                                                          ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  1:2296002/150‑1 (MQ=255)

                                                                                                                                           tAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTaga              >  2:1775106/1‑146 (MQ=255)

                                                                                                                                           tAAGCGTGTAGGATAACGTTGCGCCAGCACCGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAGTGGCCCCGAGTAAAAATCTTTCAGGTCAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAg          <  2:276401/150‑1 (MQ=255)

                                                                                                                                            aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGagtagt         <  2:6885/150‑1 (MQ=255)

                                                                                                                                            aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt         <  1:1652638/150‑1 (MQ=255)

                                                                                                                                            aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt         >  2:1157427/1‑150 (MQ=255)

                                                                                                                                            aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTGGGTAGAGCAGt         >  2:911906/1‑150 (MQ=255)

                                                                                                                                              gCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAACATTCAGGTAAC‑ACCCGTATGGGTCGTTAGCGAAGTTGGTATAGCAGTTg       >  1:1024659/1‑150 (MQ=37)

                                                                                                                                               cGTGAAGGATAACGTTGGGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  1:998438/150‑1 (MQ=255)

                                                                                                                                               cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  2:115831/150‑1 (MQ=255)

                                                                                                                                               cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  1:619515/150‑1 (MQ=255)

                                                                                                                                               cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  2:428697/150‑1 (MQ=255)

                                                                                                                                               cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  1:998360/150‑1 (MQ=255)

                                                                                                                                               cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  1:1851929/150‑1 (MQ=255)

                                                                                                                                                 tGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt    <  1:435106/150‑1 (MQ=255)

                                                                                                                                                 tGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt    >  2:1908776/1‑150 (MQ=255)

                                                                                                                                                  gAAGGATCACGGTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTACAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   <  1:455325/150‑1 (MQ=255)

                                                                                                                                                  gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   <  1:1283116/150‑1 (MQ=255)

                                                                                                                                                  gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   <  1:1890295/150‑1 (MQ=255)

                                                                                                                                                  gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   <  1:2026494/150‑1 (MQ=255)

                                                                                                                                                   aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt  <  1:1357431/150‑1 (MQ=255)

GAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781107‑781403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`