BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,546,084	G→A	5.6%	H18H (CAC→CAT)	yjiC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,546,084	0	G	A	5.6%	1081.9 / 30.5	303	H18H (CAC→CAT)	yjiC	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (127/158); new base A (9/8); total (136/167)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 322 total aligned reads displayed.

AAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTC  >  CP009273/4545937‑4546233

aaGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTg                                                                                                                                                     >  1:298084/1‑150 (MQ=255)

aaGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTg                                                                                                                                                     >  2:2111329/1‑150 (MQ=255)

 aGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGa                                                                                                                                                    <  2:3253808/150‑1 (MQ=255)

       cTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTCCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCt                                                                                                                                              <  2:4596430/150‑1 (MQ=255)

         aaTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGt                                                                                                                                            >  2:4466006/1‑150 (MQ=255)

         aaTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGt                                                                                                                                            >  2:1022620/1‑150 (MQ=255)

          aTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGta                                                                                                                                           <  1:2432093/150‑1 (MQ=255)

          aTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGta                                                                                                                                           >  1:1508155/1‑150 (MQ=255)

               gACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGCGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTaataat                                                                                                                                      >  2:1607816/1‑150 (MQ=255)

                 cAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAAttt                                                                                                                                    >  1:512212/1‑150 (MQ=255)

                  aaCGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTg                                                                                                                                   >  1:31654/1‑150 (MQ=255)

                  aaCGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTg                                                                                                                                   >  1:3419828/1‑150 (MQ=255)

                       tttttCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGActct                                                                                                                              >  2:5061014/1‑150 (MQ=255)

                        ttttCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTg                                                                                                                             <  2:1712251/150‑1 (MQ=255)

                         tttCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTAGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGa                                                                                                                            >  1:1072820/1‑150 (MQ=255)

                           tCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAg                                                                                                                          <  2:1289797/150‑1 (MQ=255)

                            ccATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAgg                                                                                                                         >  2:654562/1‑150 (MQ=255)

                                cGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                                                     >  1:4653036/1‑150 (MQ=255)

                                cGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                                                     <  2:756532/150‑1 (MQ=255)

                                  tAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTg                                                                                                                   <  1:926237/150‑1 (MQ=255)

                                  tAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTg                                                                                                                   <  1:2910885/150‑1 (MQ=255)

                                  tAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTg                                                                                                                   <  2:529514/150‑1 (MQ=255)

                                  tAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTg                                                                                                                   >  1:1349149/1‑150 (MQ=255)

                                  tAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTg                                                                                                                   >  1:3802418/1‑150 (MQ=255)

                                        tGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                                             >  2:3491256/1‑150 (MQ=255)

                                         ggTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATg                                                                                                            >  1:2259522/1‑150 (MQ=255)

                                          gTTAGTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATgg                                                                                                           <  2:4205894/150‑1 (MQ=255)

                                           ttATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGc                                                                                                          >  1:862582/1‑150 (MQ=255)

                                           ttATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGc                                                                                                          >  1:4964046/1‑150 (MQ=255)

                                            tATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCa                                                                                                         >  1:1402375/1‑150 (MQ=255)

                                             aTTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCAt                                                                                                        >  2:5092306/1‑150 (MQ=255)

                                             aTTGCCACGAGTAAGAGGCGGAAGATAGATGTTTCCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCAt                                                                                                        >  1:1147731/1‑150 (MQ=255)

                                                 ccACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGAc                                                                                                    <  1:2101623/150‑1 (MQ=255)

                                                  cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACa                                                                                                   <  2:2745808/150‑1 (MQ=255)

                                                   aCGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACat                                                                                                  <  2:408348/150‑1 (MQ=255)

                                                    cGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACata                                                                                                 <  2:1261073/150‑1 (MQ=255)

                                                      aGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATACTTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataat                                                                                               >  2:2440074/1‑150 (MQ=255)

                                                        tAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAg                                                                                             >  2:3396880/1‑150 (MQ=255)

                                                        tAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAg                                                                                             >  1:3278091/1‑150 (MQ=255)

                                                            agGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTAc                                                                                         >  1:4575073/1‑150 (MQ=255)

                                                            agGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTAc                                                                                         <  2:4245166/150‑1 (MQ=255)

                                                               cGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGTAGGGCATTGCTTAATGGCATAGACATAATAGTTACCtt                                                                                      >  2:4628729/1‑150 (MQ=255)

                                                               cGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCtt                                                                                      <  2:4531792/150‑1 (MQ=255)

                                                               cGGAAGATAGATGTTTTCCGAATCCAAGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCtt                                                                                      <  2:3325903/150‑1 (MQ=255)

                                                                ggAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCttt                                                                                     <  1:3649904/150‑1 (MQ=255)

                                                                ggAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCttt                                                                                     <  1:2335239/150‑1 (MQ=255)

                                                                 gAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCtttt                                                                                    <  2:3780018/150‑1 (MQ=255)

                                                                   agatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTAt                                                                                  >  1:155363/1‑150 (MQ=255)

                                                                   agatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTAt                                                                                  <  1:4472140/150‑1 (MQ=255)

                                                                   agatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTAt                                                                                  >  2:3838177/1‑150 (MQ=255)

                                                                    gatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATg                                                                                 >  2:4840132/1‑150 (MQ=255)

                                                                    gatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTAACATTTTATg                                                                                 >  2:2837818/1‑150 (MQ=255)

                                                                     atagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGa                                                                                >  1:4988510/1‑150 (MQ=255)

                                                                     atagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGa                                                                                <  2:1971910/150‑1 (MQ=255)

                                                                     atagatGTTTTCAGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTAAAGGGAGATGAAGTAAGTGAGCATCTGTAATAATTTGACTATGAAGTGCATTGCTTAATGGCATAGACATACTATTTACCTTATATGa                                                                                >  2:3937134/1‑150 (MQ=255)

                                                                      tagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAg                                                                               >  1:2671550/1‑150 (MQ=255)

                                                                         atGTTTTCCGAATCCATGCGTGAATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAt                                                                            <  2:4024891/150‑1 (MQ=255)

                                                                           gTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAttt                                                                          <  2:565974/150‑1 (MQ=255)

                                                                            ttttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTc                                                                         <  2:2182007/150‑1 (MQ=255)

                                                                             tttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCa                                                                        <  1:4566347/150‑1 (MQ=255)

                                                                             tgtCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCa                                                                        <  1:4172922/148‑1 (MQ=255)

                                                                              ttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCAc                                                                       >  2:292926/1‑150 (MQ=255)

                                                                               tCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACt                                                                      <  1:2414465/150‑1 (MQ=255)

                                                                                  gAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGAt                                                                   >  1:2747810/1‑150 (MQ=255)

                                                                                    aTCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGt                                                                 <  1:659561/150‑1 (MQ=255)

                                                                                    aTCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGt                                                                 >  1:3533941/1‑150 (MQ=255)

                                                                                      ccATGCGTGCATTTTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGttt                                                               <  2:4936672/150‑1 (MQ=255)

                                                                                       cATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTa                                                              <  1:1901419/150‑1 (MQ=255)

                                                                                       cATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTAGTGGGTGATGAAGGAAGTGAGCATCTGTAAGAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTa                                                              >  2:3597248/1‑150 (MQ=255)

                                                                                             tGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGaaaaa                                                        <  1:1105912/150‑1 (MQ=255)

                                                                                                    tttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaa                                                 <  2:114180/150‑1 (MQ=255)

                                                                                                    tttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAAAAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaa                                                 >  2:3647631/1‑150 (MQ=255)

                                                                                                     ttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaaa                                                <  1:3967418/150‑1 (MQ=255)

                                                                                                     ttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaaa                                                <  1:3875817/150‑1 (MQ=255)

                                                                                                     ttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaaa                                                <  1:4296050/150‑1 (MQ=255)

                                                                                                       aaTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAtt                                              <  1:1634076/150‑1 (MQ=255)

                                                                                                       aaTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAtt                                              >  1:1728192/1‑150 (MQ=255)

                                                                                                       aaTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAtt                                              <  2:1311853/150‑1 (MQ=255)

                                                                                                        aTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAttt                                             <  1:371360/150‑1 (MQ=255)

                                                                                                          ttACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAttttt                                           <  1:1289115/150‑1 (MQ=255)

                                                                                                               tgtAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgta                                      <  1:4205250/150‑1 (MQ=255)

                                                                                                               tgtAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgta                                      <  1:2396402/150‑1 (MQ=255)

                                                                                                                gtAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATACATTTTTCTgtag                                     >  2:261470/1‑150 (MQ=255)

                                                                                                                gtAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAGTTGACTCTGAAGGGCATTGCTTAATGGCAGAGACATAATAGTTACCTTTTATGAGTATTTCACGGATGTTTAGAAAAATAGATAAATTTTTCTgtag                                     >  1:2690436/1‑150 (MQ=255)

                                                                                                                 tAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgtagt                                    >  2:2276436/1‑150 (MQ=255)

                                                                                                                       cACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAag                              >  2:1537076/1‑150 (MQ=255)

                                                                                                                        aCTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAaga                             >  1:3884495/1‑150 (MQ=255)

                                                                                                                            tCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAg                         <  2:3470917/150‑1 (MQ=255)

                                                                                                                             cGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGt                        <  1:793301/150‑1 (MQ=255)

                                                                                                                             cGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGt                        <  2:5165476/150‑1 (MQ=255)

                                                                                                                                cTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGTCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTaaa                     <  2:2477293/150‑1 (MQ=255)

                                                                                                                                 tATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAAc                    <  1:1993301/150‑1 (MQ=255)

                                                                                                                                 tATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAAc                    <  1:735312/150‑1 (MQ=255)

                                                                                                                                     gggTGATGAAGGAAGTGAGCATCTGTAATAATATGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACaaaa                >  1:5241522/1‑149 (MQ=255)

                                                                                                                                           tgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatg          >  1:959874/1‑150 (MQ=255)

                                                                                                                                           tgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatg          >  1:5161272/1‑150 (MQ=255)

                                                                                                                                           tgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatg          >  1:3582604/1‑150 (MQ=255)

                                                                                                                                              aggaagTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgcat       <  1:4080362/150‑1 (MQ=255)

                                                                                                                                                   gTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTc  >  1:3625259/1‑150 (MQ=255)

                                                                                                                                                   gTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTc  >  1:1003656/1‑150 (MQ=255)

AAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTTC  >  CP009273/4545937‑4546233

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