BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,520,516	A→T	5.9%	I12N (ATT→AAT)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,520,516	0	A	T	5.9%	555.5 / 17.5	170	I12N (ATT→AAT)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand): ref base A (82/77); new base T (5/5); total (88/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 174 total aligned reads displayed.

CCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGG  >  NC_002947/5520367‑5520662

ccAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGaa                                                                                                                                                    >  1:4877986/1‑150 (MQ=255)

ccAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGaa                                                                                                                                                    >  2:3453978/1‑150 (MQ=255)

 cAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAAt                                                                                                                                                   >  1:763232/1‑150 (MQ=255)

       gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                             <  1:3648300/150‑1 (MQ=255)

       gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                             >  1:2687864/1‑150 (MQ=255)

       gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                             >  1:4613611/1‑150 (MQ=255)

         acacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGt                                                                                                                                           <  1:2222949/150‑1 (MQ=255)

          cacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTc                                                                                                                                          >  2:429071/1‑150 (MQ=255)

           acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCt                                                                                                                                         <  2:477047/150‑1 (MQ=255)

           acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCt                                                                                                                                         <  1:1615249/150‑1 (MQ=255)

           acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCt                                                                                                                                         >  1:4823837/1‑150 (MQ=255)

           acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCt                                                                                                                                         <  2:1491904/150‑1 (MQ=255)

            ccAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGACGCAGTCTg                                                                                                                                        >  1:4360702/1‑150 (MQ=255)

              agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGcc                                                                                                                                      >  1:3887025/1‑150 (MQ=255)

               gcagACCGGCCAGCAGGCTGGGCAGCAACTTCTGGGGCAATGCCAGCCGCAGCCCGGCCATGACGATGATCGGGCTGGCCAGGGGGATCTGCCGGGGGGTACAGGTCATACAGCCTCGACGGCAGACAGCGGGCCGGCAGCAGTCTGccc                                                                                                                                     >  2:1125237/1‑149 (MQ=255)

                    ccGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTACAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCtt                                                                                                                                >  2:3355778/1‑150 (MQ=255)

                    ccGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCtt                                                                                                                                <  2:4835011/150‑1 (MQ=255)

                      ggCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                              >  1:3476363/1‑150 (MQ=255)

                      ggCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGACAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                              >  1:2624506/1‑150 (MQ=255)

                       gCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggc                                                                                                                             >  2:3854803/1‑150 (MQ=255)

                           gcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGgcggcg                                                                                                                         >  2:4496713/1‑150 (MQ=255)

                              ggCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                      <  2:3481452/150‑1 (MQ=255)

                               gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAgg                                                                                                                     >  2:1079281/1‑150 (MQ=255)

                               gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAgg                                                                                                                     >  1:4221921/1‑150 (MQ=255)

                                    gcagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                                <  1:1336550/150‑1 (MQ=255)

                                    gcagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                                <  1:2047567/150‑1 (MQ=255)

                                      agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgc                                                                                                              >  2:1712788/1‑150 (MQ=255)

                                       gcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgcg                                                                                                             >  1:4089265/1‑150 (MQ=255)

                                          aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAg                                                                                                          <  2:4731458/150‑1 (MQ=255)

                                          aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAg                                                                                                          <  1:1025971/150‑1 (MQ=255)

                                            ttCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAggg                                                                                                        >  2:5167807/1‑149 (MQ=255)

                                            ttCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGt                                                                                                        >  2:4411147/1‑150 (MQ=255)

                                            ttCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGt                                                                                                        >  2:2837349/1‑150 (MQ=255)

                                              cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATTAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcg                                                                                                                            >  2:1574587/1‑128 (MQ=255)

                                              cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATTAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcg                                                                                                                            <  1:1574587/128‑1 (MQ=255)

                                               aGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCAt                                                                                                     <  2:2624506/150‑1 (MQ=255)

                                                 gggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCAttt                                                                                                   >  1:4861278/1‑150 (MQ=255)

                                                 gggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCAttt                                                                                                   >  2:3248605/1‑150 (MQ=255)

                                                    cAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                                <  1:4206972/150‑1 (MQ=255)

                                                      aaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttct                                                                                              >  1:3930945/1‑150 (MQ=255)

                                                       aGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttctt                                                                                             >  1:1405738/1‑150 (MQ=255)

                                                         ccAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                           >  1:4935149/1‑150 (MQ=255)

                                                         ccAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                           <  2:2687914/150‑1 (MQ=255)

                                                           agcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAgg                                                                                         >  1:4374941/1‑150 (MQ=255)

                                                              agcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                               <  1:497543/125‑1 (MQ=255)

                                                              agcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                               >  2:497543/1‑125 (MQ=255)

                                                                   cccGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCt                                                                                 >  2:3051850/1‑150 (MQ=255)

                                                                     cGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAg                                                                                                          <  2:4584636/123‑1 (MQ=255)

                                                                     cGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCa                                                                               >  2:2210805/1‑150 (MQ=255)

                                                                        ccATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGa                                                                            <  1:4992974/150‑1 (MQ=255)

                                                                          aTGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAg                                                                          <  1:4673320/150‑1 (MQ=255)

                                                                              cGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCt                                                                      >  1:4532917/1‑150 (MQ=255)

                                                                               gatgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCt                                                                                                                                         <  1:543122/82‑1 (MQ=37)

                                                                               gatgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCt                                                                                                                                         >  2:543122/1‑82 (MQ=37)

                                                                                 tgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                             >  1:2634642/1‑76 (MQ=255)

                                                                                 tgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                   >  2:5351188/1‑150 (MQ=255)

                                                                                 tgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                   <  2:4121569/150‑1 (MQ=255)

                                                                                 tgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                   <  1:1432823/150‑1 (MQ=255)

                                                                                   atCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgcgc                                                                 <  2:2670314/150‑1 (MQ=255)

                                                                                    tCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCt                                                                <  2:5301633/150‑1 (MQ=255)

                                                                                    tCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCt                                                                <  2:1171925/150‑1 (MQ=255)

                                                                                      aGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCa                                                              <  2:3743933/150‑1 (MQ=255)

                                                                                        gCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCa                                                                                                      >  2:1965798/1‑108 (MQ=255)

                                                                                         cTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATcc                                                           <  1:10920/150‑1 (MQ=255)

                                                                                            gCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTg                                                                                                                                        <  1:1691925/70‑1 (MQ=25)

                                                                                            gCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATcccc                                                         >  1:5400724/1‑149 (MQ=255)

                                                                                            gCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATcccc                                                         >  1:5400605/1‑149 (MQ=255)

                                                                                            gCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATcccc                                                         <  2:5400605/149‑1 (MQ=255)

                                                                                              cAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCtt                                                      <  1:3787129/150‑1 (MQ=255)

                                                                                                        tGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCt                                                       >  2:1656605/1‑139 (MQ=255)

                                                                                                        tGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCt                                                       <  1:1656605/139‑1 (MQ=255)

                                                                                                            ggggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCCTCGGCGCAGGTCATGTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCGTtggtgg                                        >  2:2517009/1‑150 (MQ=255)

                                                                                                            ggggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCCGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtgg                                        >  2:223203/1‑150 (MQ=255)

                                                                                                            ggggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtgg                                        >  2:3112958/1‑150 (MQ=255)

                                                                                                             gggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                       <  1:286286/150‑1 (MQ=255)

                                                                                                             gggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                       <  1:286211/150‑1 (MQ=255)

                                                                                                                     aGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                               >  2:1933621/1‑150 (MQ=255)

                                                                                                                          aTACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                          >  1:795471/1‑150 (MQ=255)

                                                                                                                          aTACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                          >  1:3891526/1‑150 (MQ=255)

                                                                                                                           tACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGc                         <  2:951411/150‑1 (MQ=255)

                                                                                                                           tACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGc                         <  1:5197277/150‑1 (MQ=255)

                                                                                                                            aCAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGAGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTTCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCa                        <  2:4953510/150‑1 (MQ=255)

                                                                                                                             cAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAg                       <  2:3299913/150‑1 (MQ=255)

                                                                                                                             cAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAg                       <  2:410836/150‑1 (MQ=255)

                                                                                                                             cAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAg                       <  2:1629162/150‑1 (MQ=255)

                                                                                                                                 cTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCt                   >  1:4346475/1‑150 (MQ=255)

                                                                                                                                  tCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAACATTGCAGGTCtt                  >  1:4011721/1‑150 (MQ=255)

                                                                                                                                  tCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTATGCGGCGCAGGCATCGGAGCAGGTCATTTATTCTTCAGGAAAACGATCATGAATGCCTTGCGCTCATCACCCTTCAGGGCTTTGGTGGTGGCGTAAGCATTGCAGGTCtt                  >  2:4612258/1‑150 (MQ=255)

                                                                                                                                   cGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTc                 >  1:4589370/1‑150 (MQ=255)

                                                                                                                                        gCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttt            >  1:830222/1‑150 (MQ=255)

                                                                                                                                        gCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttt            >  1:5342211/1‑150 (MQ=255)

                                                                                                                                        gCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttt            >  2:1471245/1‑150 (MQ=255)

                                                                                                                                          aGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtg                                      <  1:3389964/122‑1 (MQ=255)

                                                                                                                                           gACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTc         <  2:4877986/150‑1 (MQ=255)

                                                                                                                                           gACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTc         <  1:4394589/150‑1 (MQ=255)

                                                                                                                                            aCAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCt        >  1:4901294/1‑150 (MQ=255)

                                                                                                                                             cAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCtt       <  2:3656246/150‑1 (MQ=255)

                                                                                                                                                cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCt    <  2:4221921/150‑1 (MQ=255)

                                                                                                                                                 gggAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTg   >  1:2948731/1‑150 (MQ=255)

                                                                                                                                                  ggAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTgg  <  1:2837349/150‑1 (MQ=255)

CCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGG  >  NC_002947/5520367‑5520662

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