BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,121,729	G→C	5.1%	A307A (GCC→GCG)	PP_2736 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,121,729	0	G	C	5.1%	728.3 / 20.7	215	A307A (GCC→GCG)	PP_2736	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (98/106); new base C (2/9); total (100/115)
Fisher's exact test for biased strand distribution p-value = 6.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 220 total aligned reads displayed.

GGGCGATCACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTGGCAAGCGAA  >  NC_002947/3121581‑3121878

gggCGATCACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGACTGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTgg                                                                                                                                                      >  2:5094040/1‑150 (MQ=255)

gggCGATCACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTgg                                                                                                                                                      <  2:4827066/150‑1 (MQ=255)

    gATCACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCa                                                                                                                                                  <  2:2456769/150‑1 (MQ=255)

      tCACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAgg                                                                                                                                                <  1:7488880/150‑1 (MQ=255)

       cACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACTTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGc                                                                                                                                               <  2:2500119/150‑1 (MQ=255)

       cACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACTTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGc                                                                                                                                               <  2:7117964/150‑1 (MQ=255)

       cACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGc                                                                                                                                               <  1:6339229/150‑1 (MQ=255)

       cACCAGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGc                                                                                                                                               <  2:3799407/150‑1 (MQ=255)

          cTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAg                                                                                                                                            >  2:8565143/1‑150 (MQ=255)

           tGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGc                                                                                                                                           >  1:3594603/1‑150 (MQ=255)

              gcCGGGTCTATCAGTGAGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTgg                                                                                                                                        >  2:3890322/1‑150 (MQ=255)

               ccGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTCGCCTGGGCCAGCTgga                                                                                                                                       <  2:4378560/150‑2 (MQ=255)

               ccGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCAAGGCGGTAGTTCCAGCTCGCCTGGGCCAGCTgga                                                                                                                                       <  2:7184506/150‑2 (MQ=255)

                cGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTCGCCTGGGCCAGCTggaa                                                                                                                                      <  2:6856339/150‑3 (MQ=255)

                 gggTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTCGCCTGGGCCAGCTggaac                                                                                                                                     <  2:8113978/150‑4 (MQ=255)

                  ggTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTggcggc                                                                                                                                    <  1:2615400/150‑1 (MQ=255)

                  ggTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTggcggc                                                                                                                                    <  1:4037796/150‑1 (MQ=255)

                    tCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGt                                                                                                                                  >  1:1369374/1‑150 (MQ=255)

                    tCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGt                                                                                                                                  <  1:8082382/150‑1 (MQ=255)

                    tCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGt                                                                                                                                  <  2:7717868/150‑1 (MQ=255)

                       aTCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGgc                                                                                                                               <  2:5612253/150‑1 (MQ=255)

                       aTCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGgc                                                                                                                               >  2:3485686/1‑150 (MQ=255)

                       aTCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGgc                                                                                                                               <  1:2799357/150‑1 (MQ=255)

                          aGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGgcagc                                                                                                                            >  2:5982516/1‑150 (MQ=255)

                                     gTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATc                                                                                                                 <  1:4458924/150‑1 (MQ=255)

                                     gTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATc                                                                                                                 >  1:744570/1‑150 (MQ=255)

                                      ttGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCg                                                                                                                <  2:3292591/150‑1 (MQ=255)

                                      ttGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCg                                                                                                                >  1:6668063/1‑150 (MQ=255)

                                         aGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTa                                                                                                             >  2:1487337/1‑150 (MQ=255)

                                           gCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATg                                                                                                           >  1:7517966/1‑150 (MQ=255)

                                           gCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATg                                                                                                           >  1:4129626/1‑150 (MQ=255)

                                             tCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATgca                                                                                                         >  2:6853048/1‑150 (MQ=255)

                                                    gAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGcaccac                                                                                                  <  1:1743936/150‑1 (MQ=255)

                                                    gAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGcaccac                                                                                                  <  2:3567222/150‑1 (MQ=255)

                                                    gAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGcaccac                                                                                                  <  2:1623093/150‑1 (MQ=255)

                                                       ggTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCAcgtc                                                                                               >  2:2006797/1‑150 (MQ=255)

                                                       ggTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGGTGTGGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCAcgtc                                                                                               >  2:5363927/1‑150 (MQ=255)

                                                         ttGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCAcgtcgt                                                                                             <  1:6569928/150‑1 (MQ=255)

                                                           gTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTg                                                                                           <  1:3008592/150‑1 (MQ=255)

                                                            tCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGc                                                                                          <  1:3689574/150‑1 (MQ=255)

                                                              ggTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTg                                                                                        >  1:5924599/1‑150 (MQ=255)

                                                              ggTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTg                                                                                        >  1:8109007/1‑150 (MQ=255)

                                                              ggTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTg                                                                                        <  1:2213690/150‑1 (MQ=255)

                                                               gTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTgg                                                                                       >  2:7833711/1‑150 (MQ=255)

                                                                  ccTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGt                                                                                    <  2:2333081/150‑1 (MQ=255)

                                                                      gATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGt                                                                                                              <  2:2974440/120‑1 (MQ=255)

                                                                      gATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATc                                                                                <  2:2143922/150‑1 (MQ=255)

                                                                      gATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATc                                                                                <  2:6491222/150‑1 (MQ=255)

                                                                        tACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGa                                                                              >  1:786895/1‑150 (MQ=255)

                                                                          ccTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATg                                                                            <  1:6292522/150‑1 (MQ=255)

                                                                            ttGCAGAATGTTGATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCGCTGGCGAGCTGGCGGCGTGGCAGCAGTCGGGTATCGGTATGCAGCACGACGTCGTTGCTGGCGTAATCGATGgc                                                                          >  2:6174996/1‑150 (MQ=255)

                                                                            ttGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGgc                                                                          >  2:522126/1‑150 (MQ=255)

                                                                               cAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCAGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCt                                                                       >  2:4252416/1‑150 (MQ=255)

                                                                                 gAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGa                                                                     <  2:3908446/150‑1 (MQ=255)

                                                                                 gAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGa                                                                     <  2:7888667/150‑1 (MQ=255)

                                                                                  aaTGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAg                                                                    >  1:6409421/1‑150 (MQ=255)

                                                                                    tGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGTCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTa                                                                  <  2:4507916/150‑1 (MQ=255)

                                                                                      ttCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACg                                                                >  2:5159628/1‑150 (MQ=255)

                                                                                        cATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGgc                                                              >  1:7140658/1‑150 (MQ=255)

                                                                                        cATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCA‑CACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGgcg                                                             >  1:2084782/1‑150 (MQ=255)

                                                                                           gttgtAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGgcgcg                                                           >  2:8927806/1‑150 (MQ=255)

                                                                                             tgtAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCt                                                         >  1:3782887/1‑150 (MQ=255)

                                                                                                aGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCt                                                      <  1:3485686/150‑1 (MQ=255)

                                                                                                aGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCt                                                      <  1:4922396/150‑1 (MQ=255)

                                                                                                aGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCt                                                      <  1:3376481/150‑1 (MQ=255)

                                                                                                  gTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTgg                                                    >  2:6095064/1‑150 (MQ=255)

                                                                                                   tCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGt                                                   >  2:2551939/1‑150 (MQ=255)

                                                                                                   tCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGAAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCTATGGCGCTGAGTACGGCGCGCTCCTGGt                                                   >  1:2130365/1‑150 (MQ=255)

                                                                                                        gcggcgGGAGCCCGCTCGGGGCCGCCCAGGCGTTAGTTCCATCTGGCCCAGGCTAGCTGGCGGCTTTGCAGCAGACGTGTATCGGTATGCAGCACCTCGTCGTTTCTGGCTTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCtt                                              <  2:5150959/150‑2 (MQ=255)

                                                                                                          ggcggGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTggg                                            <  2:7726811/150‑1 (MQ=255)

                                                                                                           gcggGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTgggg                                           <  2:7465302/150‑1 (MQ=255)

                                                                                                             gggAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGtt                                         >  1:1628244/1‑150 (MQ=255)

                                                                                                                ctaccGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGtt                                                  <  1:2326796/135‑1 (MQ=255)

                                                                                                                aGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTggagaaccc                                      <  1:7756701/150‑8 (MQ=255)

                                                                                                                  cccGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCca                                    >  1:2365940/1‑150 (MQ=255)

                                                                                                                  cccGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCca                                    <  1:6149472/150‑1 (MQ=255)

                                                                                                                  cccGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCca                                    <  2:1232690/150‑1 (MQ=255)

                                                                                                                   ccGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTCTCGGTATGCAGCACCACGTCGTGGCTGGCGGCCTCGCTGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTGCTCcag                                   >  2:8304717/1‑150 (MQ=255)

                                                                                                                    cGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCTGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCcagc                                  <  1:5535879/150‑1 (MQ=255)

                                                                                                                    cGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCcagc                                  >  1:1441374/1‑150 (MQ=255)

                                                                                                                     gCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCcagca                                 <  1:4806621/150‑1 (MQ=255)

                                                                                                                       tCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGt                               <  1:5123637/150‑1 (MQ=255)

                                                                                                                          gggCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGCTGCTGGAGAACCCcagcaaccag                            <  2:6334567/150‑6 (MQ=255)

                                                                                                                               gccCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGc                       >  1:8905846/1‑150 (MQ=255)

                                                                                                                                    ggCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCACTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGAt                  >  1:675858/1‑150 (MQ=255)

                                                                                                                                     gCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATc                 <  2:128598/150‑1 (MQ=255)

                                                                                                                                     gCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATc                 <  1:5583017/150‑1 (MQ=255)

                                                                                                                                     gCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATc                 >  1:914825/1‑150 (MQ=255)

                                                                                                                                       ggTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGc               >  1:7623272/1‑150 (MQ=255)

                                                                                                                                          aGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCtgt            <  1:6201718/150‑1 (MQ=255)

                                                                                                                                           gTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCtgtg           <  1:3890322/150‑1 (MQ=255)

                                                                                                                                           gTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCtgtg           <  2:3594762/150‑1 (MQ=255)

                                                                                                                                           gTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCtgtg           <  2:3594199/150‑1 (MQ=255)

                                                                                                                                            ttCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTgg          >  1:2757575/1‑150 (MQ=255)

                                                                                                                                            ttCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTgg          >  1:2672968/1‑150 (MQ=255)

                                                                                                                                            ttCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTgg          >  1:8593430/1‑150 (MQ=255)

                                                                                                                                            cgccAGCTCGCCTGGGCCAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTCTGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCCGGGCCTGATCGCTGTgg          >  1:103795/3‑150 (MQ=255)

                                                                                                                                                aGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTGGCAAg      >  2:5740145/1‑150 (MQ=255)

                                                                                                                                                 gCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTGGCAAGc     <  2:4129626/150‑1 (MQ=255)

                                                                                                                                                    ggCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTGGCAAGCGaa  <  1:7074068/150‑1 (MQ=255)

GGGCGATCACCTGCGCCGGGTCTATCAGTGCGGTCTGGTTGAGGCTCACGCAGAAGGTTGTCGGTGCCTGGATACCTTGCAGAATGTTCATGTTGTAGGTCAGCGCGGCGGGAGCCCGCTCGGGGCCGCCCAGGCGGTAGTTCCAGCTGGCCCAGGCGAGCTGGCGGCGTGGCAGCAGACGGGTATCGGTATGCAGCACCACGTCGTTGCTGGCGTAATCGATGGCGCTGAGTACGGCGCGCTCCTGGTTGCTGGGGTTCTCCAGCAGTGCCAGGGCCTGATCGCTGTGGCAAGCGAA  >  NC_002947/3121581‑3121878

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