BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,269	C→T	7.1%	intergenic (+47/‑100)	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,269	0	C	T	7.1%	341.4 / 13.9	115	intergenic (+47/‑100)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base C (52/50); new base T (6/2); total (61/54)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 122 total aligned reads displayed.

TTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGC  >  NC_000913/781121‑781418

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

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

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

                 ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACTTCTAACTCAGATGTTAAGATTGAAGGATAACGTTGCGCCAGCAACGGCCGGTTggg                                                                                                                                        >  2:1510757/1‑150 (MQ=14)

                  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                                                                                                  >  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)

                                                           tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACTTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATACTGCACGACCCACCAATCTTAAAg                                                                                            >  2:2566831/1‑150 (MQ=255)

                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                                        >  2:435106/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:2197927/1‑150 (MQ=255)

                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                                       >  1:1721419/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)

                                                                        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)

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

                                                                                 cGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                                      <  1:571955/150‑1 (MQ=255)

                                                                                  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                                               >  2:740961/1‑150 (MQ=255)

                                                                                                       ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTc                                               >  2:104152/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‑ACCCGTATGGGTCGTTAGCTCAGt                                    >  1:1356015/5‑150 (MQ=255)

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

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

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

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

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

                                                                                                                       aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                                 <  1:1372557/101‑1 (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:2295927/150‑1 (MQ=255)

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

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

                                                                                                                             tAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTaga                             <  1:1775106/146‑1 (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:85546/150‑1 (MQ=255)

                                                                                                                              aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                        <  1:1652638/150‑1 (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                     <  1:1851929/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                                                     taGGATACCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTTCACGACCCACCAATCTTAAAGATTTGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACACGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                 <  1:1984445/149‑1 (MQ=255)

                                                                                                                                           aaCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCa           >  1:1696344/1‑150 (MQ=255)

                                                                                                                                           aaCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTGAATCa           >  1:1952719/1‑150 (MQ=255)

                                                                                                                                           aaCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAAACTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCa           >  1:394530/1‑150 (MQ=255)

                                                                                                                                            aCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa          <  1:1794511/150‑1 (MQ=255)

                                                                                                                                            aCGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa          <  1:1794512/150‑1 (MQ=255)

                                                                                                                                             cGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt         <  1:2043306/150‑1 (MQ=255)

                                                                                                                                              gTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt        <  1:450477/150‑1 (MQ=255)

                                                                                                                                               ttGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTg       >  1:885541/1‑150 (MQ=255)

                                                                                                                                                  cgcCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTc    >  1:1429086/1‑150 (MQ=255)

                                                                                                                                                   gcCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg   >  2:360121/1‑150 (MQ=255)

                                                                                                                                                    ccAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAGTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑GCCCGTCTGGGTCGTTAGAGCAGTGGGTAGAGCGGTTGACTGTTAAGCAATTGGTCGc  >  2:49970/1‑150 (MQ=37)

TTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGC  >  NC_000913/781121‑781418

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