BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,875,621	G→A	5.6%	G146R (GGG→AGG)	cigR →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,875,621	0	G	A	5.6%	518.7 / 12.5	161	G146R (GGG→AGG)	cigR	membrane protein
Reads supporting (aligned to +/- strand): ref base G (67/85); new base A (5/4); total (72/89)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCC  >  NC_002947/2875472‑2875770

gCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCg                                                                                                                                                       >  1:1817400/1‑150 (MQ=255)

    tCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCACGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCACGCGGGCGCAGACCTGATGCTGGTCGCCGTGGCGACCGGGat                                                                                                                                                   >  2:2684395/1‑150 (MQ=255)

      agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGatca                                                                                                                                                 >  1:803451/1‑150 (MQ=255)

      agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGatca                                                                                                                                                 >  1:285722/1‑150 (MQ=255)

      agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGAGTCTGGTCGCCGTGGCGACCGGGatca                                                                                                                                                 >  2:4295972/1‑150 (MQ=255)

       gaagaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGACCAGAatcag                                                                                                                                                <  2:4001653/150‑2 (MQ=255)

         agaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCt                                                                                                                                              <  2:849072/150‑1 (MQ=255)

         agaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCt                                                                                                                                              <  2:763669/150‑1 (MQ=255)

            aCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACAGGGATCATCTAtt                                                                                                                                           <  1:978171/150‑2 (MQ=255)

             ccTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGa                                                                                                                                          <  2:940353/150‑1 (MQ=255)

              cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAg                                                                                                                                         <  2:3327688/150‑1 (MQ=255)

                ggCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgt                                                                                                                                       >  1:1847681/1‑150 (MQ=255)

                 gCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgtg                                                                                                                                      <  1:1276019/150‑1 (MQ=255)

                        gCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACg                                                                                                                               >  1:1829500/1‑150 (MQ=255)

                        gCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACg                                                                                                                               >  1:2927844/1‑150 (MQ=255)

                         cAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACgg                                                                                                                              >  2:16627/1‑150 (MQ=255)

                          aaGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGTCGACCGGGATCATCTATGAGGTGTTGAACggg                                                                                                                             <  2:1579082/150‑1 (MQ=255)

                          aaGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACggg                                                                                                                             <  1:3106009/150‑1 (MQ=255)

                          aaGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACggg                                                                                                                             <  2:1866922/150‑1 (MQ=255)

                           aGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACgggg                                                                                                                            <  1:3704411/150‑1 (MQ=255)

                               ccTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCAtt                                                                                                                        >  1:2688461/1‑150 (MQ=255)

                                cTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTc                                                                                                                       >  1:572626/1‑150 (MQ=255)

                                  gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                                     <  1:1054819/150‑1 (MQ=255)

                                  gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                                     <  1:1417987/150‑1 (MQ=255)

                                  gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                                     <  2:996397/150‑1 (MQ=255)

                                   ccACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGAc                                                                                                                    <  2:2974676/150‑1 (MQ=255)

                                     aCCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACAACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTa                                                                                                                  >  1:654223/1‑150 (MQ=255)

                                      ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTaa                                                                                                                 <  2:100592/150‑1 (MQ=255)

                                      ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTaa                                                                                                                 <  1:3954376/150‑1 (MQ=255)

                                      ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTaa                                                                                                                 <  2:2296134/150‑1 (MQ=255)

                                       cTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAg                                                                                                                <  1:491577/150‑1 (MQ=255)

                                         ggTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGc                                                                                                              <  2:572626/150‑1 (MQ=255)

                                         ggTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGc                                                                                                              <  1:404360/150‑1 (MQ=255)

                                          gTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGcc                                                                                                             >  2:2307675/1‑150 (MQ=255)

                                              cGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAg                                                                                                         <  1:2197821/150‑1 (MQ=255)

                                               gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAgg                                                                                                        >  2:2551936/1‑150 (MQ=255)

                                               gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAgg                                                                                                        >  1:4025931/1‑150 (MQ=255)

                                               gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTGGAACGGGGCATTCGACTAAGGCCGCAgg                                                                                                        >  2:2550156/1‑150 (MQ=255)

                                                  aagaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGc                                                                                                     >  2:1473992/1‑150 (MQ=255)

                                                    gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGTGCGCAGACCTGATTCTGGTCGCCTTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                                   <  2:32655/150‑1 (MQ=255)

                                                    gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                                   <  1:4235419/150‑1 (MQ=255)

                                                    gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                                   <  1:765441/150‑1 (MQ=255)

                                                    gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                                   <  2:899435/150‑1 (MQ=255)

                                                    gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                                   <  1:1281538/150‑1 (MQ=255)

                                                    gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGAGGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                                   <  1:2032340/150‑1 (MQ=255)

                                                     aagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                                  <  1:1898340/150‑1 (MQ=255)

                                                     aagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                                  <  2:1777145/150‑1 (MQ=255)

                                                              ggTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCt                                                                                         <  1:2749540/150‑1 (MQ=255)

                                                                   tcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCGGACTGCTCgcac                                                                                    <  2:2853996/150‑4 (MQ=255)

                                                                   tcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTc                                                                                    >  2:3357385/1‑150 (MQ=255)

                                                                    cTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                                  >  1:3236741/1‑135 (MQ=255)

                                                                     tGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTccc                                                                                  <  2:3628575/150‑1 (MQ=255)

                                                                      gTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCt                                                                                 >  2:1607559/1‑150 (MQ=255)

                                                                             cAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGgcg                                                                          <  2:664397/150‑1 (MQ=255)

                                                                              aGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGgcgc                                                                         >  2:1761406/1‑150 (MQ=255)

                                                                              aGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGgcgc                                                                         >  2:4035960/1‑150 (MQ=255)

                                                                               gCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCt                                                                        >  1:4338733/1‑150 (MQ=255)

                                                                                    cTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGt                                                                   >  1:3091282/1‑150 (MQ=255)

                                                                                     tCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTg                                                                  >  2:2286501/1‑150 (MQ=255)

                                                                                        cTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCt                                                               <  2:2408784/150‑1 (MQ=255)

                                                                                          acgacgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGc                                                             <  2:505638/150‑1 (MQ=255)

                                                                                          acgacgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGc                                                             <  1:3655834/150‑1 (MQ=255)

                                                                                            gacgGGTATGAGTGGAGGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGa                                                           <  2:3076440/150‑1 (MQ=255)

                                                                                              cgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATg                                                         <  2:2298254/150‑1 (MQ=255)

                                                                                                 gTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGc                                                      >  2:2008106/1‑150 (MQ=255)

                                                                                                     gAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTgcg                                                  >  1:3119109/1‑150 (MQ=255)

                                                                                                       gTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgca                                                >  1:2593608/1‑150 (MQ=255)

                                                                                                         ggATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTgctgcgc                                              <  2:2305064/150‑6 (MQ=255)

                                                                                                          gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                             >  2:543980/1‑150 (MQ=255)

                                                                                                            tGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCagcagc                                           <  2:3904250/150‑1 (MQ=255)

                                                                                                               aaGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCa                                        <  1:1735345/150‑1 (MQ=255)

                                                                                                               aaGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCa                                        <  1:1574364/150‑1 (MQ=255)

                                                                                                               aaGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCa                                        >  2:3871321/1‑150 (MQ=255)

                                                                                                                 gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCatc                                      <  1:2307675/150‑3 (MQ=255)

                                                                                                                 gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAAc                                      >  1:639535/1‑150 (MQ=255)

                                                                                                                 gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCGATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCGCGGTCCCTGCAGGCGCTGGCGTGCCTGCGCTGGGCTGCGCAGCAGcccgcc                                      >  1:1068062/1‑147 (MQ=255)

                                                                                                                   gggCGCAGACCTTATTCTGGTCGCGGTGGCGACGGGTATCTTTTATGAGGTGTTGAACGGGGCTTTCGACTAAGGCCGCAGTTGCGAGCAGTCCGCTGGGTCCCTGCTGGCGCTGGGGTGCCTGCTTTGGGCTGCGCAGCAGCCCAACgg                                    <  2:1236319/150‑1 (MQ=255)

                                                                                                                   gggCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTTTTGAACGGGGCATTCGACTAAGGCCGCTGGTGCGAGCTGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACgg                                    <  2:4018752/150‑1 (MQ=255)

                                                                                                                    ggCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCATCgca                                   <  1:2044725/150‑3 (MQ=255)

                                                                                                                       tcatagatGATCCCGGTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTa                                >  1:2602281/8‑150 (MQ=255)

                                                                                                                       gcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGc                                                      <  1:1677739/128‑1 (MQ=255)

                                                                                                                       gcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGc                                                      >  2:1677739/1‑128 (MQ=255)

                                                                                                                         aGACCTTTTTCTGGTCGCGGTGGCGACGGGTATAATCTTTGTGGTTTTTAAGGGGGATTTCTACTAAGGCCGCTGGTGCGAGCAGTCCGCTCGGTCCGTGCTGGCGGTGGGGTGCCTGCTTTGGGCTGCGCAGCAGCCCAACGGCCTATc                              <  2:3584179/150‑1 (MQ=255)

                                                                                                                         aGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATc                              >  1:2317929/1‑150 (MQ=255)

                                                                                                                            ccTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAAc                           >  2:2599484/1‑150 (MQ=255)

                                                                                                                             cTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACg                          >  1:2808182/1‑150 (MQ=255)

                                                                                                                             cTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACg                          >  1:369065/1‑150 (MQ=255)

                                                                                                                                 tcccggTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc                      >  2:921172/5‑150 (MQ=255)

                                                                                                                                  tCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGcc                     >  2:1771491/1‑150 (MQ=255)

                                                                                                                                   ccggTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCa                                                                                                          <  1:2723246/62‑1 (MQ=37)

                                                                                                                                    tGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGt                   >  2:1762463/1‑150 (MQ=255)

                                                                                                                                    tGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGt                   >  2:1762523/1‑150 (MQ=255)

                                                                                                                                    cggTCGCCACGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGt                   >  1:1293362/2‑150 (MQ=255)

                                                                                                                                            gTGGCGACCGGGTTCATCTTTGAGGTGTTGACCGGGGCATTGGACTACGGCCGCCGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGGGTGCCTGGGATGGGGTGCGCAGCTGCCCAACGGCCTATCAACGCTGCCGTAACACcgc           <  2:2955418/150‑1 (MQ=255)

                                                                                                                                             tGGCGACCGTTATAATCTATGTGTTTTTGAAGGGGGAATTCTACTAAGGCCGCTGGTGCGAGAAGTCCGCTGGGTCCCTGATGGGGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACcgcg          <  2:2688461/150‑1 (MQ=255)

                                                                                                                                             gggCGACGGGCATAATCTTTGTGGTGTTGCCCGGGGCATTCGACTCCGGCCGCCGGGGCGCGCCGTCGGCGCGGTCCCTGCAGGCGGGGGGGTGCCGGCTTTGGGCTGCCCCGCCCCCACACGGCCTATCAACGCTGCCGTAACACcgcg          <  2:385802/149‑1 (MQ=255)

                                                                                                                                                cGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGcc       >  1:3458227/1‑150 (MQ=255)

                                                                                                                                                   ccGGGATCATCTTTGTGGTGTTGAACGGGGCATTCGACTAAGGCCGCCGGTGCGAGCAGGCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGc    <  2:2927844/150‑1 (MQ=255)

                                                                                                                                                   ccGGGATCATCTATGAGGTGTTGAAGGGGGCATTCGACTACGGCCGCAGGTGCTAGCTGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGc    <  2:1829409/150‑1 (MQ=255)

                                                                                                                                                    cGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTACGGCCGCCGGTGCGCGCAGTCCGCTCGGCCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGcc   <  2:285722/150‑1 (MQ=255)

GCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCC  >  NC_002947/2875472‑2875770

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