BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,054,905	T→G	8.6%	Q180H (CAA→CAC)	ygfB ←	UPF0149 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,054,905	0	T	G	8.6%	1023.4 / 23.4	337	Q180H (CAA→CAC)	ygfB	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base T (93/208); new base G (34/2); total (127/210)
Fisher's exact test for biased strand distribution p-value = 5.20e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

AAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGCTTCACC  >  NC_000913/3054756‑3055054

 aaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGtt                                                                                                                                                       <  1:1650636/150‑1 (MQ=255)

 aaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGtt                                                                                                                                                       <  1:4753198/150‑1 (MQ=255)

    cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTcgcg                                                                                                                                                    <  2:1708331/150‑1 (MQ=255)

    cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                                    <  1:4218934/150‑1 (MQ=255)

    cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                                    <  1:1224848/150‑1 (MQ=255)

    cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                                    <  2:3716807/150‑1 (MQ=255)

       cgcgGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGAt                                                                                                                                                 <  1:4471674/150‑1 (MQ=255)

         cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                                               <  2:4803186/150‑1 (MQ=255)

         cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                                               <  1:5340236/150‑1 (MQ=255)

         cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                                               <  2:1115503/150‑1 (MQ=255)

             gcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCAGGCGAGGTCGGTTGCGGCTGAGGaa                                                                                                                                           >  2:351593/1‑150 (MQ=255)

              cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaaa                                                                                                                                          <  2:3156954/150‑1 (MQ=255)

               tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                                         <  1:3643865/150‑1 (MQ=255)

                   gggTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGGGCGGTCGGTTGCGGGTGAGTAAAGgtgt                                                                                                                                     >  1:2643199/1‑150 (MQ=255)

                   gggTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGGCGGTTGCGGATGAGTAAAGgttt                                                                                                                                     >  1:1357458/1‑148 (MQ=255)

                         cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTCTGGCGCGGTCGGTTGCGGATTAGTAAAGGTGTCGTGGc                                                                                                                               >  2:742521/1‑150 (MQ=255)

                         cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCTGGCGCGGGCGGTTGCGGATGTGTAACGGGGGCGTGGc                                                                                                                               >  2:5231496/1‑150 (MQ=255)

                           tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTAATTCTTGCGCGGTGGGTTTCGGATGGGTAAAGGTGTGGTGGCAt                                                                                                                             >  2:3807280/1‑150 (MQ=255)

                            ttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACGTCTGGCGCGGTCGGTTGCGGGTGAGTTAAGGTGGCGTGGCATa                                                                                                                            >  2:2173839/1‑150 (MQ=255)

                                tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTCTTCTGGCGCGGGCGGGTGCGGATGGGTAAAGGTGTCGTGGCCTaaacg                                                                                                                        >  1:5140658/1‑147 (MQ=255)

                                 ccaccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAgc                                                                                                                       <  1:5237234/150‑1 (MQ=255)

                                   accaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCTTGCGCGGTAGGTTGCGGGTGTGTAAAGGTGTAGTGGCATAAAAgcgc                                                                                                                     >  2:412191/1‑150 (MQ=255)

                                    ccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTCGTCTGGCGCGGGCGGTTTCGGATGGGGAAAGGGGGGGGGGCAGAACAGcgcg                                                                                                                    >  1:2102179/1‑150 (MQ=255)

                                      aGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTACTTCTGGCGCGGGCGGGTGCGGATGAGGAAAGGTGTCGGGGCATAAAAGCGCGGc                                                                                                                  >  2:162020/1‑150 (MQ=255)

                                        ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCaa                                                                                                                <  2:4907867/150‑1 (MQ=255)

                                         gCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAATGGTGTCGTGGCATAACAGCGCGGCAac                                                                                                               <  1:5439563/150‑1 (MQ=255)

                                            tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTGTGGCGCGGTCGGGTGGGGGTGAGTAAAGGTGTGGTGGCATAACAGCGCGGCAaca‑cc                                                                                                           >  2:1514212/1‑149 (MQ=255)

                                               cGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAAc                                                                                                        <  2:4098982/150‑1 (MQ=255)

                                                     cGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTc                                                                                                  <  2:1031362/150‑1 (MQ=255)

                                                     cGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTc                                                                                                  <  2:1818222/150‑1 (MQ=255)

                                                       cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGG‑TTTTGTTCGTCTGGCGCGGTGGGTTTCTGATTATTAAAGGTGTCGTGGCATAAAACCGCGGCAACA‑CTAAAGAATACgat                                                                                               >  1:3110425/1‑150 (MQ=255)

                                                        tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgat                                                                                               <  2:2268423/150‑1 (MQ=255)

                                                        tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgat                                                                                               <  1:1006666/150‑1 (MQ=255)

                                                        tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgat                                                                                               <  1:1210994/150‑1 (MQ=255)

                                                         gAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGGACTTCAGGCGCGGGCGGGTGCGGGTGGGTTAATGGGTCGGGGCATAACACCGCGGCAACA‑CCAACGTATGCgagg                                                                                              >  1:2236376/1‑148 (MQ=255)

                                                          aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgatga                                                                                             <  1:4272962/150‑1 (MQ=255)

                                                          aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgatga                                                                                             <  1:2409274/150‑1 (MQ=255)

                                                          aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgatga                                                                                             <  1:850982/150‑1 (MQ=255)

                                                           aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCgatgat                                                                                            <  2:3103700/150‑1 (MQ=255)

                                                              tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTGTTCTTGCGCGGGGGGGTGCGGGTGAGTAAATGTGGGGTGGCATAAAAGCGCGGCCACACCCAACG‑ATTCgatgttata                                                                                         >  1:2344112/1‑145 (MQ=255)

                                                              tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTATTATGGCGCGGTCGGTTGCGGAGGAGTAAATGTTTGGTGGCCAAAAACCCCGGCAACA‑CGAAAGTATTCGATGAtatc                                                                                         >  2:4701382/1‑150 (MQ=255)

                                                                ttGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCtt                                                                                       <  1:1991048/150‑1 (MQ=255)

                                                                ttGCCGGGATATCTCACTCATAACACTCTCCGTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCtt                                                                                       <  1:2919287/150‑1 (MQ=255)

                                                                   ccGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTGTTGCGCGGTGGGTTTCGGATTTGTAAATGTGTGGTTGCAATAAAGCCCGGCAAAA‑CGAAACGATTAGCTGATATCTTaaa                                                                                    >  2:3359480/1‑150 (MQ=255)

                                                                     gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGc                                                                                  <  1:4982949/150‑1 (MQ=255)

                                                                         tatCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGCGGTCGGTTGCGGCTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CCAACCTATTAGATGATCTCTTCAAGCGACa                                                                              >  2:2840350/1‑150 (MQ=255)

                                                                             tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTc                                                                          <  1:3524403/150‑1 (MQ=255)

                                                                             tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTc                                                                          <  1:2967258/150‑1 (MQ=255)

                                                                               actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                                        <  1:2774187/150‑1 (MQ=255)

                                                                                ctcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAg                                                                       <  1:576287/150‑1 (MQ=255)

                                                                                ctcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAg                                                                       <  1:3374029/150‑1 (MQ=255)

                                                                                ctcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAg                                                                       <  1:4149596/150‑1 (MQ=255)

                                                                                 tcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGc                                                                      <  2:562540/150‑1 (MQ=255)

                                                                                 tcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCCGCAACA‑CCAACGTATTAGATGATCACATCAAGCGACAATTCaaac                                                                      >  2:3937557/1‑148 (MQ=255)

                                                                                   aTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGctc                                                                    <  2:2010943/150‑1 (MQ=255)

                                                                                    tAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCTGGCGCGGTCGGTTGCGGATGAGTAAATGGGTCGTGGCATAACAGCGCGGCAACA‑CCAAAGTATTAGATGATCTCTTAAAGCGACATTTCAAGctat                                                                   >  2:1186217/1‑148 (MQ=255)

                                                                                        acTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCt                                                               <  2:5097380/150‑1 (MQ=255)

                                                                                         ctctCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTg                                                              <  2:840300/150‑1 (MQ=255)

                                                                                          tctcCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGa                                                             <  2:1067107/150‑1 (MQ=255)

                                                                                          tctcCTGACGTCTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGCGGTCGGGTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTTTTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTgt                                                             <  2:4627326/150‑2 (MQ=255)

                                                                                             ccTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtct                                                          <  2:273034/150‑1 (MQ=255)

                                                                                               ttACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttc                                                        <  1:3207602/150‑1 (MQ=255)

                                                                                               ttACGTTTTTTGTTTTTAGTGTAGAGTCGGATTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttc                                                        <  1:3844780/150‑1 (MQ=255)

                                                                                                     tttttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtcttc                                                  <  1:3786916/150‑1 (MQ=255)

                                                                                                       tttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATcttcgtcttcgt                                                <  2:5152558/150‑1 (MQ=255)

                                                                                                        ttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCTGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCttcctcttcgtc                                               <  1:3633104/150‑1 (MQ=255)

                                                                                                        ttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCttcgtcttcgtc                                               <  1:2162675/150‑1 (MQ=255)

                                                                                                          gTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTcgtcgt                                             <  2:3306977/150‑1 (MQ=255)

                                                                                                            ttttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTaa                                           <  1:2504607/150‑1 (MQ=255)

                                                                                                             tttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAc                                          <  2:2964687/150‑1 (MQ=255)

                                                                                                               tAGTGTAGAGTCGGTTTTTTTACTTATTGCGCGGTCGGGTGCGGATGAGTAAAGGTGTCGTGGCATAACACCGCGGCAACA‑AGAAACTATTCGATGTTCTCTTAAAGCGACATTTCAAGCTCTTACGGATCTTCGTCTTCGTCGTAAccc                                        >  2:781938/1‑150 (MQ=255)

                                                                                                                  tgtAGAGTCGGTTTTTGTAATTAAGTCGCGGTCGGGTGAGGATGAGTAAAGGTGTCGTGGCAAAAAAACCCGGCAACA‑CGAACCTATTCGATTATCACTTCAAGCGACATTTAAATCAATTCCTGATCATCGGCATCGTCGTAAcacact                                     >  2:4585741/1‑148 (MQ=255)

                                                                                                                     agagTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgc                                  <  2:3475811/150‑1 (MQ=255)

                                                                                                                         tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa                              >  2:1857242/1‑150 (MQ=255)

                                                                                                                              ttttGTACTTCTGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAATCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTAGTGACCCAGTTGCGGAATGttt                         >  1:1011443/1‑149 (MQ=255)

                                                                                                                                 tGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCACTGTTACGc                      >  2:1334086/1‑150 (MQ=255)

                                                                                                                                 tGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGc                      >  1:3441298/1‑150 (MQ=255)

                                                                                                                                  gTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCa                     >  2:4416491/1‑150 (MQ=255)

                                                                                                                                     cTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAt                  <  1:1714198/150‑1 (MQ=255)

                                                                                                                                     cTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAt                  <  2:1617540/150‑1 (MQ=255)

                                                                                                                                     cTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAt                  <  2:423304/150‑1 (MQ=255)

                                                                                                                                     cTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAt                  >  2:5318372/1‑150 (MQ=255)

                                                                                                                                     cTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAt                  <  2:2771964/150‑1 (MQ=255)

                                                                                                                                     cTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAt                  <  1:2814115/150‑1 (MQ=255)

                                                                                                                                      ttCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAtc                 >  1:4884805/1‑150 (MQ=255)

                                                                                                                                        cTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAtcgt               <  2:2663041/150‑1 (MQ=255)

                                                                                                                                        cTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAtcgt               <  2:2356268/150‑1 (MQ=255)

                                                                                                                                        cTGGCGCGGTCGGGTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAtcgt               <  2:2570305/150‑1 (MQ=255)

                                                                                                                                         tGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGAtcgtc              <  2:1446351/150‑1 (MQ=255)

                                                                                                                                         tGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCATCCAGATCTTCGTCTTCGTCGTACCCCAGTTGCGCAATGTTACGCAGAtcgtc              >  1:2569543/1‑150 (MQ=255)

                                                                                                                                           gcgcGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGa            <  1:4954288/150‑1 (MQ=255)

                                                                                                                                           gcgcGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGa            >  2:1380262/1‑150 (MQ=255)

                                                                                                                                            cgcgGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACATCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCACTTCCAGATATTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTAAGCAGATAGTCGAt           >  1:4037033/1‑150 (MQ=255)

                                                                                                                                            cgcgGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGAt           <  2:4353861/150‑1 (MQ=255)

                                                                                                                                             gcgGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATa          <  2:89102/150‑1 (MQ=255)

                                                                                                                                               ggtcggtTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGc        <  1:4382808/150‑1 (MQ=255)

                                                                                                                                                   ggtTGCGGATGAGTAATGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGCTTCa    >  2:5296327/1‑150 (MQ=255)

                                                                                                                                                   ggtTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGCTTCa    >  1:183767/1‑150 (MQ=255)

                                                                                                                                                   ggtTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGCTTCa    >  1:151986/1‑150 (MQ=255)

                                                                                                                                                     ttGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGCTTCAcc  >  2:4037584/1‑150 (MQ=255)

AAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACA‑CGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGATCGTCGATAGCTTCACC  >  NC_000913/3054756‑3055054

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