BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,310,453	G→C	5.8%	H362Q (CAC→CAG)	PP_1144 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,310,453	0	G	C	5.8%	386.6 / 11.0	121	H362Q (CAC→CAG)	PP_1144	membrane protein
Reads supporting (aligned to +/- strand): ref base G (58/56); new base C (5/2); total (63/58)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 124 total aligned reads displayed.

GTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAC  >  NC_002947/1310307‑1310602

gtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTgg                                                                                                                                                     <  2:1974916/150‑1 (MQ=255)

gtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTgg                                                                                                                                                     <  1:514108/150‑1 (MQ=255)

  ggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGtt                                                                                                                                                   <  1:1561854/150‑1 (MQ=255)

    tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGc                                                                                                                                                 >  1:1945280/1‑150 (MQ=255)

      aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                               >  1:851224/1‑150 (MQ=255)

      aTGACAACGCGCGGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                               <  2:1358341/150‑1 (MQ=255)

         aCCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTa                                                                                                                                            <  1:246249/150‑1 (MQ=255)

            aTGCGCAGCCCGCACAGGTCGGTGGGGTGTTTAAGAAGTTGCATGATGCGGTGGGACAGGAGGCGGGCTTCGTCGCGGCTGTTTAGGCTGGTGACCACACAGAACTCGTCGCCACCACAGCGCGCAACTAGGGCGGGGTTGCGGGTGGCg                                                                                                                                         >  2:1028266/1‑150 (MQ=255)

                gccgccCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCt                                                                                                                                     >  1:896029/1‑150 (MQ=255)

                 ccgcccgcCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCtt                                                                                                                                    >  2:1950528/1‑150 (MQ=255)

                        cccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTgg                                                                                                                             <  2:1268123/150‑1 (MQ=255)

                        cccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGATCAGGCTGGTGACCACACAGGACTCGTCGCCCCCAAAGCGCGCAACCCGGTCGGGGTTGCGGGTGGCGGCCTTGGTGTgg                                                                                                                             >  2:2148571/1‑150 (MQ=255)

                              tCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCg                                                                                                                       <  2:2044312/150‑1 (MQ=255)

                                 aTGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATc                                                                                                                    <  1:1930475/150‑1 (MQ=255)

                                  tGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGGTGGCGaaga                                                                                                                   >  2:961066/1‑147 (MQ=255)

                                   ggggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAc                                                                                                                  >  2:138884/1‑150 (MQ=255)

                                    gggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAcc                                                                                                                 <  1:19634/150‑1 (MQ=255)

                                     ggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAccg                                                                                                                >  1:2430156/1‑149 (MQ=255)

                                      gTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                               <  1:127497/150‑1 (MQ=255)

                                       tCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTg                                                                                                              >  2:321301/1‑150 (MQ=255)

                                        cTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGa                                                                                                             <  1:321301/150‑1 (MQ=255)

                                             ccccGTTCCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                        <  1:1028266/148‑1 (MQ=255)

                                             aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                        <  2:278610/150‑1 (MQ=255)

                                             aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGAGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                        <  2:2421971/150‑1 (MQ=255)

                                              tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgc                                                                                                                                               >  1:351781/1‑109 (MQ=255)

                                              tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgc                                                                                                                                               <  2:351781/110‑2 (MQ=255)

                                                    gCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCAc                                                                                                 >  1:1134072/1‑150 (MQ=255)

                                                     cATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCAcc                                                                                                <  2:2175000/150‑1 (MQ=255)

                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                               <  2:1945280/150‑1 (MQ=255)

                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                               >  2:2217508/1‑150 (MQ=255)

                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                               >  2:2120167/1‑150 (MQ=255)

                                                       tgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCgg                                                                                              >  1:2176236/1‑150 (MQ=255)

                                                         atgCGCTGGGCCAGGTGGCGGGCTTGGTCTCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGcg                                                                                            >  1:1638665/1‑150 (MQ=255)

                                                              cTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑g                                                                                      <  2:851224/150‑1 (MQ=255)

                                                               tGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑gg                                                                                     <  2:2142815/150‑1 (MQ=255)

                                                               tGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑gg                                                                                     >  1:2002133/1‑150 (MQ=255)

                                                                 ggCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGcc                                                                                   >  1:442786/1‑150 (MQ=255)

                                                                   ccAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGc                                                                                                                                       <  2:1265881/97‑1 (MQ=255)

                                                                   ccAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCaa                                                                                 >  2:2329936/1‑150 (MQ=255)

                                                                   ccAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCaa                                                                                 >  1:2003792/1‑150 (MQ=255)

                                                                     aGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAgg                                                                               <  2:1134072/150‑1 (MQ=255)

                                                                             ggCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGtt                                                                       <  2:1554057/150‑1 (MQ=255)

                                                                                    tcgCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGc                                                                <  1:138884/150‑1 (MQ=255)

                                                                                      gcgGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCtt                                                              <  1:2207196/150‑1 (MQ=255)

                                                                                         gCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGaa                                                           <  2:94665/150‑1 (MQ=255)

                                                                                         gCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTAGTTGCGGGTAGCCGCCTTGATGTAGTTGGCGATCTCCTTGAGCTGCTCGGCAGCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGaa                                                           <  2:605378/150‑1 (MQ=255)

                                                                                            gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATg                                                        <  2:1638665/150‑1 (MQ=255)

                                                                                            gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGAGCCTCAGGCTTGCGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATg                                                        <  2:1638598/150‑1 (MQ=255)

                                                                                             ttCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATgg                                                       >  1:176049/1‑150 (MQ=255)

                                                                                              tCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                      <  2:442786/150‑1 (MQ=255)

                                                                                                aGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCg                                                    >  1:2451963/1‑150 (MQ=255)

                                                                                                   cTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATg                                                 >  2:297689/1‑150 (MQ=255)

                                                                                                       tGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGa                                             <  1:297689/150‑1 (MQ=255)

                                                                                                        gACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCgtg                                            <  1:2302265/150‑3 (MQ=255)

                                                                                                          ccACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGa                                          <  2:425232/150‑1 (MQ=255)

                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGa                                                                     <  1:463552/122‑1 (MQ=255)

                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGaa                                         >  2:1962356/1‑150 (MQ=255)

                                                                                                            acacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGCGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGGCGT‑TGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                        <  1:1022887/150‑1 (MQ=255)

                                                                                                              acaGAACTCGTCGCCCCCAAGGCGCGCGACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACAt                                      <  2:389305/150‑1 (MQ=255)

                                                                                                              acaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACAt                                      >  1:38221/1‑150 (MQ=255)

                                                                                                                aGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGa                                    >  2:892682/1‑150 (MQ=255)

                                                                                                                     tcgtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                      <  1:964529/127‑1 (MQ=255)

                                                                                                                      cgtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCa                              <  2:2003792/150‑1 (MQ=255)

                                                                                                                          gCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                        <  1:1016481/73‑1 (MQ=255)

                                                                                                                          gCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                        >  2:1016481/1‑73 (MQ=255)

                                                                                                                          gCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGAAATGCATCCCGGCAGCCTTGAACAGCGCTTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCg                          >  1:1486089/1‑150 (MQ=255)

                                                                                                                           ccaccaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTGCTTCGC‑TAGTAGACTAAGACTTGCGTGGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCgg                         >  1:1861629/1‑150 (MQ=255)

                                                                                                                           ccaccaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCgg                         >  1:1533646/1‑150 (MQ=255)

                                                                                                                            caccaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                        <  2:1684255/134‑1 (MQ=255)

                                                                                                                            caccaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                        >  1:1684255/1‑134 (MQ=255)

                                                                                                                             accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAcc                                                                                                                 >  2:1496959/1‑61 (MQ=255)

                                                                                                                             accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAcc                                                                                                                 <  1:1496959/61‑1 (MQ=255)

                                                                                                                             accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt                       >  2:183496/1‑150 (MQ=255)

                                                                                                                               caAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCg                                              <  1:1376287/125‑1 (MQ=255)

                                                                                                                               caAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCg                                              >  2:1376287/1‑125 (MQ=255)

                                                                                                                               caAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtttt                     <  1:183496/150‑1 (MQ=255)

                                                                                                                                aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc                    <  1:588550/150‑1 (MQ=255)

                                                                                                                                aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc                    >  1:1829872/1‑150 (MQ=255)

                                                                                                                                  aGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc                  <  1:815091/150‑1 (MQ=255)

                                                                                                                                  aGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc                  <  2:1943091/150‑1 (MQ=255)

                                                                                                                                    cccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTg                >  2:2338077/3‑150 (MQ=255)

                                                                                                                                     ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTgg               >  2:204019/2‑150 (MQ=255)

                                                                                                                                     ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTgg               >  2:1229574/2‑150 (MQ=255)

                                                                                                                                      cgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                      >  2:964173/1‑110 (MQ=255)

                                                                                                                                      cgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                      <  1:964173/110‑1 (MQ=255)

                                                                                                                                       gcAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCt                                                               <  1:1273544/100‑1 (MQ=255)

                                                                                                                                        cAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCt            >  2:185767/1‑150 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGAt                                                                    <  1:119887/93‑1 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGAt                                                                    >  2:119887/1‑93 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAg                                            <  1:930957/117‑1 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAg                                            >  2:930957/1‑117 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTg           >  1:1801984/1‑150 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTg           >  2:1815275/1‑150 (MQ=255)

                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGAAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGc                                                                >  1:257292/1‑97 (MQ=255)

                                                                                                                                           ccAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGAt         <  1:2148571/150‑1 (MQ=255)

                                                                                                                                                tCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCg    >  2:51741/1‑150 (MQ=255)

                                                                                                                                                tCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCg    <  2:38221/150‑1 (MQ=255)

                                                                                                                                                  gTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAc  <  2:2451963/150‑1 (MQ=255)

                                                                                                                                                  gTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAc  <  1:2296114/150‑1 (MQ=255)

                                                                                                                                                  gTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAc  <  2:1533646/150‑1 (MQ=255)

GTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGT‑GGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGAC  >  NC_002947/1310307‑1310602

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