BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	311,762	T→A	10.8%	S206C (AGC→TGC) ‡	PP_0256 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	311,762	0	T	A	10.8%	583.9 / 46.9	194	S206C (AGC→TGC) ‡	PP_0256	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (70/103); new base A (6/15); total (76/118)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 203 total aligned reads displayed.

CACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGT  >  NC_002947/311614‑311911

 acacCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCttt                                                                                                                                                     >  1:1042701/1‑150 (MQ=255)

  caccacCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTg                                                                                                                                                    >  1:3241060/1‑150 (MQ=255)

   accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                                                   <  1:5233092/150‑1 (MQ=255)

   accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                                                   <  2:2306805/150‑1 (MQ=255)

   accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                                                   <  2:2306846/150‑1 (MQ=255)

         tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                                                             <  2:63729/150‑1 (MQ=255)

         tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                                                             <  2:3829754/150‑1 (MQ=255)

              gcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACCAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCCacgccac                                                                                                                                        >  2:4342832/1‑150 (MQ=255)

              gcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccac                                                                                                                                        >  2:1879979/1‑150 (MQ=255)

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

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

                    cTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCt                                                                                                                                  <  2:3756026/150‑1 (MQ=255)

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

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

                     tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                                                                 >  1:420534/1‑150 (MQ=255)

                       tACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGc                                                                                                                               <  1:3753003/150‑1 (MQ=255)

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

                               cACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGc                                                                                                                       <  1:673100/150‑1 (MQ=255)

                                  gCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGcaca                                                                                                                    >  1:2395689/1‑150 (MQ=255)

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

                                        ggTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCg                                                                                                              <  2:3009440/150‑1 (MQ=255)

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

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

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

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

                                          cccaacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGaa                                                                                                            <  2:2088191/145‑2 (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)

                                             aGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                                                                  <  2:3589351/109‑1 (MQ=255)

                                              gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                                                        >  1:361604/1‑150 (MQ=255)

                                              gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                                                        <  2:3126731/150‑1 (MQ=255)

                                              gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                                                        <  2:3819626/150‑1 (MQ=255)

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

                                                 tGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAg                                                                                                     <  1:2803231/150‑1 (MQ=255)

                                                   cgcgGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAggg                                                                                                   <  2:2399080/150‑1 (MQ=255)

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

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

                                                     cgGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGaa                                                                                                 <  2:1102996/150‑1 (MQ=255)

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

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

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

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

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

                                                            ttGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttc                                                                                          <  1:2149738/150‑2 (MQ=255)

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

                                                             tGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttcg                                                                                         <  1:1847330/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)

                                                                  ccGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCg                                                                                    >  1:3681481/1‑150 (MQ=255)

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

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

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

                                                                         tCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGgcgc                                                                             <  1:2220821/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                          aTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTAc                                                            <  1:3163121/150‑1 (MQ=255)

                                                                                            ggCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAg                                                          <  1:2721429/150‑1 (MQ=255)

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

                                                                                                      tGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCt                                                <  1:2033693/150‑1 (MQ=255)

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

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

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

                                                                                                           gAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGtt                                           >  1:405256/1‑150 (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                                       >  2:1045586/1‑150 (MQ=255)

                                                                                                                 gtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc                                     >  1:4744636/1‑150 (MQ=255)

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

                                                                                                                       ggTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTcc                               >  2:4347772/1‑150 (MQ=255)

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

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

                                                                                                                          g‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGa                            <  1:150501/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                            <  2:654267/145‑1 (MQ=255)

                                                                                                                                    tttGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGaa                   <  1:3451392/150‑1 (MQ=255)

                                                                                                                                    tttGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGaa                   <  1:1826932/150‑1 (MQ=255)

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

                                                                                                                                         cGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCgg              >  1:2669334/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                                gACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTc       >  1:3921349/1‑150 (MQ=255)

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

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

CACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGGGT  >  NC_002947/311614‑311911

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