BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	311,734	G→C	8.9%	T215S (ACC→AGC) ‡	PP_0256 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	311,734	0	G	C	8.9%	651.0 / 41.8	201	T215S (ACC→AGC) ‡	PP_0256	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (92/91); new base C (13/5); total (105/96)
Fisher's exact test for biased strand distribution p-value = 8.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 203 total aligned reads displayed.

CCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCC  >  NC_002947/311585‑311882

 cAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGt                                                                                                                                                      >  1:2217523/1‑150 (MQ=255)

    aCCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑Ac                                                                                                                                                  <  2:2712177/150‑1 (MQ=255)

    aCCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑Ac                                                                                                                                                  >  1:4549779/1‑150 (MQ=255)

    aCCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑Ac                                                                                                                                                  <  1:4734001/150‑1 (MQ=255)

      cACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTg                                                                                                                                                <  1:4145270/150‑1 (MQ=255)

       aCGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTTGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGt                                                                                                                                               <  2:3221025/150‑1 (MQ=255)

       aCGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGt                                                                                                                                               <  2:1153929/150‑1 (MQ=255)

             cTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCttt                                                                                                                                         >  2:2957935/1‑150 (MQ=255)

               tgtgCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGc                                                                                                                                       >  2:510447/4‑150 (MQ=255)

               tGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGc                                                                                                                                       >  2:35785/1‑150 (MQ=255)

                gAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGcc                                                                                                                                      >  1:4493426/1‑150 (MQ=255)

                  gCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGa                                                                                                                                    >  1:4379783/1‑150 (MQ=255)

                  gCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGa                                                                                                                                    >  1:1801441/1‑150 (MQ=255)

                           tGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGc                                                                                                                           >  1:2049256/1‑150 (MQ=255)

                           tGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGc                                                                                                                           <  1:3951269/150‑1 (MQ=255)

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

                                accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                      >  1:3203779/1‑150 (MQ=255)

                                 ccaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                                     <  1:5116464/150‑1 (MQ=255)

                                     cTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCaa                                                                                                                 <  1:1573443/150‑1 (MQ=255)

                                     cTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCaa                                                                                                                 >  2:4441832/1‑150 (MQ=255)

                                     cTGGGCGCCACGCTGTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCaa                                                                                                                 <  2:714537/150‑1 (MQ=255)

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

                                      tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                                >  1:1371672/1‑150 (MQ=255)

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

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

                                             cACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacgc                                                                                                         >  1:3812465/1‑150 (MQ=255)

                                                gatttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCAcgccgc                                                                                                      >  1:75939/3‑150 (MQ=255)

                                                gCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCAcgccgc                                                                                                      <  1:1067670/150‑1 (MQ=255)

                                                gCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCAcgccgc                                                                                                      >  1:76235/1‑150 (MQ=255)

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

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

                                                      ccGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCt                                                                                                >  2:5004616/1‑150 (MQ=255)

                                                       cGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTc                                                                                               >  2:1890110/1‑150 (MQ=255)

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

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

                                                           tCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                           >  2:510246/1‑150 (MQ=255)

                                                             aCGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGca                                                                                         <  2:4912244/150‑1 (MQ=255)

                                                              cGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGcac                                                                                        <  1:3780863/150‑1 (MQ=255)

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

                                                                 aGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATg                                                                                     >  1:2630921/1‑150 (MQ=255)

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

                                                                        ccaacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGaac                                                                              >  1:4851817/5‑148 (MQ=255)

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

                                                                             aTGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCa                                                                         >  2:1201987/1‑150 (MQ=255)

                                                                             aTGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGTAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCAGATGTTCGCGCCGTCa                                                                         >  2:3705254/1‑150 (MQ=255)

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

                                                                              tGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCAGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAg                                                                        <  2:4535202/150‑1 (MQ=255)

                                                                                 ggggatGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGa                                                                     >  1:2024516/6‑150 (MQ=255)

                                                                                 gcgGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGa                                                                     <  1:35785/150‑1 (MQ=255)

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

                                                                                   gggTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAg                                                                   >  1:633397/1‑150 (MQ=255)

                                                                                     gTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgtt                                                                 <  2:4822750/150‑1 (MQ=255)

                                                                                        aTTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgttgtt                                                              <  2:3266216/150‑1 (MQ=255)

                                                                                        aTTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgttgtt                                                              >  2:891336/1‑150 (MQ=255)

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

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

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

                                                                                           ggTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGt                                                           >  2:4071296/1‑150 (MQ=255)

                                                                                             tGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGc                                                         >  1:1078636/1‑150 (MQ=255)

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

                                                                                                  gggTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAg                                                    <  1:5050399/150‑1 (MQ=255)

                                                                                                       cTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGgcgcg                                               <  2:2444557/150‑1 (MQ=255)

                                                                                                            cccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGa                                          >  2:5037213/1‑150 (MQ=255)

                                                                                                             ccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGaa                                         <  1:1179115/150‑1 (MQ=255)

                                                                                                             cacgACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAGCGCGCACGaa                                         >  2:1769160/4‑150 (MQ=255)

                                                                                                              cAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACTCCGACGCATTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAAc                                        >  2:4563007/1‑150 (MQ=255)

                                                                                                               aGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACa                                       <  1:1450130/150‑1 (MQ=255)

                                                                                                                gACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAGCGCGCACGAAAAg                                      >  2:1981830/1‑150 (MQ=255)

                                                                                                                gACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAGCGCGCACGAAAAg                                      >  1:2628648/1‑150 (MQ=255)

                                                                                                                gACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAg                                      <  1:1287589/150‑1 (MQ=255)

                                                                                                                   aaaaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTg                                   >  2:2550110/1‑150 (MQ=255)

                                                                                                                   aaaaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTg                                   >  2:1261873/1‑150 (MQ=255)

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

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

                                                                                                                     aaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGt                                 >  1:4055499/1‑150 (MQ=255)

                                                                                                                      aagggCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTa                                <  2:4304950/147‑1 (MQ=255)

                                                                                                                      aaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTa                                >  1:36465/1‑150 (MQ=255)

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

                                                                                                                          gCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAgg                            <  1:2125883/150‑1 (MQ=255)

                                                                                                                                  gTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTAcgccgc                    >  1:2688211/1‑150 (MQ=255)

                                                                                                                                   tGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCt                   >  2:839397/1‑150 (MQ=255)

                                                                                                                                    gCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCtt                  <  1:1838409/150‑1 (MQ=255)

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

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

                                                                                                                                       cGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGt               <  2:3945599/150‑1 (MQ=255)

                                                                                                                                       cGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGt               <  2:1942285/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                            tcgtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTGCGCCGCTTCGTTGCTg          <  2:578264/150‑1 (MQ=255)

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

                                                                                                                                              gtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc        <  2:3053388/150‑1 (MQ=255)

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

                                                                                                                                              gtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc        <  2:234054/150‑1 (MQ=255)

                                                                                                                                               tcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCa       >  2:4018821/1‑150 (MQ=255)

                                                                                                                                                cgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCAc      <  2:1426981/150‑1 (MQ=255)

                                                                                                                                                  aaGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGt    <  1:738258/150‑1 (MQ=255)

                                                                                                                                                   aGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTc   <  2:336186/150‑1 (MQ=255)

                                                                                                                                                    ggTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTcc  <  1:2533710/150‑1 (MQ=255)

CCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCC  >  NC_002947/311585‑311882

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