BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,595,071	A→G	6.0%	Q434R (CAA→CGA)	PP_2270 →	DNA primase/helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,595,071	0	A	G	6.0%	420.9 / 10.0	133	Q434R (CAA→CGA)	PP_2270	DNA primase/helicase
Reads supporting (aligned to +/- strand): ref base A (67/58); new base G (2/6); total (69/64)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 138 total aligned reads displayed.

TCACCTGACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGAATTCT  >  NC_002947/2594922‑2595220

tcACCTGACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCa                                                                                                                                                       >  2:4801137/1‑150 (MQ=255)

  aCCTGACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCaaa                                                                                                                                                     >  2:2408517/1‑150 (MQ=255)

  aCCTGACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCaaa                                                                                                                                                     >  2:2408543/1‑150 (MQ=255)

   ccTGACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAac                                                                                                                                                    >  1:557969/1‑150 (MQ=255)

      gACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccac                                                                                                                                                 >  1:1104106/1‑150 (MQ=255)

      gACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccac                                                                                                                                                 >  1:1104134/1‑150 (MQ=255)

      gACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccac                                                                                                                                                 <  1:103589/150‑1 (MQ=255)

      gACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccac                                                                                                                                                 >  1:4568626/1‑150 (MQ=255)

        ccaccCCCGAGGGCCACCCCCCTGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAACGGCTCCCGAGCCATCGGCTTCGGGAGCGACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTTGGCGAaccacgg                                                                                                                                               <  2:4025214/150‑3 (MQ=255)

        ccaccCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTTGGCGAaccacgg                                                                                                                                               <  2:4413374/150‑3 (MQ=255)

        ccaccCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccacca                                                                                                                                               <  2:3774640/150‑1 (MQ=255)

         caccCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCACAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTTGGCGAaccacggg                                                                                                                                              <  2:2771239/150‑4 (MQ=255)

         caccCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTTGGCGAaccacggg                                                                                                                                              <  2:3041784/150‑4 (MQ=255)

         caccCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccaccac                                                                                                                                              >  1:72189/1‑150 (MQ=255)

         caccCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAaccaccac                                                                                                                                              >  2:3248050/1‑150 (MQ=255)

           cccccGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCt                                                                                                                                            >  1:3620120/1‑150 (MQ=255)

            ccccGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCtt                                                                                                                                           >  2:713876/1‑150 (MQ=255)

              ccGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCAGGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTGcc                                                                                                                                         >  2:4625979/1‑150 (MQ=255)

                gAGGGCAAACCCCATGAGGAAGGCGGGTGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGc                                                                                                                                       <  2:1684559/150‑1 (MQ=255)

                      aaaCCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTc                                                                                                                                 >  1:1788560/1‑150 (MQ=255)

                        aCCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCaa                                                                                                                               <  1:3554446/150‑1 (MQ=255)

                        aCCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCaa                                                                                                                               <  2:2164774/150‑1 (MQ=255)

                           ccATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGa                                                                                                                            >  1:4761288/1‑150 (MQ=255)

                           ccATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGa                                                                                                                            >  1:23930/1‑150 (MQ=255)

                            cATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGAc                                                                                                                           >  2:4011678/1‑150 (MQ=255)

                             aTGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGAcc                                                                                                                          >  2:1434187/1‑150 (MQ=255)

                                      gcgggcggGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCg                                                                                                                 <  2:1815562/150‑1 (MQ=255)

                                      gcgggcggGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCg                                                                                                                 <  2:1815678/150‑1 (MQ=255)

                                        gggcgggTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGc                                                                                                               <  2:23930/150‑1 (MQ=255)

                                           cgggTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCaa                                                                                                            <  2:4330072/150‑1 (MQ=255)

                                              gTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCg                                                                                                         >  1:659480/1‑150 (MQ=255)

                                              gTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCg                                                                                                         >  1:659564/1‑150 (MQ=255)

                                                   gATTCGCCACTTCAACGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTGCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCgg                                                                                                    <  2:1718465/150‑1 (MQ=255)

                                                   gATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCgg                                                                                                    <  1:5109429/150‑1 (MQ=255)

                                                   gATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCgg                                                                                                    <  1:2081915/150‑1 (MQ=255)

                                                       cGCCAATTCAAAGGCTCCCGAGCCATCGGATTCTGGAGCTACTTCTTGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCACGGACCGCTACACCGGCCAAGCGACCGGGGaa                                                                                                <  1:1656305/150‑1 (MQ=255)

                                                        gCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAg                                                                                               >  1:2828035/1‑150 (MQ=255)

                                                         ccACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGt                                                                                              >  2:3413629/1‑150 (MQ=255)

                                                         ccACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGt                                                                                              >  1:2288215/1‑150 (MQ=255)

                                                         ccACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGt                                                                                              >  1:3270237/1‑150 (MQ=255)

                                                          cACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTg                                                                                             >  1:2659867/1‑150 (MQ=255)

                                                           aCTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGc                                                                                            <  1:3093842/150‑1 (MQ=255)

                                                            cTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCGGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCt                                                                                           <  1:2843822/150‑1 (MQ=255)

                                                            cTTCAAAGGCTCCCGAGCCATCGGCGTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCt                                                                                           <  2:5008862/150‑1 (MQ=255)

                                                               cAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTa                                                                                        <  2:4160470/150‑1 (MQ=255)

                                                                aaaGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTAc                                                                                       >  2:1802020/1‑150 (MQ=255)

                                                                aaaGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTAc                                                                                       >  1:3742590/1‑150 (MQ=255)

                                                                aaaGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTAc                                                                                       >  1:4202214/1‑150 (MQ=255)

                                                                  aGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGaa                                                                                                >  2:4195801/1‑139 (MQ=255)

                                                                   ggCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTACGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTg                                                                                    <  1:1521471/150‑1 (MQ=255)

                                                                       cccGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCAt                                                                                <  2:4642435/150‑1 (MQ=255)

                                                                        ccGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCAGTACCTGGCATa                                                                               >  2:4194092/1‑150 (MQ=255)

                                                                               aTCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCACCAGGCCGAAGACCCCGTGGTTCGCAAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACcgc                                                                        <  1:31006/150‑1 (MQ=255)

                                                                               aTCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACcgc                                                                        <  2:5061892/150‑1 (MQ=255)

                                                                                     ttCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCggtgta                                                                                                >  1:2260501/1‑116 (MQ=255)

                                                                                     ttCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACAcc                                                                  >  2:3161586/1‑150 (MQ=255)

                                                                                       cTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCgg                                                                >  2:4255193/1‑150 (MQ=255)

                                                                                        tGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGa                                                               >  1:4107579/1‑150 (MQ=255)

                                                                                          gAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGAtt                                                             <  1:3759192/150‑1 (MQ=255)

                                                                                          gAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGAtt                                                             >  2:4754544/1‑150 (MQ=255)

                                                                                            gCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGc                                                           >  2:3451502/1‑150 (MQ=255)

                                                                                                cTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCtctc                                                       >  1:3365055/1‑150 (MQ=255)

                                                                                                cTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCtctc                                                       >  2:1254995/1‑150 (MQ=255)

                                                                                                cTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCtctc                                                       >  2:1255114/1‑150 (MQ=255)

                                                                                                cTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCtctc                                                       <  1:392212/150‑1 (MQ=255)

                                                                                                 ttCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCg                                                      >  1:667554/1‑150 (MQ=255)

                                                                                                 ttCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCg                                                      >  2:322261/1‑150 (MQ=255)

                                                                                                 ttCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCg                                                      >  2:3711103/1‑150 (MQ=255)

                                                                                                         cGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGa                                              >  2:4138114/1‑150 (MQ=255)

                                                                                                           gCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAg                                            >  2:727133/1‑150 (MQ=255)

                                                                                                            ccTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAgc                                           <  2:3035815/150‑1 (MQ=255)

                                                                                                            ccTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAgc                                           <  2:3086196/150‑1 (MQ=255)

                                                                                                                gAGCGTGACCAACATGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAACGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGccc                                       <  2:743230/150‑1 (MQ=255)

                                                                                                                   cGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAg                                    <  2:1459021/150‑1 (MQ=255)

                                                                                                                   cGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAg                                    <  2:1459016/150‑1 (MQ=255)

                                                                                                                      gACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCg                                 <  1:697101/150‑1 (MQ=255)

                                                                                                                         cAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTcc                              >  2:77870/1‑150 (MQ=255)

                                                                                                                          aaCAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCt                             >  2:1865429/1‑150 (MQ=255)

                                                                                                                          aaCAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCt                             <  1:3122678/150‑1 (MQ=255)

                                                                                                                             aGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCAc                          >  2:3641372/1‑150 (MQ=255)

                                                                                                                               gCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCAccc                        <  2:3356781/150‑1 (MQ=255)

                                                                                                                               gCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCAccc                        <  1:713876/150‑1 (MQ=255)

                                                                                                                                ccGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCt                       <  1:4801137/150‑1 (MQ=255)

                                                                                                                                 cGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCg                                                      >  1:40127/1‑118 (MQ=255)

                                                                                                                                  gAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGACGCGCCCGAGCCGACCTCACCCtta                     >  1:2739679/1‑149 (MQ=255)

                                                                                                                                  gAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCttt                     >  1:1517206/1‑150 (MQ=255)

                                                                                                                                     gACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTaaa                  <  1:4278324/150‑1 (MQ=255)

                                                                                                                                        cccGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGAt               <  1:3248050/150‑1 (MQ=255)

                                                                                                                                         tggTGGTTTGGCGAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATg              >  2:1919775/3‑150 (MQ=255)

                                                                                                                                         tggTGGTTTGGCGAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATg              >  2:1920065/3‑150 (MQ=255)

                                                                                                                                         ccGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATg              <  2:4616315/150‑1 (MQ=255)

                                                                                                                                          cGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGa             <  2:3620120/150‑1 (MQ=255)

                                                                                                                                              gTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATcc         <  2:2828026/150‑1 (MQ=255)

                                                                                                                                              gTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATcc         <  2:2828035/150‑1 (MQ=255)

                                                                                                                                                tCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGa       <  2:2132597/150‑1 (MQ=255)

                                                                                                                                                 cGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGaa      <  2:557969/150‑1 (MQ=255)

                                                                                                                                                  gCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGAAt     >  2:472163/1‑150 (MQ=255)

                                                                                                                                                  gCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGAAt     <  1:4254082/150‑1 (MQ=255)

                                                                                                                                                  gCCAAACCACCACATTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGAAt     >  2:3061091/1‑150 (MQ=255)

                                                                                                                                                     aaaCCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGAATTCt  >  1:2303764/1‑150 (MQ=255)

TCACCTGACCACCCCCGAGGGCAAACCCCATGAGGAAGGCGGGCGGGTAACGATTCGCCACTTCAAAGGCTCCCGAGCCATCGGCTTCTGGAGCTACTTCATGTTCGGCCTTGAGCGTGACCAACAGGCCGAAGACCCCGTGGTTCGCCAAACCACCACCTTCCGCATCCTCAAGGACCGCTACACCGGCCAAGCGACCGGGGAAGTGCTGTACCTGGCATACGACCGCGACACCGGATTGCTCTCGCTCACCGAAGCGCCCGAGCCGTCCTCACCCTTTAAAGATGAATCCGAATTCT  >  NC_002947/2594922‑2595220

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