BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,135,947	T→A	5.4%	Q350L (CAG→CTG)	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,135,947	0	T	A	5.4%	1274.9 / 38.2	372	Q350L (CAG→CTG)	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (218/133); new base A (8/12); total (226/146)
Fisher's exact test for biased strand distribution p-value = 5.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 391 total aligned reads displayed.

ACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGTCG  >  NC_002947/4135798‑4136096

  accacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTgg                                                                                                                                                     >  1:6918710/1‑150 (MQ=255)

  accacTGGGGCAAAGGATGCTCACGTCCTCACCTGGTGACAGCTGCAGGTCGAACTGCGCCAGCACATGCTGCCAACCACCCTCACGCGGGTAGCCCAATCGGATGTCACGGGCCTCCAGCAACACCGCGCGCATGCATCGCATGCCTgg                                                                                                                                                     >  1:2934528/1‑150 (MQ=255)

      cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCCCGCGGGTAGCCCCAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGtt                                                                                                                                                 <  2:6621941/150‑1 (MQ=255)

        ggggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGc                                                                                                                                               <  1:5116454/150‑1 (MQ=255)

         gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGCACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGACGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCAGGCAGGCGTTGCa                                                                                                                                              <  2:3596904/150‑1 (MQ=255)

          ggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCAGGCAGGCGTTGCat                                                                                                                                             <  1:1428376/150‑2 (MQ=255)

          ggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCAGGCAGGCGTTGCat                                                                                                                                             <  1:7215219/150‑2 (MQ=255)

           gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTg                                                                                                                                                <  2:1439435/146‑1 (MQ=255)

            ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGtt                                                                                                                                           >  2:4524844/1‑150 (MQ=255)

             cAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTc                                                                                                                                          >  1:1788465/1‑150 (MQ=255)

             cAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCAGGCAGGCGTTGCatgag                                                                                                                                          <  1:1272559/150‑5 (MQ=255)

             cAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCACGCAACGCCAGGCAGGCGTTGCatgag                                                                                                                                          <  1:8264808/150‑5 (MQ=255)

              aaGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCg                                                                                                                                         <  2:2451712/150‑1 (MQ=255)

                ggATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGgc                                                                                                                                       >  2:2362247/1‑150 (MQ=255)

                     cTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCa                                                                                                                                  >  1:6488620/1‑150 (MQ=255)

                      tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCaa                                                                                                                                 >  2:2748666/1‑150 (MQ=255)

                      tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCaa                                                                                                                                 >  1:256911/1‑150 (MQ=255)

                      tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCaa                                                                                                                                 >  1:87120/1‑150 (MQ=255)

                        aCGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAAtt                                                                                                                               <  1:456327/150‑1 (MQ=255)

                              tCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCa                                                                                                                         >  2:172393/1‑150 (MQ=255)

                                  cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCt                                                                                                                     >  2:1357654/1‑150 (MQ=255)

                                   tGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTg                                                                                                                    >  2:8024760/1‑150 (MQ=255)

                                   tGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTg                                                                                                                    <  2:7610035/150‑1 (MQ=255)

                                      tGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                                 <  2:109461/150‑1 (MQ=255)

                                        ccAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGt                                                                                                               >  2:5703052/1‑150 (MQ=255)

                                        ccAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGt                                                                                                               >  1:8989351/1‑150 (MQ=255)

                                           gCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGAc                                                                                                            >  2:885730/1‑150 (MQ=255)

                                                  gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                     <  1:815042/150‑1 (MQ=255)

                                                   tCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCa                                                                                                    <  1:5333315/150‑1 (MQ=255)

                                                      aaCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGa                                                                                                 >  2:6578893/1‑150 (MQ=255)

                                                       aCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGaa                                                                                                >  2:8721825/1‑150 (MQ=255)

                                                          gcgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCACTTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGc                                                                                             <  1:2331847/150‑1 (MQ=255)

                                                          gcgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGc                                                                                             <  1:9163712/150‑1 (MQ=255)

                                                          gcgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGc                                                                                             <  1:947270/150‑1 (MQ=255)

                                                           cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                            >  1:1002446/1‑150 (MQ=255)

                                                           cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                            >  1:1753080/1‑150 (MQ=255)

                                                             ccAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGa                                                                                          <  2:3140941/150‑1 (MQ=255)

                                                              cAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGAc                                                                                         >  1:5036276/1‑150 (MQ=255)

                                                              cAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGAc                                                                                         >  1:808027/1‑150 (MQ=255)

                                                               aGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACt                                                                                        >  1:5073147/1‑150 (MQ=255)

                                                                   ccTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCg                                                                                    <  2:8551783/150‑1 (MQ=255)

                                                                    atgctgcCAACCACCCACACGCGGGTAGACCAAGCGGATGTCACGGGCCTCCAGCAACAACGCGCTCAAGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGc                                                                                   >  2:5123721/2‑150 (MQ=255)

                                                                           cAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                            <  2:1522759/150‑1 (MQ=255)

                                                                           cAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                            >  2:1144325/1‑150 (MQ=255)

                                                                           cAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                            <  1:3775501/150‑1 (MQ=255)

                                                                            aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                           <  2:3442151/150‑1 (MQ=255)

                                                                            aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                           >  1:7083519/1‑150 (MQ=255)

                                                                            aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                           >  1:1950565/1‑150 (MQ=255)

                                                                             accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                          <  2:3290087/150‑1 (MQ=255)

                                                                              ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                         >  1:4771926/1‑150 (MQ=255)

                                                                                 cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaa                                                                      >  1:4673608/1‑150 (MQ=255)

                                                                                  ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCACTTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                     >  2:7364193/1‑150 (MQ=255)

                                                                                   cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCGATGGCgggagtcg                                                                    <  1:6427990/150‑6 (MQ=255)

                                                                                      aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCt                                                                 >  2:3028530/1‑150 (MQ=255)

                                                                                        gcgGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTc                                                               >  2:8209198/1‑150 (MQ=255)

                                                                                           ggTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                                            >  1:2741791/1‑150 (MQ=255)

                                                                                           ggTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                                            >  2:6353313/1‑150 (MQ=255)

                                                                                                 ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTc                                                      <  1:7432676/150‑1 (MQ=255)

                                                                                                 ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTc                                                      >  2:5890452/1‑150 (MQ=255)

                                                                                                 ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACACCCGCCATCGGCGGGCAAAACCTTCGCCGTATTc                                                      >  1:7773114/1‑150 (MQ=255)

                                                                                                  cAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                                     >  2:400759/1‑150 (MQ=255)

                                                                                                  cAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                                     >  2:6065630/1‑150 (MQ=255)

                                                                                                   aaGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                                   <  1:3961264/87‑1 (MQ=255)

                                                                                                      cGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                                 <  2:6444352/150‑1 (MQ=255)

                                                                                                       ggATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGTAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                                >  2:6228583/1‑150 (MQ=255)

                                                                                                       ggATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                                <  1:1143736/150‑1 (MQ=255)

                                                                                                        gATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAgg                                               >  1:3284699/1‑150 (MQ=255)

                                                                                                        gATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAgg                                               >  1:6183086/1‑150 (MQ=255)

                                                                                                            tCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAg                                           <  1:577381/150‑1 (MQ=255)

                                                                                                             cACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                          >  1:5192448/1‑150 (MQ=255)

                                                                                                                 ggCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTg                                      <  1:2748666/150‑1 (MQ=255)

                                                                                                                  gCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgc                                     >  1:2105457/1‑150 (MQ=255)

                                                                                                                    cTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccg                                   >  2:2539709/1‑150 (MQ=255)

                                                                                                                     tCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCACCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                                  >  2:5562520/1‑150 (MQ=255)

                                                                                                                       cAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCa                                >  1:6791004/1‑150 (MQ=255)

                                                                                                                        aGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                               <  2:6637924/150‑1 (MQ=255)

                                                                                                                         gCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACt                              >  2:2880750/1‑150 (MQ=255)

                                                                                                                         gCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACt                              >  1:8540185/1‑150 (MQ=255)

                                                                                                                          cAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                             <  1:4931583/150‑1 (MQ=255)

                                                                                                                            acacCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCATAGGTAAGCCTTGCCGCCACTGGc                           >  2:8459186/1‑150 (MQ=255)

                                                                                                                             cacCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGcc                          <  2:4834063/150‑1 (MQ=255)

                                                                                                                                cgcgCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc                       >  1:6659841/1‑150 (MQ=255)

                                                                                                                                cgcgCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc                       >  1:5882084/1‑150 (MQ=255)

                                                                                                                                cgcgCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc                       >  1:8414950/1‑150 (MQ=255)

                                                                                                                                     tCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCCAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc                  >  2:7802898/1‑150 (MQ=255)

                                                                                                                                     tCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc                  >  2:3187617/1‑150 (MQ=255)

                                                                                                                                      cATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTcc                 <  2:4723045/150‑1 (MQ=255)

                                                                                                                                       aTGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCATAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa                >  1:310406/1‑150 (MQ=255)

                                                                                                                                       aTGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa                >  1:1566781/1‑150 (MQ=255)

                                                                                                                                        ttgcgcGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg               >  2:4148665/6‑150 (MQ=255)

                                                                                                                                        tGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg               <  2:2626372/150‑1 (MQ=255)

                                                                                                                                         gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGTTCGGCACGAAGGCCGACACCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt              <  1:4513773/150‑1 (MQ=255)

                                                                                                                                         gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt              <  2:1788465/150‑1 (MQ=255)

                                                                                                                                          cAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGtt             >  2:3850298/1‑150 (MQ=255)

                                                                                                                                           aaCGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTg            >  2:6983549/1‑150 (MQ=255)

                                                                                                                                              gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                         >  1:3606299/1‑134 (MQ=255)

                                                                                                                                                 tgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                          >  2:4252824/1‑82 (MQ=255)

                                                                                                                                                 tgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                          <  1:4252824/82‑1 (MQ=255)

                                                                                                                                                 tgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACcgc      >  2:4931936/1‑150 (MQ=255)

                                                                                                                                                    ctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGTc   <  1:6594130/150‑1 (MQ=255)

ACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGTCG  >  NC_002947/4135798‑4136096

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