BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,810,513	G→T	5.5%	noncoding (412/434 nt)	PP_mr35 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,810,513	0	G	T	5.5%	459.8 / 10.8	147	noncoding (412/434 nt)	PP_mr35	cspA regulatory region
Reads supporting (aligned to +/- strand): ref base G (62/76); new base T (6/3); total (68/79)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

GCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810364‑2810658

gcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAg                                                                                                                                                    <  1:958632/150‑1 (MQ=255)

    aaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCg                                                                                                                                                <  2:2038542/150‑1 (MQ=255)

     aTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                               <  1:1753704/150‑1 (MQ=255)

     aTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                               <  2:4855737/150‑1 (MQ=255)

      tgtgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCggg                                                                                                                                              <  1:3057744/150‑1 (MQ=255)

       gtgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑t                                                                                                                                            <  2:758245/150‑1 (MQ=255)

        tgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑tt                                                                                                                                           <  1:22055/150‑1 (MQ=255)

         gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑ttt                                                                                                                                          <  1:187380/150‑1 (MQ=255)

         gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑ttt                                                                                                                                          >  2:1896618/1‑150 (MQ=255)

           caaGCCTGCGCGGTCGCCTTGTAGGAGCGGGGTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGG‑ttttt                                                                                                                                        <  2:2972198/149‑1 (MQ=255)

           acaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGG‑ttttt                                                                                                                                        <  1:4985363/148‑1 (MQ=39)

           aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑ttttt                                                                                                                                        <  1:2677281/150‑1 (MQ=255)

            aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑tttttt                                                                                                                                       >  2:1225394/1‑150 (MQ=255)

               ccTGCGCGGTCGCCGTGGAGGAGCGGGGTCACCCGCGACCACGGGCAACGCCCGTGCCAGTCACCCTACAGTCGGGTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCCGAAAAAAACCCGCCCGTGTGGCGGG‑TTTTTTTat                                                                                                                                    <  2:685119/150‑2 (MQ=255)

                 tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAg                                                                                                                                  >  1:1984607/1‑150 (MQ=255)

                       gTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATcc                                                                                                                            <  2:3777954/150‑1 (MQ=255)

                       gTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATcc                                                                                                                            <  2:3777907/150‑1 (MQ=255)

                         cGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGa                                                                                                                          <  2:5314503/150‑1 (MQ=255)

                          gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAGAACCCGCCACGAGGGCGGG‑TTGTTTTACAGCTATCCGaa                                                                                                                         >  2:5042498/1‑150 (MQ=255)

                          gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACa                                                                                                                                   <  1:77517/140‑1 (MQ=255)

                          gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACa                                                                                                                                   >  2:77517/1‑140 (MQ=255)

                          gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGaa                                                                                                                         <  2:1622000/150‑1 (MQ=255)

                             ttGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTGCAGAGGATGCATCACCCATAATAAACCCCGCCACAAGGGCGGG‑TTTGTTGACAGCTCTCCgcagtc                                                                                                                      >  2:4049711/1‑145 (MQ=38)

                              tGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATACAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTAc                                                                                                                     <  2:1489831/150‑1 (MQ=255)

                                   gAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCgg                                                                                                                >  1:1003096/1‑150 (MQ=255)

                                    aGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGa                                                                                                               >  2:862550/1‑150 (MQ=255)

                                        ggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCCTAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGACTTACTCCGGACGaa                                                                                                           <  1:3326659/150‑1 (MQ=255)

                                          ttCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCGTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGGCGAGcc                                                                                                         >  1:3536298/1‑150 (MQ=255)

                                           tCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCt                                                                                                        <  1:1550721/150‑1 (MQ=255)

                                               ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCGTCAGCTGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    >  1:967716/1‑150 (MQ=255)

                                                      cacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATcc                                                                                                                            <  1:3108583/119‑1 (MQ=255)

                                                      cacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATcc                                                                                                                            >  2:3108583/1‑119 (MQ=255)

                                                       acGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCa                                                                                            >  2:3610853/1‑150 (MQ=255)

                                                       acGGGAAAAGCCCGTGCCATTCACCCAAAAGTCGGGTTCGCAGCCAAACCCGCCCCATCAGATGAAGCAACACCCAAAAAAAAACCCGCCACAAGCGCTTG‑TATATTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCa                                                                                            >  2:4115183/1‑150 (MQ=255)

                                                        cGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCCGATGATGCATCACCCAGACAAAAACCCGCCACAAGGGCGGG‑TTTTATTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCAt                                                                                           <  2:366604/150‑1 (MQ=255)

                                                        cGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCAt                                                                                           <  2:1821488/150‑1 (MQ=255)

                                                            cAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGccc                                                                                       <  1:4641616/150‑1 (MQ=255)

                                                             aaaGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                                                                                                                   <  2:1562926/90‑1 (MQ=255)

                                                             aaaGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                                                                                                                   >  1:1562926/1‑90 (MQ=255)

                                                                gCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCt                                                                                   >  2:3305036/1‑150 (MQ=255)

                                                                 cccGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TGTTTTTACAGCTATCCGAATTACTCCGGACGAACCGGTGCAGCCTGCCTGCCCTTCTg                                                                                  >  2:4233910/1‑150 (MQ=255)

                                                                  ccGTGCCCTTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTgg                                                                                 <  2:867942/150‑1 (MQ=255)

                                                                   cGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGc                                                                                >  2:1050885/1‑150 (MQ=255)

                                                                           tCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGAAGCAACACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGCATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCACg                                                                        >  1:128343/1‑150 (MQ=255)

                                                                                cTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGTTGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgac                                                                   <  1:3359952/150‑1 (MQ=255)

                                                                                  aCAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgacga                                                                 >  1:4028886/1‑150 (MQ=255)

                                                                                   cAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTCTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGaa                                                                <  2:2781438/150‑1 (MQ=255)

                                                                                    aGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAg                                                               >  2:971715/1‑150 (MQ=255)

                                                                                      tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGa                                                             >  2:3458208/1‑150 (MQ=255)

                                                                                      tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTGCTGGCCGCGCTCGGCGACGAAGGa                                                             >  1:1936872/1‑150 (MQ=255)

                                                                                        ggTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                           <  1:2695960/150‑1 (MQ=255)

                                                                                            tCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTGCTGGCCGCGCTCGGCGACGAAGGAAACAGt                                                       >  1:3975775/1‑150 (MQ=255)

                                                                                              gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  1:3245615/1‑150 (MQ=255)

                                                                                               cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  1:5079540/150‑1 (MQ=255)

                                                                                               cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  1:5079333/150‑1 (MQ=255)

                                                                                               cAGGCAAGCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCGGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  2:465786/150‑1 (MQ=255)

                                                                                                 ggCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGc                                                  <  2:2031640/150‑1 (MQ=255)

                                                                                                  gCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGcc                                                 >  1:4284578/1‑150 (MQ=255)

                                                                                                  gCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTATTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGcc                                                 <  1:5145915/150‑1 (MQ=255)

                                                                                                        ccGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgac                                                                   <  1:346559/126‑1 (MQ=255)

                                                                                                         cGCCCCTTCAGATGCTGCATCACCCATCAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAAGTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGACACAGTCTGGCCTTCTTTc                                          <  2:4321680/150‑1 (MQ=255)

                                                                                                         cGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTc                                          <  2:608565/150‑1 (MQ=255)

                                                                                                         cGCCCCTTCAGATGATGCATCACCCAAAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTc                                          <  2:2141884/150‑1 (MQ=255)

                                                                                                         cGCCCCTTCAGATGATGCATCACCCAAAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTc                                          <  2:1912025/150‑1 (MQ=255)

                                                                                                           ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        >  1:3617990/1‑150 (MQ=255)

                                                                                                              cTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCgg                                                                                                                <  2:2180023/75‑1 (MQ=255)

                                                                                                              cTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCgg                                                                                                                >  1:2180383/1‑75 (MQ=255)

                                                                                                                   gatgatGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGCCGGACACCGTCTGGCCGTCTTTCAGGCTCTTGa                                >  2:4352215/1‑150 (MQ=255)

                                                                                                                    atgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               <  2:2993492/150‑1 (MQ=255)

                                                                                                                        tgCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCg                           >  2:3488656/1‑150 (MQ=255)

                                                                                                                         gCATCACCCATACAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                          <  1:560368/150‑1 (MQ=255)

                                                                                                                         gCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGCAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                          >  2:4458663/1‑150 (MQ=255)

                                                                                                                         gCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                          <  2:2796555/150‑1 (MQ=255)

                                                                                                                         gCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                          <  1:560456/150‑1 (MQ=255)

                                                                                                                           aTCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTGTCAGGCTCTTGAAGCCGTCa                        >  1:1160726/1‑150 (MQ=255)

                                                                                                                            atagctgtAAAAAAACCCGCCCTTGTGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCGGGCCGCGCTCGGCGACGAAGGCAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAg                       >  1:917289/8‑150 (MQ=38)

                                                                                                                                ccATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTc                   >  2:2496048/1‑150 (MQ=255)

                                                                                                                                 tgtAAAAAAACCCGCCGGTGTGGCGGG‑TTTTTTTACAGCTATCCGAAGTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCg                  >  2:3032337/3‑150 (MQ=39)

                                                                                                                                 tgtAAAAAAACCCGCCCTTGTGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCg                  >  1:1264226/3‑150 (MQ=39)

                                                                                                                                 tgtAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               <  2:2630913/135‑1 (MQ=255)

                                                                                                                                 tgtAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               >  1:2630913/3‑137 (MQ=255)

                                                                                                                                  aTAAAAAATCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGa                 >  1:105161/1‑150 (MQ=255)

                                                                                                                                  aTAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCTCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGa                 <  1:2942703/150‑1 (MQ=255)

                                                                                                                                  aTAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGa                 <  2:1984607/150‑1 (MQ=255)

                                                                                                                                   tAACAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGAt                <  1:2777186/150‑1 (MQ=255)

                                                                                                                                   tAAAAAAACCCGCCCTTGTGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACACTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGAt                >  1:4034499/1‑150 (MQ=39)

                                                                                                                                    aaaaaaaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATg               <  1:3147272/150‑1 (MQ=255)

                                                                                                                                     aaaaaaCCCGCCCTTGTGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATgg              >  1:4065759/1‑150 (MQ=39)

                                                                                                                                      aaaaaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  1:3694920/1‑150 (MQ=255)

                                                                                                                                      aaaaaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  2:3401734/1‑150 (MQ=255)

                                                                                                                                      aaaaaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             <  1:4939325/150‑1 (MQ=255)

                                                                                                                                         aaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCttt          <  2:4797814/150‑1 (MQ=255)

                                                                                                                                            ccGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       <  2:3694920/150‑1 (MQ=255)

                                                                                                                                            ccGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       >  2:502362/1‑150 (MQ=255)

                                                                                                                                            ccGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAAGTCCTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa       <  2:4554979/150‑1 (MQ=255)

                                                                                                                                             cGCCACAAGGGCGGGTTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCcgc                                                                             >  1:5140374/1‑80 (MQ=255)

                                                                                                                                             cGCCACAAGGGCGGGTTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCcgc                                                                             <  2:5140374/80‑1 (MQ=255)

                                                                                                                                               ccACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAgtg    >  1:1125527/1‑150 (MQ=255)

                                                                                                                                                cacaAGGGCGGG‑TATTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGTAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAgtgt   <  1:2883566/150‑1 (MQ=255)

                                                                                                                                                 acacGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  2:3258350/150‑1 (MQ=255)

GCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810364‑2810658

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