BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,889,129	C→A	5.6%	A4S (GCA→TCA)	ygcQ ←	putative flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,889,129	0	C	A	5.6%	966.2 / 22.6	251	A4S (GCA→TCA)	ygcQ	putative flavoprotein
Reads supporting (aligned to +/- strand): ref base C (124/113); new base A (10/4); total (134/117)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 255 total aligned reads displayed.

TGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCGGCAGC  >  CP009273/2888981‑2889278

    aCTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGatat                                                                                                                                                  >  2:4363163/1‑150 (MQ=255)

        ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCAt                                                                                                                                              <  2:4088935/150‑1 (MQ=255)

        ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCAt                                                                                                                                              <  2:311353/150‑1 (MQ=255)

        ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCAt                                                                                                                                              <  2:2591999/150‑1 (MQ=255)

        ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCAt                                                                                                                                              >  2:105999/1‑150 (MQ=255)

            accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCt                                                                                                                                          >  1:3108288/1‑150 (MQ=255)

            accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCt                                                                                                                                          >  2:3443158/1‑150 (MQ=255)

             ccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTg                                                                                                                                         >  2:884121/1‑150 (MQ=255)

              caACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGc                                                                                                                                        >  2:4472276/1‑150 (MQ=255)

               aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCa                                                                                                                                       <  1:3929065/150‑1 (MQ=255)

               aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCa                                                                                                                                       >  1:837269/1‑150 (MQ=255)

               aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCa                                                                                                                                       <  1:517457/150‑1 (MQ=255)

                aCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCAt                                                                                                                                      <  1:298646/150‑1 (MQ=255)

                   cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCt                                                                                                                                   <  1:4241150/150‑1 (MQ=255)

                      ccAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTg                                                                                                                                <  1:1621049/150‑1 (MQ=255)

                             cTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                         >  1:1334490/1‑150 (MQ=255)

                                ttCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAg                                                                                                                      <  1:776523/150‑1 (MQ=255)

                                ttCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAg                                                                                                                      >  1:126651/1‑150 (MQ=255)

                                 tCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGt                                                                                                                     <  1:911937/150‑1 (MQ=255)

                                   cGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTcc                                                                                                                   >  1:1173497/1‑150 (MQ=255)

                                    gCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCt                                                                                                                  >  1:2001624/1‑150 (MQ=255)

                                    gCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCt                                                                                                                  >  2:4583891/1‑150 (MQ=255)

                                       aCTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGcc                                                                                                               <  1:3749960/150‑1 (MQ=255)

                                       aCTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCCGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGcc                                                                                                               <  2:2178556/150‑1 (MQ=255)

                                        cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCACTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCa                                                                                                              <  2:525253/150‑1 (MQ=255)

                                         tACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCAt                                                                                                             >  2:2731385/1‑150 (MQ=255)

                                         tACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCAt                                                                                                             >  2:2731395/1‑150 (MQ=255)

                                           cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                           <  1:707527/150‑1 (MQ=255)

                                           cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                           <  2:377163/150‑1 (MQ=255)

                                           cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                           <  1:1428024/150‑1 (MQ=255)

                                             ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCt                                                                                                         <  2:1721228/150‑1 (MQ=255)

                                              gCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                                        >  2:3550815/1‑150 (MQ=255)

                                               cTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCttt                                                                                                       >  2:1245648/1‑150 (MQ=255)

                                                   ggTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAg                                                                                                   <  2:1358489/150‑1 (MQ=255)

                                                   ggTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAg                                                                                                   <  1:2500745/150‑1 (MQ=255)

                                                      tCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCtt                                                                                                >  1:1758533/1‑150 (MQ=255)

                                                           cTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCt                                                                                           <  1:3479154/150‑1 (MQ=255)

                                                             gCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGc                                                                                         <  2:1919110/150‑1 (MQ=255)

                                                              ccAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCt                                                                                        <  2:1652433/150‑1 (MQ=255)

                                                              ccAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCt                                                                                        >  1:4627856/1‑150 (MQ=255)

                                                                aGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTAc                                                                                      >  1:3609141/1‑150 (MQ=255)

                                                                  tGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCg                                                                                    >  2:557794/1‑150 (MQ=255)

                                                                    cgcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGtt                                                                                  >  1:2884041/1‑150 (MQ=255)

                                                                    cgcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGtt                                                                                  <  2:1432232/150‑1 (MQ=255)

                                                                     gcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGttt                                                                                 >  1:4331614/1‑150 (MQ=255)

                                                                           gTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                           <  1:2615700/150‑1 (MQ=255)

                                                                               aaCCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaat                                                                       <  2:1098512/150‑1 (MQ=255)

                                                                               aaCCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaat                                                                       <  1:1430362/150‑1 (MQ=255)

                                                                                 ccGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatca                                                                     <  1:2073190/150‑1 (MQ=255)

                                                                                   gCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGc                                                                   >  2:2013703/1‑150 (MQ=255)

                                                                                   gCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGc                                                                   >  1:183994/1‑150 (MQ=255)

                                                                                   gCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGc                                                                   >  1:2417706/1‑150 (MQ=255)

                                                                                    cTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCg                                                                  <  1:2362755/150‑1 (MQ=255)

                                                                                     tGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGt                                                                 >  2:3728758/1‑150 (MQ=255)

                                                                                     tGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATAGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGt                                                                 >  2:2226977/1‑150 (MQ=255)

                                                                                      gAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTc                                                                >  2:794867/1‑150 (MQ=255)

                                                                                        aaTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGc                                                              >  1:259069/1‑150 (MQ=255)

                                                                                         aTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCt                                                             <  1:1072778/150‑1 (MQ=255)

                                                                                          tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTc                                                            >  2:3111582/1‑150 (MQ=255)

                                                                                           ccTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTcc                                                           <  1:1143660/150‑1 (MQ=255)

                                                                                           ccTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTcc                                                           >  1:3250597/1‑150 (MQ=255)

                                                                                            cTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTCTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCcg                                                          >  1:3508738/1‑150 (MQ=255)

                                                                                              gcgcCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCcgcg                                                        >  1:4413273/1‑150 (MQ=255)

                                                                                                 ccgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcg                                                     <  2:807978/150‑1 (MQ=255)

                                                                                                  cgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcgt                                                    >  2:492735/1‑150 (MQ=255)

                                                                                                   gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGtt                                                   <  2:2089738/150‑1 (MQ=255)

                                                                                                   gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGtt                                                   <  1:3350809/150‑1 (MQ=255)

                                                                                                   gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGtt                                                   <  1:1978997/150‑1 (MQ=255)

                                                                                                   gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCAGCCATAGCTTATGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGACGCTCCGcgtcgat                                                   >  1:1045931/1‑148 (MQ=255)

                                                                                                     cagccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgc                                                 >  1:3190291/1‑150 (MQ=255)

                                                                                                     cagccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTATGCCCATCAATCAGCGTCGCTCCGCGTCGTtgc                                                 >  1:4209124/1‑150 (MQ=255)

                                                                                                          agcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgc                                            >  2:2466321/1‑150 (MQ=255)

                                                                                                           gcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgct                                           <  2:3484133/150‑1 (MQ=255)

                                                                                                             ttGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTgctctgctct                                         >  1:3032417/1‑150 (MQ=255)

                                                                                                              tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCtt                                        >  1:3864961/1‑150 (MQ=255)

                                                                                                                  aTTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCa                                    >  2:2300284/1‑150 (MQ=255)

                                                                                                                      cGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTg                                <  2:3661136/150‑1 (MQ=255)

                                                                                                                         cATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAg                             >  2:1661487/1‑150 (MQ=255)

                                                                                                                         cATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAg                             >  2:1301941/1‑150 (MQ=255)

                                                                                                                          aTTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                            >  1:94038/1‑150 (MQ=255)

                                                                                                                            tctcCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGaca                          >  2:2911877/1‑150 (MQ=255)

                                                                                                                             ctcCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGacac                         <  1:3146322/150‑1 (MQ=255)

                                                                                                                              tcCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                        <  1:2757015/150‑1 (MQ=255)

                                                                                                                              tcCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                        <  2:1184109/150‑1 (MQ=255)

                                                                                                                              tcCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                        >  1:380584/1‑150 (MQ=255)

                                                                                                                                 tGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa                     >  2:3602899/1‑150 (MQ=255)

                                                                                                                                 tGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGATCTTGCCAGTTGCAGACACAGCa                     >  1:3650825/1‑150 (MQ=255)

                                                                                                                                   aTTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc                   <  2:3916401/150‑1 (MQ=255)

                                                                                                                                   aTTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc                   <  1:2773177/150‑1 (MQ=255)

                                                                                                                                   aTTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc                   >  1:739426/1‑150 (MQ=255)

                                                                                                                                    tttATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCg                  >  2:2072610/1‑150 (MQ=255)

                                                                                                                                        ttGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTgg              >  1:2830061/1‑150 (MQ=255)

                                                                                                                                           tCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAt           <  2:2941251/150‑1 (MQ=255)

                                                                                                                                            cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAtt          >  1:4058291/1‑150 (MQ=255)

                                                                                                                                             acaATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAttt         <  1:2178170/150‑1 (MQ=255)

                                                                                                                                              caATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTc        <  2:3366030/150‑1 (MQ=255)

                                                                                                                                              caATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTc        <  1:449336/150‑1 (MQ=255)

                                                                                                                                               aaTTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCg       >  1:436106/1‑150 (MQ=255)

                                                                                                                                                aTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCgg      >  2:3816600/1‑150 (MQ=255)

                                                                                                                                                  tGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCGgca    >  2:4199267/1‑150 (MQ=255)

TGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCGGCAGC  >  CP009273/2888981‑2889278

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