BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,802,498	A→T	6.6%	L163Q (CTG→CAG)	vdh ←	vanillin dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,802,498	0	A	T	6.6%	1060.4 / 46.0	335	L163Q (CTG→CAG)	vdh	vanillin dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (170/140); new base T (20/2); total (193/142)
Fisher's exact test for biased strand distribution p-value = 6.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 338 total aligned reads displayed.

GACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGGTCATG  >  NC_002947/3802350‑3802647

gacgaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAg                                                                                                                                                      >  2:6634855/1‑150 (MQ=255)

gacgaTGGCGGGGGCATCCTGCGGCGCATTGCTGAGGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAg                                                                                                                                                      >  2:416731/1‑150 (MQ=255)

     tGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCa                                                                                                                                                 <  1:3847078/150‑1 (MQ=255)

      ggCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCat                                                                                                                                                <  2:327908/150‑1 (MQ=255)

        cGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCCCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatgg                                                                                                                                              <  1:3729961/150‑1 (MQ=255)

          ggggCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggca                                                                                                                                            <  2:5168122/150‑1 (MQ=255)

            ggCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaat                                                                                                                                          <  1:4292/150‑1 (MQ=255)

            ggCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaat                                                                                                                                          <  2:1752994/150‑1 (MQ=255)

              cATGCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaatgg                                                                                                                                        >  2:6518892/1‑150 (MQ=255)

                 ccTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAc                                                                                                                                     <  2:123492/150‑1 (MQ=255)

                  cTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAcg                                                                                                                                    >  1:5876432/1‑150 (MQ=255)

                   tGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAcgc                                                                                                                                   >  2:3348896/1‑150 (MQ=255)

                      ggcgCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTg                                                                                                                                <  2:1119601/150‑1 (MQ=255)

                       gcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTgg                                                                                                                               <  1:3401902/150‑1 (MQ=255)

                          cATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCa                                                                                                                            <  1:3789818/150‑1 (MQ=255)

                            ttGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGt                                                                                                                          <  2:1056843/150‑1 (MQ=255)

                            ttGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGt                                                                                                                          <  2:5282151/150‑1 (MQ=255)

                             tGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTa                                                                                                                         <  2:5043052/150‑1 (MQ=255)

                                tgatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCa                                                                                                                      >  2:2667501/1‑150 (MQ=255)

                                tgatgaCATTGACCACCCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCa                                                                                                                      >  2:5784627/1‑150 (MQ=255)

                                 gatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAc                                                                                                                     >  2:3401695/1‑150 (MQ=255)

                                     aCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                                                                                                                 <  2:4894661/150‑1 (MQ=255)

                                       aTTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCAAGGAAATGGAACGCGTGGCCAGTATCACCGGGGc                                                                                                               >  1:5069514/1‑150 (MQ=255)

                                          gACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGtt                                                                                                            >  2:1943541/1‑150 (MQ=255)

                                           aCCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTc                                                                                                           >  2:4778201/1‑150 (MQ=255)

                                            ccACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTcc                                                                                                          >  1:66174/1‑150 (MQ=255)

                                            acACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGACGATCAGACCATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGAACCACGGTGTTGCAGAAGGCCAGAGGCATGGCAATGGCACGCGTGGAAAGTATCAAAGGGGCGTTcc                                                                                                          >  2:398267/2‑150 (MQ=255)

                                             cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCa                                                                                                         <  2:4664900/150‑1 (MQ=255)

                                             cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCa                                                                                                         >  2:5652845/1‑150 (MQ=255)

                                               cGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACg                                                                                                       <  2:6321000/150‑1 (MQ=255)

                                               cGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACg                                                                                                       >  2:1573061/1‑150 (MQ=255)

                                               cGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGACCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACg                                                                                                       <  2:696376/150‑1 (MQ=255)

                                                  cGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTg                                                                                                    <  2:3760657/150‑1 (MQ=255)

                                                  cGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTg                                                                                                    >  1:2058058/1‑150 (MQ=255)

                                                  cGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTg                                                                                                    >  1:4820799/1‑150 (MQ=255)

                                                   gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGc                                                                                                   <  2:295591/150‑1 (MQ=255)

                                                      gCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGAt                                                                                                <  2:4373698/150‑1 (MQ=255)

                                                        cGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGc                                                                                              >  2:2381188/1‑150 (MQ=255)

                                                         gATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCATTATCACCGGGGCGTTCCACGGTGCGATGcc                                                                                             <  1:1967932/150‑1 (MQ=255)

                                                         gATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGcc                                                                                             <  1:873832/150‑1 (MQ=255)

                                                              cTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCg                                                                                                                                                    <  2:22174/90‑1 (MQ=255)

                                                               tGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCCGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAacac                                                                                       >  1:4671756/1‑150 (MQ=255)

                                                                  aTCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAAcaccac                                                                                    >  2:4307132/1‑150 (MQ=255)

                                                                  aTCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAAcaccac                                                                                    >  1:6333985/1‑150 (MQ=255)

                                                                   tCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCAcg                                                                                   >  1:168249/1‑150 (MQ=255)

                                                                      tGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCAcgccg                                                                                >  2:6203412/1‑150 (MQ=255)

                                                                             ccTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcg                                                                         <  1:4193269/150‑1 (MQ=255)

                                                                                 gCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgcg                                                                     >  2:3126690/1‑150 (MQ=255)

                                                                                    gATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAg                                                                  <  2:2587026/150‑1 (MQ=255)

                                                                                     actaGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAgg                                                                 <  2:3924790/147‑1 (MQ=255)

                                                                                     aTCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAgg                                                                 >  2:3215589/1‑150 (MQ=255)

                                                                                         gCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCa                                                             >  2:4244464/1‑150 (MQ=255)

                                                                                           cGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCGCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                           <  2:3496073/150‑1 (MQ=255)

                                                                                            gATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTg                                                          >  1:629016/1‑150 (MQ=255)

                                                                                              tGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCg                                                        >  2:3014623/1‑150 (MQ=255)

                                                                                               ggACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGa                                                       <  1:4556563/150‑1 (MQ=255)

                                                                                               ggACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGa                                                       >  1:4067859/1‑150 (MQ=255)

                                                                                                gACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGaa                                                      <  2:781950/150‑1 (MQ=255)

                                                                                                gACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGGGGTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGaa                                                      <  2:6354958/150‑1 (MQ=255)

                                                                                                      cgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGcc                                                >  2:1443974/1‑150 (MQ=255)

                                                                                                       gggCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCg                                               <  2:2903819/150‑1 (MQ=255)

                                                                                                       gggCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCg                                               <  1:3523259/150‑1 (MQ=255)

                                                                                                        ggCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCgg                                              <  1:6220516/150‑1 (MQ=255)

                                                                                                        ggCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCgg                                              <  1:5007244/150‑1 (MQ=255)

                                                                                                          cTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGa                                            <  1:2674310/150‑1 (MQ=255)

                                                                                                           tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGac                                           >  1:4000052/1‑149 (MQ=255)

                                                                                                           tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGaa                                           >  1:1435963/1‑150 (MQ=255)

                                                                                                           tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGaa                                           >  1:6431280/1‑150 (MQ=255)

                                                                                                            cAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAAc                                          <  1:5043962/150‑1 (MQ=255)

                                                                                                                tCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCGGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGACCGTcg                                      >  1:1902027/1‑149 (MQ=255)

                                                                                                                  cGAGGCCTTGAGCACCCCGGTGGTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGGTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGa                                    <  2:190915/150‑1 (MQ=255)

                                                                                                                  cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGa                                    <  2:1630388/150‑1 (MQ=255)

                                                                                                                   gAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAg                                   >  1:2377192/1‑150 (MQ=255)

                                                                                                                       ccTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGa                               <  1:3617767/150‑1 (MQ=255)

                                                                                                                        cTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGAt                              >  2:547831/1‑150 (MQ=255)

                                                                                                                               accacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCa                       >  2:4564722/1‑150 (MQ=255)

                                                                                                                                ccacGGTGTTGCCGCAGGCCAGCGGCATGGCAAGGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAc                      <  1:5189788/150‑1 (MQ=255)

                                                                                                                                 cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                     >  1:5602398/1‑150 (MQ=255)

                                                                                                                                 cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                     >  1:2859316/1‑150 (MQ=255)

                                                                                                                                 cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                     >  2:6440095/1‑150 (MQ=255)

                                                                                                                                 cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                     >  1:3252064/1‑150 (MQ=255)

                                                                                                                                 cacGGGGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGGTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                     >  2:3632384/1‑150 (MQ=255)

                                                                                                                                    ggTGTTGGCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGt                  <  1:6045797/150‑1 (MQ=255)

                                                                                                                                      tgtTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGa                >  1:5928977/1‑150 (MQ=255)

                                                                                                                                       gtTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGAt               <  1:1185014/150‑1 (MQ=255)

                                                                                                                                       catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGACCGTCCGAGGGGATCACTTCACCGGTGAt               >  1:3468123/3‑150 (MQ=255)

                                                                                                                                       catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTggt               >  1:2197069/3‑148 (MQ=255)

                                                                                                                                       catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGAt               >  2:444142/3‑150 (MQ=255)

                                                                                                                                        atGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATc              >  1:5101581/2‑150 (MQ=255)

                                                                                                                                        atGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATc              >  1:2526645/2‑150 (MQ=255)

                                                                                                                                         tGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCt             >  2:6124359/1‑150 (MQ=255)

                                                                                                                                         tGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCt             >  2:2767898/1‑150 (MQ=255)

                                                                                                                                           ccGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGc           >  1:6793147/1‑150 (MQ=255)

                                                                                                                                            cGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCg          >  2:3558875/1‑150 (MQ=255)

                                                                                                                                             gCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTTTCACCGGGGCGTTCCACGGTGCGATGCCCAACACCCCGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGCTCTGCGt         <  2:3832476/150‑1 (MQ=255)

                                                                                                                                              cAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGATGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTg        <  2:3531032/150‑1 (MQ=255)

                                                                                                                                                 gCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGGTc     >  1:310748/1‑150 (MQ=255)

                                                                                                                                                   cAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGGTCAt   <  2:4912120/150‑1 (MQ=255)

                                                                                                                                                   cAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGGTCAt   <  2:5971235/150‑1 (MQ=255)

                                                                                                                                                    aGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGGTCATg  >  2:1180511/1‑150 (MQ=255)

GACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGGTCATG  >  NC_002947/3802350‑3802647

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