BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	311,778	G→A	10.1%	G200G (GGC→GGT)	PP_0256 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	311,778	0	G	A	10.1%	555.8 / 39.3	198	G200G (GGC→GGT)	PP_0256	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (74/103); new base A (5/15); total (80/118)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

CCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTGTTCAG  >  NC_002947/311629‑311927

ccACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacg                                                                                                                                                       >  2:830706/1‑150 (MQ=255)

  aCGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacgcc                                                                                                                                                     >  1:3475558/1‑150 (MQ=255)

      tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                                                                                 <  2:1284917/150‑1 (MQ=255)

      tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                                                                                 <  1:5181672/150‑1 (MQ=255)

        tACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTAGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGc                                                                                                                                               <  2:3829389/150‑1 (MQ=255)

         aCCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGcc                                                                                                                                              >  1:1766072/1‑150 (MQ=255)

            gCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCa                                                                                                                                           >  2:2858569/1‑150 (MQ=255)

             cATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCAt                                                                                                                                          >  1:284610/1‑150 (MQ=255)

                      gTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGt                                                                                                                                 >  1:412278/1‑150 (MQ=255)

                        gggTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                                                                               <  2:3603045/150‑1 (MQ=255)

                          gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                                                                             <  1:1185687/150‑1 (MQ=255)

                          gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                                                                             <  2:2790975/150‑1 (MQ=255)

                           cccaacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGaa                                                                                                                            <  2:4885896/145‑2 (MQ=255)

                             gAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCa                                                                                                                          >  1:4515812/1‑150 (MQ=255)

                             gAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCa                                                                                                                          <  1:5366498/150‑1 (MQ=255)

                              aacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAg                                                                                                                         >  2:3384367/3‑150 (MQ=255)

                                 aTGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCa                                                                                                                      <  2:1680453/150‑1 (MQ=255)

                                     gcgGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGa                                                                                                                  >  1:5320372/1‑150 (MQ=255)

                                     gcgGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGa                                                                                                                  >  1:2594523/1‑150 (MQ=255)

                                       ggatGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGaaa                                                                                                                <  1:2939829/147‑2 (MQ=255)

                                        ggTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                                                                        <  1:2296565/125‑1 (MQ=255)

                                        ggTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgt                                                                                                               >  2:4015479/1‑150 (MQ=255)

                                         atGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAatt                                                                                                              <  1:2798487/149‑1 (MQ=255)

                                         atGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAatt                                                                                                              <  2:1112865/149‑1 (MQ=255)

                                           gATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgttgt                                                                                                            <  2:1913405/150‑1 (MQ=255)

                                            aTTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgtt                                                                                                           <  2:2295193/150‑1 (MQ=255)

                                             ttGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGttgttg                                                                                                          <  1:379824/150‑1 (MQ=255)

                                              tGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                                                                               <  2:187281/128‑1 (MQ=255)

                                              tGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                                                                               >  1:187281/1‑128 (MQ=255)

                                              tGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttcg                                                                                                         <  1:1847329/150‑3 (MQ=255)

                                               ggTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTACTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGt                                                                                                        <  1:541350/150‑1 (MQ=255)

                                               ggTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGt                                                                                                        <  2:1789431/150‑1 (MQ=255)

                                                gTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTg                                                                                                       >  1:3339559/1‑150 (MQ=255)

                                                  gCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGcc                                                                                                     <  2:5064580/150‑1 (MQ=255)

                                                    cGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGt                                                                                                   <  1:2676809/150‑1 (MQ=255)

                                                      gggTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAg                                                                                                 >  1:1333434/1‑150 (MQ=255)

                                                      gggTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAg                                                                                                 >  1:1096766/1‑150 (MQ=255)

                                                       ggTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAgg                                                                                                >  2:3872207/1‑150 (MQ=255)

                                                            tGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGTCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGgcgcgc                                                                                           >  2:3096652/1‑150 (MQ=255)

                                                             ggCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCa                                                                                          >  1:3503488/1‑150 (MQ=255)

                                                              gCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCAc                                                                                         <  2:5335527/150‑1 (MQ=255)

                                                              gCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCAc                                                                                         <  1:694469/150‑1 (MQ=255)

                                                               ccccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCAcgccgc                                                                                                                                                   <  2:3796016/91‑1 (MQ=255)

                                                                cccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCCAAGGTG‑ACTGTACATTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGACTCATGGCACATGATCGAGCAGTCAGGGAAGATGATGGAGCCGTAGGCGCGCACGa                                                                                       >  1:436899/1‑150 (MQ=255)

                                                                ccacgACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAGCGCGCACGa                                                                                       >  1:4986358/5‑150 (MQ=255)

                                                                 ccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGaa                                                                                      >  2:4420022/1‑150 (MQ=255)

                                                                  cAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAAc                                                                                     <  2:1494475/150‑1 (MQ=255)

                                                                     aCAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGt                                                                                  >  1:2511125/1‑150 (MQ=255)

                                                                      cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                                                                                 >  2:1315288/1‑150 (MQ=255)

                                                                      cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                                                                                 <  2:2308561/150‑1 (MQ=255)

                                                                        aaaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTgg                                                                               <  2:1092420/150‑1 (MQ=255)

                                                                         aaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGt                                                                              <  1:4886850/150‑1 (MQ=255)

                                                                          aaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTa                                                                             <  1:3178297/150‑1 (MQ=255)

                                                                            tGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACa                                                                           <  2:2120215/150‑1 (MQ=255)

                                                                               cGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGt                                                                        <  1:4568972/150‑1 (MQ=255)

                                                                                gTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTa                                                                       >  1:4706220/1‑150 (MQ=255)

                                                                                     cGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTAcgccg                                                                  >  1:484852/1‑150 (MQ=255)

                                                                                        tcTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCtt                                                               <  2:2716548/149‑1 (MQ=255)

                                                                                         cTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTc                                                              <  2:3457583/150‑1 (MQ=255)

                                                                                          tCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCg                                                             >  2:1109370/1‑150 (MQ=255)

                                                                                            gAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGtt                                                           >  1:405256/1‑150 (MQ=255)

                                                                                            gAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGtt                                                           >  1:2081504/1‑150 (MQ=255)

                                                                                              aGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGc                                                         <  1:1265063/150‑1 (MQ=255)

                                                                                               gtcgtcGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCt                                                        <  2:1208242/150‑1 (MQ=255)

                                                                                               gtcgtcGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCt                                                        <  1:4420022/150‑1 (MQ=255)

                                                                                                tcgtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTcacccagc                                                       <  1:4265612/150‑8 (MQ=255)

                                                                                                tcgtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTg                                                       <  1:3271511/150‑1 (MQ=255)

                                                                                                   tcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCa                                                    <  1:3993225/150‑1 (MQ=255)

                                                                                                         atG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTccc                                              <  1:3788868/149‑1 (MQ=255)

                                                                                                           g‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGa                                            <  2:1714586/150‑1 (MQ=255)

                                                                                                           g‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGa                                            <  1:150501/150‑1 (MQ=255)

                                                                                                             aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                                           <  1:2170066/150‑1 (MQ=255)

                                                                                                             aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                                           <  1:2031952/150‑1 (MQ=255)

                                                                                                             aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                                           <  2:484852/150‑1 (MQ=255)

                                                                                                             aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                                           <  2:3212009/150‑1 (MQ=255)

                                                                                                               tGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGt                                         >  2:723154/1‑150 (MQ=255)

                                                                                                                 tACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGa                                            >  1:654267/1‑145 (MQ=255)

                                                                                                                  aCCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGtata                                      >  2:647013/1‑150 (MQ=255)

                                                                                                                         ccGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCg                               >  1:226955/1‑150 (MQ=255)

                                                                                                                           gACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGc                             <  1:537889/150‑1 (MQ=255)

                                                                                                                               ccGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc                                                     >  1:4212444/1‑122 (MQ=255)

                                                                                                                               ccGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc                                                     <  2:4212444/122‑1 (MQ=255)

                                                                                                                               ccGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTgg                         >  1:5357320/1‑150 (MQ=255)

                                                                                                                                cGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGt                        >  2:3338493/1‑150 (MQ=255)

                                                                                                                                    gCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAgg                    <  1:830706/150‑1 (MQ=255)

                                                                                                                                    gCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAgg                    <  1:4343128/150‑1 (MQ=255)

                                                                                                                                      tttGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGt                  >  2:4384395/1‑150 (MQ=255)

                                                                                                                                        tGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTg                <  1:3384574/150‑1 (MQ=255)

                                                                                                                                        tGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGAAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTg                <  1:1879979/150‑1 (MQ=255)

                                                                                                                                            ccAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgtt                                                                                                              <  2:3640003/52‑1 (MQ=255)

                                                                                                                                            ccAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgtt                                                                                                              >  1:3640003/1‑52 (MQ=255)

                                                                                                                                            ccAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCgg                              >  2:4547410/1‑132 (MQ=255)

                                                                                                                                            ccAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCa            <  2:1338981/150‑1 (MQ=255)

                                                                                                                                             cAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAg           >  1:3156222/1‑150 (MQ=255)

                                                                                                                                              aaCGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGc          >  1:3755972/1‑150 (MQ=255)

                                                                                                                                                 gccacgccGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTg       >  1:3741535/1‑150 (MQ=255)

                                                                                                                                                   cacgccGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTGtt     <  2:2653133/150‑1 (MQ=255)

                                                                                                                                                     cgccgcTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTGTTCa   <  2:5357177/150‑1 (MQ=255)

                                                                                                                                                     cgccgcTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTGTTCa   <  1:3687036/150‑1 (MQ=255)

                                                                                                                                                      gccgcTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTGTTCAg  <  1:2399458/150‑1 (MQ=255)

CCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGTTGGCCAGCTTGTTCAG  >  NC_002947/311629‑311927

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