BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,154,429	T→G	6.7%	Q129P (CAG→CCG)	PP_1909 ←	septum formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,154,429	0	T	G	6.7%	335.5 / 15.5	120	Q129P (CAG→CCG)	PP_1909	septum formation protein
Reads supporting (aligned to +/- strand): ref base T (64/48); new base G (4/4); total (68/52)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 121 total aligned reads displayed.

ACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGA  >  NC_002947/2154280‑2154577

aCCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCt                                                                                                                                                      <  2:568598/150‑1 (MQ=255)

  ccGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTgg                                                                                                                                                    >  1:1389088/1‑150 (MQ=255)

            tctGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCgg                                                                                                                                          >  1:1349355/1‑150 (MQ=255)

             ctGGAACAAGCTCACCCACAGCCCCTCCGCCTTGGAACTCCCTTCTCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGt                                                                                                                                         <  2:843229/150‑1 (MQ=255)

              tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                        >  2:1277768/1‑150 (MQ=255)

              tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                        >  1:340901/1‑150 (MQ=255)

              tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                        >  2:562093/1‑150 (MQ=255)

               ggAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTgg                                                                                                                                       <  1:562093/150‑1 (MQ=255)

                 aaCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                     >  2:242378/1‑150 (MQ=255)

                 aaCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCCCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                     >  2:109484/1‑150 (MQ=255)

                   ctaGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCGCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCt                                                                                                                                   <  1:1081423/148‑1 (MQ=255)

                      gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTTCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                         >  2:231769/1‑141 (MQ=255)

                      gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTTCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                         <  1:231769/141‑1 (MQ=255)

                      gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCCTCGCCAGCTGACGCTCtt                                                                                                                                >  2:1081423/1‑150 (MQ=255)

                      gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                         <  1:231659/141‑1 (MQ=255)

                      gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                         >  2:231659/1‑141 (MQ=255)

                        tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTca                                                                                                                              >  2:1951818/1‑150 (MQ=255)

                        tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTca                                                                                                                              >  1:2209683/1‑150 (MQ=255)

                         cACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCCAGGTCAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                                                             <  1:540957/150‑1 (MQ=255)

                           cccccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcagca                                                                                                                           <  1:1153960/150‑1 (MQ=255)

                             cccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                                                         <  2:2209683/150‑1 (MQ=255)

                               taGCCCCTCCTCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                                                                       <  1:1257991/149‑1 (MQ=255)

                                 gCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCa                                                                                                                     <  1:16119/150‑1 (MQ=255)

                                  ccccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCCGCTCCCTGCTCAGTTCCCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCACTGCCCGGTGGCACTGCTGAACAGCGCCac                                                                                                                    <  1:683820/150‑2 (MQ=255)

                                    ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGcc                                                                                                                  >  1:907956/1‑150 (MQ=255)

                                    ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGcc                                                                                                                  <  2:1389088/150‑1 (MQ=255)

                                    ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCGCGGTTAATGCGCTGGCCCTGCAGGGGGTGGGCAGGCAACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGcc                                                                                                                  >  1:860718/1‑150 (MQ=255)

                                     cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGccc                                                                                                                 >  2:1015432/1‑150 (MQ=255)

                                       ccGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCCCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                               <  1:728948/150‑1 (MQ=255)

                                       ccGCCTTGAAACTCCCTGCGCAATACGGCTTCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                               <  1:1170462/150‑1 (MQ=255)

                                         gCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTca                                                                                                                              <  1:2293896/133‑1 (MQ=255)

                                            ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                                                                       <  1:399533/137‑1 (MQ=255)

                                            ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                                                                       >  2:399533/1‑137 (MQ=255)

                                              gAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                        >  2:1568216/1‑150 (MQ=255)

                                                    cccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCAc                                                                                                  <  2:1849732/150‑1 (MQ=255)

                                                           cAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGCCCGGTAAGAAAGCTCACCtgctgc                                                                                           >  1:278349/1‑150 (MQ=255)

                                                           cAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCtgctgc                                                                                           <  1:1957322/150‑1 (MQ=255)

                                                                 aGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTg                                                                                     <  2:340901/150‑1 (MQ=255)

                                                                     gCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACt                                                                                 >  1:2257959/1‑150 (MQ=255)

                                                                       tGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTcc                                                                               >  1:271804/1‑150 (MQ=255)

                                                                        gTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                              <  1:1951818/150‑1 (MQ=255)

                                                                        gTTCTGCTGCCACGTAACGCCCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                              <  1:923993/150‑1 (MQ=255)

                                                                            tgctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                >  1:1207437/1‑144 (MQ=255)

                                                                            tgctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                <  2:1207437/144‑1 (MQ=255)

                                                                              ctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGTAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                        >  1:2223090/1‑150 (MQ=255)

                                                                               tgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                       <  1:109484/150‑1 (MQ=255)

                                                                               tgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                       >  2:376812/1‑150 (MQ=255)

                                                                                  cACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                                                                    <  1:242378/150‑1 (MQ=255)

                                                                                     gTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                                                         >  2:834411/1‑94 (MQ=255)

                                                                                      tAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCaa                                                                <  1:1015432/150‑1 (MQ=255)

                                                                                          gCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAgcgc                                                            >  2:110702/1‑150 (MQ=255)

                                                                                              gACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTc                                                        >  1:1417454/1‑150 (MQ=255)

                                                                                              gACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTc                                                        <  1:1568216/150‑1 (MQ=255)

                                                                                                cctgctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGa                                                      >  1:1041539/1‑150 (MQ=255)

                                                                                                    ctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGt                                                  >  1:2200274/1‑150 (MQ=255)

                                                                                                     tccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                        <  2:1406295/95‑1 (MQ=255)

                                                                                                     tccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                        >  1:1406295/1‑95 (MQ=255)

                                                                                                      ccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTAt                                                >  1:608815/1‑150 (MQ=255)

                                                                                                       cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                               >  1:1597062/1‑150 (MQ=255)

                                                                                                       cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                               >  2:1218859/1‑150 (MQ=255)

                                                                                                       cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                               >  2:396644/1‑150 (MQ=255)

                                                                                                         tgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCg                                             >  1:1648092/1‑150 (MQ=255)

                                                                                                           ctcAGTTCCCGCCGTGCCACGGCAAATGACTCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCACTATTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGc                                           <  1:514330/150‑1 (MQ=255)

                                                                                                           ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGc                                           >  1:1750148/1‑150 (MQ=255)

                                                                                                           ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCCCGAACGTATGCGGc                                           >  2:54916/1‑150 (MQ=255)

                                                                                                            tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                        <  1:747408/120‑1 (MQ=255)

                                                                                                            tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                        >  2:747408/1‑120 (MQ=255)

                                                                                                            tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCt                                          >  1:1402399/1‑150 (MQ=255)

                                                                                                            tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCt                                          >  2:1069376/1‑150 (MQ=255)

                                                                                                             cAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCtt                                         <  1:1218859/150‑1 (MQ=255)

                                                                                                             cAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCtt                                         <  1:396644/150‑1 (MQ=255)

                                                                                                               gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGc                                       <  2:271804/150‑1 (MQ=255)

                                                                                                                ttCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCGCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGcc                                      <  2:2257959/150‑1 (MQ=255)

                                                                                                                ttCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGcc                                      >  2:452116/1‑150 (MQ=255)

                                                                                                                  cACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCa                                    <  2:1402399/150‑1 (MQ=255)

                                                                                                                   aCGCAGTGTCACGGTCACAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAg                                   >  2:208852/1‑150 (MQ=255)

                                                                                                                   aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCCCTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAg                                   >  2:1699000/1‑150 (MQ=255)

                                                                                                                   aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAg                                   >  2:1410950/1‑150 (MQ=255)

                                                                                                                       aGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAcc                               >  1:1006075/1‑150 (MQ=255)

                                                                                                                       aGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAcc                               <  1:1620406/150‑1 (MQ=255)

                                                                                                                          gtCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGc                            <  2:1417454/150‑1 (MQ=255)

                                                                                                                            cACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTc                          <  1:376812/150‑1 (MQ=255)

                                                                                                                               ggTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCACGGGCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGcc                       >  1:651202/1‑150 (MQ=255)

                                                                                                                                        cACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                              <  2:491870/86‑1 (MQ=255)

                                                                                                                                        cACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                              >  1:491870/1‑86 (MQ=255)

                                                                                                                                        cACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGc              <  2:1692316/150‑1 (MQ=255)

                                                                                                                                         aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCa             >  2:2305209/1‑150 (MQ=255)

                                                                                                                                         aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCa             >  1:2023868/1‑150 (MQ=255)

                                                                                                                                         aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGACCACAACAGTCCAGTAGCTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCa             >  2:1134920/1‑150 (MQ=255)

                                                                                                                                          cGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAc                                >  2:940521/1‑130 (MQ=255)

                                                                                                                                          cGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAc                                <  1:940521/130‑1 (MQ=255)

                                                                                                                                             aaTCCACCTGGCATTTCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                                                                      >  2:338562/1‑73 (MQ=255)

                                                                                                                                             aaTCCACCTGGCATTTCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                                                                      <  1:338562/73‑1 (MQ=255)

                                                                                                                                             aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt         >  2:2279702/1‑150 (MQ=255)

                                                                                                                                             aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt         >  2:177244/1‑150 (MQ=255)

                                                                                                                                             aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt         >  2:1109974/1‑150 (MQ=255)

                                                                                                                                              aTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTg        <  2:2223090/150‑1 (MQ=255)

                                                                                                                                               gccACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                                                                    >  2:1520763/2‑89 (MQ=255)

                                                                                                                                                  aCCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTcc    >  2:1758612/1‑150 (MQ=255)

                                                                                                                                                    cTGGCATTGCCCGGTGGCGCGGTTCAGCAGTGCCATCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGa  <  2:1643997/150‑1 (MQ=255)

ACCCGTGCGTGCTCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGA  >  NC_002947/2154280‑2154577

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