BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,810,512	A→G	6.0%	noncoding (413/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,512	0	A	G	6.0%	448.4 / 12.0	150	noncoding (413/434 nt)	PP_mr35	cspA regulatory region
Reads supporting (aligned to +/- strand): ref base A (64/74); new base G (6/3); total (70/80)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 156 total aligned reads displayed.

CGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810363‑2810658

cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                     >  2:3781564/1‑150 (MQ=255)

cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                     <  2:2888743/150‑1 (MQ=255)

cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                     >  2:560368/1‑150 (MQ=255)

cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCGTCAGATGATGCATCACCCATAAAAAAACCCGCCCCaa                                                                                                                                                     >  2:2795089/1‑150 (MQ=255)

   aaaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTAAGATGATGCATCACCCATAAAAAAACCCGCCCCAAggg                                                                                                                                                  >  2:4154369/1‑150 (MQ=255)

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

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

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

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

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

          gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCCTAAAAAAACCCGCCACAAGGGCGGG‑ttt                                                                                                                                          >  1:2993492/1‑150 (MQ=255)

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

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

           aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑tttt                                                                                                                                         <  2:2893196/150‑1 (MQ=255)

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

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

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

                      cgGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTAt                                                                                                                              <  1:1904315/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                                                                                                                                   >  2:77517/1‑140 (MQ=255)

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

                                 tAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCGATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTc                                                                                                                   <  2:37366/150‑1 (MQ=255)

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

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

                                        gggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGa                                                                                                            <  2:1003096/150‑1 (MQ=255)

                                           ttCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAATAACCCGGCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAAcc                                                                                                         >  1:3115152/1‑150 (MQ=255)

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

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

                                              aCCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTGTTTACAGCTATCCGAATTACTCCGGACGAACCtgt                                                                                                      >  2:398105/1‑150 (MQ=255)

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

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

                                                       cacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATcc                                                                                                                            <  1:3108583/119‑1 (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)

                                                          gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTa                                                                                                                                     <  1:5341663/107‑1 (MQ=255)

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

                                                                aGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTACGGACGAACCTGTGCAGCCTGCATGCCCTTc                                                                                    >  1:498906/1‑150 (MQ=255)

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

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

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

                                                                              aCCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCCACGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGc                                                                      <  2:2027806/150‑1 (MQ=255)

                                                                               cccTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcg                                                                     >  1:2213729/1‑150 (MQ=255)

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

                                                                                  tACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgacg                                                                  >  1:3306837/1‑150 (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)

                                                                                         ggTTTCGCAGGCAATCCCGCCCCATCCGATGATGCATCCCCCATAACCAAACCCGCCACAAGGGCGGG‑TTTTATTACAGCTACCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTACTGGCCGCGCTCGGCGACGAAGGaaa                                                           <  1:2695954/150‑1 (MQ=255)

                                                                                          gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAAc                                                          >  2:1642349/1‑150 (MQ=255)

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

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

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

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

                                                                                                 aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                   >  1:512285/1‑150 (MQ=255)

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

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

                                                                                                          cGCCCCTTCAGCTTTTGCATCACCCATAAAAAAACCCGCCACAAGGGAGTT‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGGCGGGCTCGGCTTCTTAGGAAAAAGTCTGGCCTTCTTTc                                          <  2:4321833/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                ttCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGctc                                    >  1:634916/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)

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

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

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

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

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

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

                                                                                                                                  cATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCg                  >  1:4797814/1‑150 (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)

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

                                                                                                                                     aaaaaaaCCCGCCACACGGGCGGG‑TTTTTTTACAGCTCTCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATg               <  1:94420/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                        aaaaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCt            >  2:1250400/1‑150 (MQ=255)

                                                                                                                                         aaaCCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           >  2:1554127/1‑150 (MQ=255)

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

                                                                                                                                            cccGCCACAAGGGCGGG‑TTTTTTTACAGCGATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGa        <  1:2795089/150‑1 (MQ=255)

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

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

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

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

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

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

CGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGG‑TTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810363‑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: `‑`