BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,793,940	A→T	12.5%	I209N (ATT→AAT)	PP_3355 ←	beta‑ketothiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,793,940	0	A	T	12.5%	323.6 / 32.9	112	I209N (ATT→AAT)	PP_3355	beta‑ketothiolase
Reads supporting (aligned to +/- strand): ref base A (46/52); new base T (7/7); total (53/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 119 total aligned reads displayed.

GCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCAT  >  NC_002947/3793793‑3794087

gccgccGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGcaatgcaatg                                                                                                                                                   <  2:2504894/150‑1 (MQ=255)

gccgccGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGcaatgcaatg                                                                                                                                                   <  1:1827343/150‑1 (MQ=255)

gccgccGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGcaatgcaatg                                                                                                                                                   <  1:1774469/150‑1 (MQ=255)

  cgccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGcc                                                                                                                                                 <  1:2192650/150‑1 (MQ=255)

   gccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCcg                                                                                                                                                >  2:1599430/1‑150 (MQ=255)

        gcCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGAACGCCACCAGCGTGC‑CCGGCATCTAATCGAGCCAGTGAGTCCGCGTCTGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATTCCGCGtgg                                                                                                                                           >  1:375833/1‑150 (MQ=255)

          cATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGtggct                                                                                                                                         >  2:2491025/1‑150 (MQ=255)

            tCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGtggctgg                                                                                                                                       >  2:282755/1‑150 (MQ=255)

              gACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCa                                                                                                                                     >  1:894396/1‑150 (MQ=255)

                 caccGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCcgcg                                                                                                                                              >  2:686535/1‑138 (MQ=255)

                 caccGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCcgcg                                                                                                                                              <  1:686535/138‑1 (MQ=255)

                  accGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGCCGCGtg                                                                                                                                            >  2:1323755/1‑139 (MQ=255)

                  accGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGCCGCGtg                                                                                                                                            <  1:1323755/139‑1 (MQ=255)

                  accGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCt                                                                                                                                 >  2:1688846/1‑150 (MQ=255)

                   ccGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGCCGCGtg                                                                                                                                            <  1:2696898/138‑1 (MQ=255)

                    cgcgCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGTCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTc                                                                                                                               <  2:1641927/150‑1 (MQ=255)

                    cgcgCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTc                                                                                                                               <  1:2523151/150‑1 (MQ=255)

                    cgcgCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTc                                                                                                                               <  2:2167656/150‑1 (MQ=255)

                    cgcgCAACTGTTGCCACCCGTCTGCACGCCACCGGCGTGCACCGGTCAC‑AGTCGGGACCGTTCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGTCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTc                                                                                                                               >  2:176878/1‑150 (MQ=255)

                       gcAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTAc                                                                                                                            <  2:1020930/150‑1 (MQ=255)

                         aaCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTAc                                                                                                                            <  2:2248017/148‑1 (MQ=255)

                         aaCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTAc                                                                                                                            >  1:2248017/1‑148 (MQ=255)

                           cTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTcc                                                                                                                        >  1:1394587/1‑150 (MQ=255)

                            tGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCa                                                                                                                       <  2:655324/150‑1 (MQ=255)

                              ttGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGt                                                                                                                     >  2:214917/1‑150 (MQ=255)

                                gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                                   <  1:1269280/150‑1 (MQ=255)

                                 ccAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGcc                                                                                                                  >  2:1743377/1‑150 (MQ=255)

                                  cAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTAGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCt                                                                                                                 >  2:1769276/1‑150 (MQ=255)

                                     ccGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTgcg                                                                                                              >  2:1038273/1‑150 (MQ=255)

                                          tGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTc                                                                                                                         <  2:2799514/134‑1 (MQ=255)

                                          tGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTc                                                                                                                         >  1:2799514/1‑134 (MQ=255)

                                             aCGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACg                                                                                                      <  2:1665005/150‑1 (MQ=255)

                                                ccACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTg                                                                                                   >  2:1470699/1‑150 (MQ=255)

                                                   ccgccgTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACg                                                                                                <  1:2820199/150‑1 (MQ=255)

                                                    cGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGAc                                                                                                       >  1:2053098/1‑142 (MQ=255)

                                                      gCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGAtt                                                                                             <  1:1688846/150‑1 (MQ=255)

                                                        gTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTc                                                                                           <  1:2752425/150‑1 (MQ=255)

                                                        gTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTc                                                                                           >  2:1378356/1‑150 (MQ=255)

                                                         tGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTcc                                                                                          <  1:1407388/150‑1 (MQ=255)

                                                         tGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTcc                                                                                          <  1:169011/150‑1 (MQ=255)

                                                         tGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTcc                                                                                          <  2:2944486/150‑1 (MQ=255)

                                                           cACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCtt                                                                                                                                >  2:279624/1‑110 (MQ=255)

                                                           cACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCtt                                                                                                                                <  1:279624/110‑1 (MQ=255)

                                                           cACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtc                                                                                        <  1:1524394/150‑1 (MQ=255)

                                                            aCCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcg                                                                                       >  1:1720184/1‑150 (MQ=255)

                                                            aCCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcg                                                                                       >  2:2761747/1‑150 (MQ=255)

                                                            aCCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTACCTAACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcg                                                                                       <  2:2652437/150‑1 (MQ=255)

                                                              cGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGGCGCTGACGATTTCCtcgtc                                                                                     <  2:720586/150‑1 (MQ=255)

                                                              cGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtc                                                                                     <  1:214917/150‑1 (MQ=255)

                                                                  cgc‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGGGTGGCTGGCAGCCTTCTCCCGCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAAc                                                                                 <  2:2411273/150‑1 (MQ=255)

                                                                   gc‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgt                                                                                      >  1:913810/1‑144 (MQ=255)

                                                                        tCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGcc                                                                            <  1:1769276/150‑1 (MQ=255)

                                                                        tCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGcc                                                                            >  2:1584979/1‑150 (MQ=255)

                                                                           ggCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTg                                                                         >  2:2520237/1‑150 (MQ=255)

                                                                                 tGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTg                                                                   >  1:1031902/1‑150 (MQ=255)

                                                                                 tGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTACCCTCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGCCGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTg                                                                   <  1:2491025/150‑1 (MQ=255)

                                                                                    ggCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGc                                                                <  1:682855/150‑1 (MQ=255)

                                                                                         cgTCGGTGACCCTTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCa                                                           <  1:2433033/150‑1 (MQ=255)

                                                                                             ggTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGc                                                       >  2:2964125/1‑150 (MQ=255)

                                                                                              gTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGcc                                                      >  2:258755/1‑150 (MQ=255)

                                                                                                  gCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCt                                                  <  2:2722567/150‑1 (MQ=255)

                                                                                                     gtgtgtggCTGTCCTGCCTGACCTCACTCACGCCACGCGCCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTCACGATTTCCTCCTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGAt                                               <  1:176878/150‑1 (MQ=255)

                                                                                                     gtgggtggCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgt                                                                                      <  2:1542789/111‑1 (MQ=255)

                                                                                                     gtgggtggCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgt                                                                                      >  1:1542789/1‑111 (MQ=255)

                                                                                                             cTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTgcgc                                       >  2:1387654/1‑150 (MQ=255)

                                                                                                               gTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAg                                     <  1:1038273/150‑1 (MQ=255)

                                                                                                                   tGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTAATGGCTGCGCAGGGcc                                 <  1:1107950/150‑1 (MQ=255)

                                                                                                                     gCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAgg                                                         >  2:1726150/1‑124 (MQ=255)

                                                                                                                     gCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAgg                                                         <  1:1726150/124‑1 (MQ=255)

                                                                                                                     gCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCAGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCa                               >  1:1700029/1‑150 (MQ=255)

                                                                                                                       tGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAg                                                                              <  2:2741367/101‑1 (MQ=255)

                                                                                                                           ctcactcaCGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCa                                                           >  2:391175/1‑116 (MQ=255)

                                                                                                                           ctcactcaCGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCa                                                           <  1:391175/116‑1 (MQ=255)

                                                                                                                           ctcactcaCGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTc                         <  2:1031902/150‑1 (MQ=255)

                                                                                                                            tcactcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTcc                        >  2:761309/1‑150 (MQ=255)

                                                                                                                               ctcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCAcc                     >  1:2508396/1‑150 (MQ=255)

                                                                                                                                tcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCAcgc                            >  2:286534/1‑142 (MQ=255)

                                                                                                                                tcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCAcgc                            <  1:286534/142‑1 (MQ=255)

                                                                                                                                tcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCg                    >  1:2766248/1‑150 (MQ=255)

                                                                                                                                     ccacgcGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTCCCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCa               <  2:375833/150‑1 (MQ=255)

                                                                                                                                     ccacgcGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCa               <  1:2897614/150‑1 (MQ=255)

                                                                                                                                      cacgcGGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAgg                                    <  2:227668/128‑1 (MQ=255)

                                                                                                                                      cacgcGGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAgg                                    >  1:227668/1‑128 (MQ=255)

                                                                                                                                      cacgcGGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt                  <  2:13883/146‑1 (MQ=255)

                                                                                                                                      cacgcGGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt                  >  1:13883/1‑146 (MQ=255)

                                                                                                                                      cacgcGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAc              >  2:1309901/1‑150 (MQ=255)

                                                                                                                                       acgcGGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTgg                                                                  <  2:2916585/97‑1 (MQ=255)

                                                                                                                                       acgcGGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTgg                                                                  >  1:2916585/1‑97 (MQ=255)

                                                                                                                                       acgcGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCAcc                     <  1:2705924/142‑1 (MQ=255)

                                                                                                                                       acgcGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCAcc                     >  2:2705924/1‑142 (MQ=255)

                                                                                                                                        cgcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcgc            >  2:912415/1‑150 (MQ=255)

                                                                                                                                         gcgGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATg                                              >  2:3046793/1‑115 (MQ=255)

                                                                                                                                         gcgGCATTGCATTGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATg                                              <  1:3046793/115‑1 (MQ=255)

                                                                                                                                         gcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGACCCAGCCCTGCACCTGGGCTTCCAGGGCCCGCTGATGGCTGCGCACGGCCCAAGCTTCCCCAGTTCCAcccg           >  2:490808/1‑150 (MQ=255)

                                                                                                                                           ggCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGcc         >  1:674169/1‑150 (MQ=255)

                                                                                                                                           ggCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGcc         >  1:134516/1‑150 (MQ=255)

                                                                                                                                                tgcaatgcCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGcc    >  1:792753/1‑150 (MQ=255)

                                                                                                                                                 gcaatgcCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTgc                                         <  1:49724/112‑1 (MQ=255)

                                                                                                                                                 gcaatgcCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTgc                                         >  2:49724/1‑112 (MQ=255)

                                                                                                                                                  caatgcCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCAt  >  1:2459212/1‑150 (MQ=255)

GCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGC‑AGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCAT  >  NC_002947/3793793‑3794087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`