BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,282,339	G→A	6.0%	H258H (CAC→CAT)	ampC ←	beta‑lactamase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,282,339	0	G	A	6.0%	779.8 / 26.2	248	H258H (CAC→CAT)	ampC	beta‑lactamase
Reads supporting (aligned to +/- strand): ref base G (105/127); new base A (9/6); total (115/133)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 259 total aligned reads displayed.

TGTTGCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAG  >  NC_002947/3282191‑3282488

 gTTGCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCAtgtg                                                                                                                                                     >  2:354545/1‑150 (MQ=255)

  ttGCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGc                                                                                                                                                    <  2:3204607/150‑1 (MQ=255)

     ccATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAAc                                                                                                                                                 >  1:2231061/1‑150 (MQ=255)

      cATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAAcc                                                                                                                                                >  2:3782623/1‑150 (MQ=255)

      cATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGGTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAAcc                                                                                                                                                <  2:4837082/150‑1 (MQ=255)

       aTCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCg                                                                                                                                               >  2:4014858/1‑150 (MQ=255)

        tCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGc                                                                                                                                              >  2:1136909/1‑150 (MQ=255)

         cGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGTTGCACATGCAGCCGca                                                                                                                                             <  1:5175494/150‑1 (MQ=255)

          gATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGTTGCACATGCAGCCGcag                                                                                                                                            <  1:1000488/150‑2 (MQ=255)

          gATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGcac                                                                                                                                            >  2:1965545/1‑150 (MQ=255)

           aTCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGcaca                                                                                                                                           <  1:4367236/150‑1 (MQ=255)

           aTCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGcaca                                                                                                                                           <  2:2833063/150‑1 (MQ=255)

                    cTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                  <  2:3738838/150‑1 (MQ=255)

                    cTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                  <  2:3107691/150‑1 (MQ=255)

                    cTCAGCGTGACAGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                  >  1:3970005/1‑150 (MQ=255)

                      cAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGcagcag                                                                                                                                >  2:2323783/1‑150 (MQ=255)

                       aGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAgg                                                                                                                               <  1:5368747/150‑1 (MQ=255)

                       aGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAgg                                                                                                                               <  2:3490933/150‑1 (MQ=255)

                         cGTGACCGGGTAGGGGTAGCGTTCCCAGCCCGGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTc                                                                                                                             <  2:5110642/150‑1 (MQ=255)

                         cGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTc                                                                                                                             >  1:1405730/1‑150 (MQ=255)

                           tGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgc                                                                                                                           >  1:2357589/1‑150 (MQ=255)

                            gACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgct                                                                                                                          <  1:5042755/150‑1 (MQ=255)

                             aCCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctg                                                                                                                         >  2:4601520/1‑150 (MQ=255)

                               cGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggc                                                                                                                       <  1:1617319/150‑1 (MQ=255)

                                 ggtaggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgc                                                                                                                     <  2:1704713/150‑1 (MQ=255)

                                  gtaggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgct                                                                                                                    <  1:891364/150‑1 (MQ=255)

                                   taggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgctg                                                                                                                   <  1:4587277/150‑1 (MQ=255)

                                   taggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgctg                                                                                                                   <  1:1211905/150‑1 (MQ=255)

                                     ggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGt                                                                                                                 <  2:1086351/150‑1 (MQ=255)

                                      gggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtt                                                                                                                <  1:985582/150‑1 (MQ=255)

                                      gggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtt                                                                                                                >  1:5392892/1‑150 (MQ=255)

                                      gggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtt                                                                                                                >  2:5444603/1‑150 (MQ=255)

                                       ggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGttt                                                                                                               <  2:280241/150‑1 (MQ=255)

                                        gtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtttt                                                                                                              <  2:862199/150‑1 (MQ=255)

                                        gtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtttt                                                                                                              >  1:4688403/1‑150 (MQ=255)

                                         tagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTg                                                                                                             <  2:5252844/150‑1 (MQ=255)

                                           gCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGAt                                                                                                           <  1:5403055/150‑1 (MQ=255)

                                           gCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCGTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGAt                                                                                                           <  1:4813985/150‑1 (MQ=255)

                                             gTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATgt                                                                                                         <  2:2476035/150‑2 (MQ=255)

                                              ttCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGcc                                                                                                        >  2:5004303/1‑150 (MQ=255)

                                               tCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCg                                                                                                       >  1:3278605/1‑150 (MQ=255)

                                                cccagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGt                                                                                                      >  1:307526/1‑150 (MQ=255)

                                                 ccagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTa                                                                                                     <  2:73759/150‑1 (MQ=255)

                                                 ccagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTa                                                                                                     >  2:2553457/1‑150 (MQ=255)

                                                  cagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAg                                                                                                    >  1:4451669/1‑150 (MQ=255)

                                                      ccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCt                                                                                                >  1:4998609/1‑150 (MQ=255)

                                                       cagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTc                                                                                               >  1:5426669/1‑150 (MQ=255)

                                                           ccTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCg                                                                                           <  2:1562423/150‑1 (MQ=255)

                                                            cTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCgg                                                                                          >  1:1256883/1‑150 (MQ=255)

                                                             ttGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGc                                                                                         <  1:4198362/150‑1 (MQ=255)

                                                                cGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTa                                                                                      <  2:1932089/150‑1 (MQ=255)

                                                                 gTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTaa                                                                                     <  2:2924252/150‑1 (MQ=255)

                                                                  tCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAg                                                                                    <  1:1751963/150‑1 (MQ=255)

                                                                  tCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGAAACCGCACAAAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGTCGTAAg                                                                                    >  2:5628759/1‑150 (MQ=255)

                                                                    aTGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAggg                                                                                  >  1:2365473/1‑150 (MQ=255)

                                                                         gcCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAgggccggg                                                                             >  1:2295823/1‑150 (MQ=255)

                                                                            aaCCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGccc                                                                          >  1:49538/1‑150 (MQ=255)

                                                                            aaCCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGccc                                                                          >  1:591023/1‑150 (MQ=255)

                                                                             aCCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCt                                                                         <  2:4341732/150‑1 (MQ=255)

                                                                              ccTGGAAGTAGCTCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTa                                                                        <  1:4893964/150‑1 (MQ=255)

                                                                               cTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTac                                                                       <  2:4810825/150‑1 (MQ=255)

                                                                               cTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTac                                                                       <  2:3266864/150‑1 (MQ=255)

                                                                               cTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTac                                                                       >  1:1370121/1‑150 (MQ=255)

                                                                                tGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTaca                                                                      <  1:1966663/150‑1 (MQ=255)

                                                                                     gTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCAc                                                                 >  2:3919757/1‑150 (MQ=255)

                                                                                      gaGCCCTGCTGGGTAAGGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCAcc                                                                <  2:4513321/149‑1 (MQ=255)

                                                                                         cccTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGc                                                             <  1:1793482/150‑1 (MQ=255)

                                                                                         cccTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGc                                                             >  1:5568329/1‑150 (MQ=255)

                                                                                         cccTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGc                                                             <  2:5355915/150‑1 (MQ=255)

                                                                                          ccTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCt                                                            <  2:4451669/150‑1 (MQ=255)

                                                                                            tgctgGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTg                                                          >  2:5583263/1‑150 (MQ=255)

                                                                                             gctgGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTgt                                                         >  1:2930086/1‑150 (MQ=255)

                                                                                              ctgGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTgtc                                                        >  2:2210673/1‑150 (MQ=255)

                                                                                                   tAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCgg                                                   <  1:3664428/150‑1 (MQ=255)

                                                                                                       ggCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTc                                               <  2:3860436/150‑1 (MQ=255)

                                                                                                                tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                      <  2:5426096/150‑1 (MQ=255)

                                                                                                                tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                      >  1:3860993/1‑150 (MQ=255)

                                                                                                                tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                      <  1:2466621/150‑1 (MQ=255)

                                                                                                                tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                      >  2:679946/1‑150 (MQ=255)

                                                                                                                tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                      >  1:4550280/1‑150 (MQ=255)

                                                                                                                 ttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTgg                                     <  2:1370121/150‑1 (MQ=255)

                                                                                                                 ttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTgg                                     <  1:3146603/150‑1 (MQ=255)

                                                                                                                   cGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGc                                   <  2:555651/150‑1 (MQ=255)

                                                                                                                    gCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATCCCCTTGGGCa                                  >  2:229591/1‑150 (MQ=255)

                                                                                                                     ccAGGGCCGGTGGCAGGCCCTGCGGCTGCATGAGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCAt                                 <  2:4122116/150‑1 (MQ=255)

                                                                                                                       aGGGCCGGGGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATaa                               <  2:1480023/150‑1 (MQ=255)

                                                                                                                           ccGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGcc                           >  2:4665428/1‑150 (MQ=255)

                                                                                                                                ggCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTgg                      <  1:155566/150‑1 (MQ=255)

                                                                                                                                 gCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTggg                     >  2:5296403/1‑150 (MQ=255)

                                                                                                                                 gCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTggg                     >  2:3934637/1‑150 (MQ=255)

                                                                                                                                  cgctatgtGCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCACTGGgc                    >  2:268783/8‑150 (MQ=255)

                                                                                                                                  cAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGgc                    <  1:3782623/150‑1 (MQ=255)

                                                                                                                                         tGCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGcc             >  1:2152455/1‑150 (MQ=255)

                                                                                                                                         tGCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCACTGGGCGCTGGcc             <  1:268783/150‑1 (MQ=255)

                                                                                                                                         tGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTGTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGcc             >  2:4826653/1‑150 (MQ=255)

                                                                                                                                           cGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCgg           <  2:1857380/150‑1 (MQ=255)

                                                                                                                                                  aTGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGa    >  1:3405244/1‑150 (MQ=255)

                                                                                                                                                  aTGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATACCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGa    >  1:787818/1‑150 (MQ=255)

                                                                                                                                                   tgtgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGaa   <  1:3930666/150‑1 (MQ=255)

                                                                                                                                                    gtgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAag  <  2:2979693/150‑1 (MQ=255)

TGTTGCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAG  >  NC_002947/3282191‑3282488

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