BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,838,027	A→T	6.5%	R527W (AGG→TGG)	PP_3388 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,838,027	0	A	T	6.5%	322.8 / 11.6	108	R527W (AGG→TGG)	PP_3388	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (59/42); new base T (4/3); total (63/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

Only 100 of 108 total aligned reads displayed.

CAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACA  >  NC_002947/3837879‑3838174

cAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAg                                                                                                                                                    <  2:5488/150‑1 (MQ=40)

          gCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGcc                                                                                                                                             <  1:1848041/147‑1 (MQ=255)

          gCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGcc                                                                                                                                             >  2:1848041/1‑147 (MQ=255)

          gCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCAt                                                                                                                                          >  1:2313529/1‑150 (MQ=255)

          gCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGATATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCAt                                                                                                                                          >  1:1842830/1‑150 (MQ=255)

            aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCTGGCAAGcgcttcc                                                                                                                                        <  1:1808278/150‑7 (MQ=255)

            aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCTGGCAAGcgcttcc                                                                                                                                        >  2:1808278/1‑144 (MQ=255)

            aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGc                                                                                                                                        <  1:1032423/150‑1 (MQ=255)

            aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGc                                                                                                                                        >  1:30156/1‑150 (MQ=255)

              aCCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCg                                                                                                                                      >  1:517921/1‑150 (MQ=255)

                  gggCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGc                                                                                                                                  <  2:517921/150‑1 (MQ=255)

                         cTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGcc                                                                                                                           >  1:1564065/1‑150 (MQ=255)

                         cTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGACATAATTAAGACAACGTCAGCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGcc                                                                                                                           >  1:1223853/1‑150 (MQ=14)

                          tGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCTGGCAAGcgc                                                                                                                                            >  1:727944/1‑130 (MQ=255)

                          tGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCTGGCAAGcgc                                                                                                                                            <  2:727944/132‑3 (MQ=255)

                          tGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCaa                                                                                                                               <  2:2073556/145‑1 (MQ=255)

                          tGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCt                                                                                                                          >  1:1782597/1‑150 (MQ=255)

                            gCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCttt                                                                                                                        >  1:588248/1‑150 (MQ=255)

                             cAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTg                                                                                                                       <  1:1156181/150‑1 (MQ=255)

                              aaTACGCGGCCTTACGTCGACCCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTgg                                                                                                                      <  2:1321537/150‑1 (MQ=14)

                              aaTACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCAt                                                                                                                                          <  1:177148/130‑1 (MQ=255)

                              aaTACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCAt                                                                                                                                          >  2:177148/1‑130 (MQ=255)

                              aaTACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTgg                                                                                                                      <  1:1387256/150‑1 (MQ=255)

                              aaTACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTgg                                                                                                                      <  1:2427711/150‑1 (MQ=255)

                               aTACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTggg                                                                                                                     <  2:1564065/150‑1 (MQ=255)

                               aTACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTggg                                                                                                                     <  2:1350657/150‑1 (MQ=255)

                                tACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGgc                                                                                                                    <  2:237458/150‑1 (MQ=255)

                                   gcgGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCa                                                                                                                 >  1:2385617/1‑150 (MQ=255)

                                     ggCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATc                                                                                                               <  1:2307852/150‑1 (MQ=255)

                                        aTCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAac                                                                                                            >  2:611687/1‑150 (MQ=255)

                                          cACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAacac                                                                                                          <  1:1298060/150‑1 (MQ=255)

                                            ccTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACAccc                                                                                                        >  1:599350/1‑150 (MQ=255)

                                             cTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCt                                                                                                       >  2:736162/1‑150 (MQ=255)

                                             cTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCt                                                                                                       >  2:2398037/1‑150 (MQ=255)

                                             cTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCt                                                                                                       >  1:1566562/1‑150 (MQ=255)

                                             cTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATAAAAAAAGGCGCATCAACACCCt                                                                                                       >  2:2275531/1‑150 (MQ=14)

                                              tACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTg                                                                                                      >  2:526461/1‑150 (MQ=255)

                                                 aaCCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGGCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTAc                                                                                                   <  2:30156/150‑1 (MQ=255)

                                                 aaCCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTAc                                                                                                   >  1:2286036/1‑150 (MQ=255)

                                                    cAGGAGGACCTGCGGGAAGACCTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTg                                                                                                      <  2:390219/144‑1 (MQ=255)

                                                    cAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGaa                                                                                                <  2:2286036/150‑1 (MQ=255)

                                                    cAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGaa                                                                                                <  1:1718727/150‑1 (MQ=255)

                                                      ggaggaCCTGCGGGAAGAACTAGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACa                                                                                              >  2:602339/1‑150 (MQ=255)

                                                      ggaggaCATGCGGGAAAAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAACGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCGACACCCTGTACGAACa                                                                                              <  2:1223853/150‑1 (MQ=255)

                                                       gaggaCCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAg                                                                                             >  1:341888/1‑150 (MQ=255)

                                                           aCCGGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGc                                                                                         <  1:1095505/150‑1 (MQ=255)

                                                            ccTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAac                                                                                                            >  1:557649/1‑130 (MQ=255)

                                                            ccTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAac                                                                                                            <  2:557649/130‑1 (MQ=255)

                                                            ccTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCa                                                                                        >  2:297148/1‑150 (MQ=255)

                                                             cTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAg                                                                                       <  1:526461/150‑1 (MQ=255)

                                                             cTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAg                                                                                       >  1:982084/1‑150 (MQ=255)

                                                                   gaagaaCTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAACGTCGTCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACag                                                                                 >  1:2272764/1‑150 (MQ=255)

                                                                   gaagaaCTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACAGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACag                                                                                 >  2:1769012/1‑150 (MQ=255)

                                                                      gaaCTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACAGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGagc                                                                              <  1:1769012/150‑1 (MQ=255)

                                                                          tGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAgcagca                                                                          <  1:736162/150‑1 (MQ=255)

                                                                          tGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAgcagca                                                                          >  2:1799380/1‑150 (MQ=255)

                                                                          tGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAgcagca                                                                          >  2:705830/1‑150 (MQ=255)

                                                                          tGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCATGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCAGCAACACCCTGTACGAACAGATGCAGGCACAGAgcagca                                                                          >  2:303701/1‑150 (MQ=255)

                                                                             cggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACgg                                                                       <  1:506895/150‑1 (MQ=255)

                                                                              ggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCTGGCAAGcgc                                                                                                                                            >  1:1079348/1‑78 (MQ=255)

                                                                              ggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCTGGCAAGcgc                                                                                                                                            <  2:1079348/80‑3 (MQ=255)

                                                                                   tgctgGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCa                                                                 >  1:517441/1‑150 (MQ=255)

                                                                                      tgGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCAtgg                                                              >  1:501531/1‑148 (MQ=255)

                                                                                       ggAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCgg                                                             >  2:1155232/1‑150 (MQ=255)

                                                                                          aGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCt                                                          <  2:1212227/150‑1 (MQ=255)

                                                                                             cgccAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCCATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTa                                                       >  1:155469/1‑150 (MQ=255)

                                                                                              gccAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTAAGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTAt                                                      >  1:1958551/1‑150 (MQ=255)

                                                                                                cAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGc                                                    >  1:1250299/1‑150 (MQ=255)

                                                                                                cAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGc                                                    <  2:341888/150‑1 (MQ=255)

                                                                                                 aaGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGcc                                                   <  1:611687/150‑1 (MQ=255)

                                                                                                 aaGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGACCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGGATGcc                                                   >  2:1238453/1‑150 (MQ=255)

                                                                                                  atggcttaCCGCCGTCACCGTCATCCTTCGCCCTGCCGGATGCTCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCt                                                  <  2:2229794/143‑1 (MQ=11)

                                                                                                    gCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTTGGCGCATCAACACCCTGTACGAACAGATGTAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTgg                                                <  2:1982824/150‑1 (MQ=255)

                                                                                                    gCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTgg                                                >  2:245300/1‑150 (MQ=40)

                                                                                                      cGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTa                                              <  1:1571000/150‑1 (MQ=39)

                                                                                                       gAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCACGCAAGCCCGTGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTAc                                             >  1:55387/1‑150 (MQ=255)

                                                                                                           cgccgTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTcg                                         >  2:1305146/1‑150 (MQ=38)

                                                                                                           cgccgTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTcg                                         >  1:1968136/1‑150 (MQ=38)

                                                                                                            gccgTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCCTCAACACCCTGTACGACCAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTcgc                                        <  1:2126728/150‑1 (MQ=255)

                                                                                                                 cAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGcc                                                                                                                           <  2:28051/62‑1 (MQ=37)

                                                                                                                 cAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGcc                                                                                                                           >  1:28051/1‑62 (MQ=37)

                                                                                                                  aGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCAACTGCCAAGGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTggg                                  >  2:2432664/1‑150 (MQ=255)

                                                                                                                            ttCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGtg                        >  1:620893/1‑150 (MQ=255)

                                                                                                                               gCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGtgctg                     >  2:1400968/1‑150 (MQ=255)

                                                                                                                                 ccTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTgcg                   >  1:1966744/1‑150 (MQ=255)

                                                                                                                                  cTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTgcgc                  >  2:1850422/1‑150 (MQ=255)

                                                                                                                                  cTGCCGGAAGCGCGTGCCAGGCAAGCCGTTGGCGCTGGCATTACCTTAGGGCGCATCAACACCCTTTACGAGCAGATGCATGCACAGAGCTGCACAGCCAGCATCGGCCTGTATGCCTGGTACCTCGTTCTGGGGCGGCAAGtgctttgc                  >  2:1098653/1‑146 (MQ=11)

                                                                                                                                       ggAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGc             <  1:1305146/150‑1 (MQ=40)

                                                                                                                                       ggAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGc             <  2:1568838/150‑1 (MQ=255)

                                                                                                                                         atgggcTTGCCTGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCa           >  2:1314447/5‑150 (MQ=255)

                                                                                                                                            cgcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCAccc        <  2:155469/150‑1 (MQ=255)

                                                                                                                                             gcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACAcc                                                                                                         >  2:2276754/1‑52 (MQ=38)

                                                                                                                                             gcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACAcc                                                                                                         <  1:2276754/52‑1 (MQ=35)

                                                                                                                                             gcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCAcccc       >  1:1691697/1‑150 (MQ=255)

                                                                                                                                             gcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCAcccc       >  1:698695/1‑150 (MQ=255)

                                                                                                                                               ttGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCa     >  1:2194737/1‑150 (MQ=255)

                                                                                                                                                taccAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTa                                              >  2:2355595/3‑108 (MQ=255)

                                                                                                                                                taccAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTa                                              <  1:2355595/106‑1 (MQ=255)

                                                                                                                                                  ccAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACa  >  1:1314741/1‑150 (MQ=255)

                                                                                                                                                  ccAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACa  >  2:515761/1‑150 (MQ=255)

CAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACA  >  NC_002947/3837879‑3838174

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