BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,776,798	T→G	8.0%	T353T (ACT→ACG)	betA‑II →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,776,798	0	T	G	8.0%	263.6 / 14.6	100	T353T (ACT→ACG)	betA‑II	choline dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (42/50); new base G (5/3); total (47/53)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

Only 100 of 101 total aligned reads displayed.

CTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAA  >  NC_002947/5776649‑5776944

cTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACt                                                                                                                                                     >  1:1784590/1‑150 (MQ=255)

    aCCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTATGGTGGAACCAGCCGGCGATCTGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTc                                                                                                                                                 <  1:877371/150‑1 (MQ=255)

     ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTcc                                                                                                                                                <  2:2237415/150‑1 (MQ=255)

      cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCg                                                                                                                                               >  2:2221712/1‑150 (MQ=255)

       tGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCgg                                                                                                                                              >  2:4848/1‑150 (MQ=255)

        gCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGa                                                                                                                                             >  1:2233127/1‑150 (MQ=255)

         cAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACGAGTCCGGat                                                                                                                                            <  1:2290651/150‑2 (MQ=255)

              cGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       <  2:1351016/150‑1 (MQ=255)

               gCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAg                                                                                                                                      <  1:2425740/150‑1 (MQ=255)

                   gCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGc                                                                                                                                  <  2:2012983/150‑1 (MQ=255)

                        cAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAAc                                                                                                                             >  1:239912/1‑150 (MQ=255)

                             ggTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCa                                                                                                                        >  1:584760/1‑150 (MQ=255)

                             ggTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGAGCCTAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCa                                                                                                                        >  2:1908091/1‑150 (MQ=255)

                              gTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       >  1:972987/1‑150 (MQ=255)

                                         cccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCt                                                                                                            <  1:4848/150‑1 (MQ=255)

                                         cccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCt                                                                                                            >  1:1883543/1‑150 (MQ=255)

                                         cccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCt                                                                                                            >  1:1516754/1‑150 (MQ=255)

                                          ccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTg                                                                                                           <  2:239912/150‑1 (MQ=255)

                                           cGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAAAGATCGGAAGGCCGGGGGGGTGATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGc                                                                                                          >  1:2252328/1‑150 (MQ=255)

                                            gTCGCTGCGGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGcc                                                                                                         >  1:2038811/1‑150 (MQ=255)

                                                ctgctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGtt                                                                                                     <  2:2164213/150‑1 (MQ=255)

                                                  gctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGc                                                                                                   <  2:1688273/150‑1 (MQ=255)

                                                    tgtgGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                                                                 <  1:590429/150‑1 (MQ=255)

                                                      tggtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAtt                                                                                               >  1:1335900/1‑150 (MQ=255)

                                                        gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTa                                                                                                                     <  1:818100/126‑1 (MQ=255)

                                                        gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTa                                                                                                                     >  2:818100/1‑126 (MQ=255)

                                                        gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTc                                                                                                              >  1:753011/1‑133 (MQ=255)

                                                        gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTc                                                                                                              <  2:753011/133‑1 (MQ=255)

                                                        gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTaa                                                                                             >  2:2388996/1‑150 (MQ=255)

                                                        gtggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTCCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTaa                                                                                             <  1:553312/150‑1 (MQ=255)

                                                         tggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                                                            >  1:494171/1‑150 (MQ=255)

                                                         tggAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                                                            <  2:741819/150‑1 (MQ=255)

                                                              ccAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACaa                                                                                       >  1:1889129/1‑150 (MQ=255)

                                                              ccAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACaa                                                                                       <  2:1889129/150‑1 (MQ=255)

                                                              ccAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCC‑GGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAAc                                                                                      <  1:689978/150‑1 (MQ=255)

                                                               cAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTANAAc                                                                                      <  2:1516754/150‑1 (MQ=255)

                                                                aGCCGGCGATCGGTGCCGAGAGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGTCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACg                                                                                     >  2:742652/1‑150 (MQ=255)

                                                                   cGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCt                                                                                  >  1:1524669/1‑150 (MQ=255)

                                                                    ggCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTc                                                                                 <  1:239890/150‑1 (MQ=255)

                                                                       gATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCaa                                                                              >  2:2330696/1‑150 (MQ=255)

                                                                          cGGTGCCGAGCGGCAGTCCCATCGCACGGGGATCATCGCCAGCAAGCAGTTCTAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATgg                                                                           <  2:856756/150‑1 (MQ=255)

                                                                           ggTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGc                                                                          <  2:1524669/150‑1 (MQ=255)

                                                                            gTGCGGAGTGGTTGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCg                                                                         <  1:1216313/150‑1 (MQ=255)

                                                                             tGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGt                                                                        <  2:1614786/150‑1 (MQ=255)

                                                                                cGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       >  2:1491281/1‑100 (MQ=255)

                                                                                cGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       <  1:1491281/100‑1 (MQ=255)

                                                                                cGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     <  2:2038811/150‑1 (MQ=255)

                                                                                 gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     >  1:816454/1‑149 (MQ=255)

                                                                                 gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     <  2:816454/149‑1 (MQ=255)

                                                                                 gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    >  2:18129/1‑150 (MQ=255)

                                                                                 gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    <  1:2449777/150‑1 (MQ=255)

                                                                                  aGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGTCGGGGGGGTTCATGCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAgg                                                                   >  2:929957/1‑150 (MQ=255)

                                                                                  aGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATCAACTACAACGGCTCCAATGGCGTGAAgg                                                                   <  1:423861/150‑1 (MQ=255)

                                                                                     ggATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAAc                                                                <  2:1784590/150‑1 (MQ=255)

                                                                                      gATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACa                                                               >  1:2006933/1‑150 (MQ=255)

                                                                                        tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             >  1:1996756/1‑150 (MQ=255)

                                                                                        tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             <  1:1285113/150‑1 (MQ=255)

                                                                                        tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             <  2:1883543/150‑1 (MQ=255)

                                                                                          ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGc                                                           <  1:67418/150‑1 (MQ=255)

                                                                                             aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       >  2:2155317/1‑71 (MQ=255)

                                                                                             aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       <  1:2155317/71‑1 (MQ=255)

                                                                                             aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc                                                        >  1:1006839/1‑150 (MQ=255)

                                                                                               tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCa                                                      >  1:1032988/1‑150 (MQ=255)

                                                                                               tGGCGCCGATACCGGTGCCAGCAACCAGATCGAGGCCGGGGGGGGCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACAACCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCa                                                      >  2:588741/1‑150 (MQ=255)

                                                                                                ggCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCCCTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAg                                                     <  2:2252328/150‑1 (MQ=255)

                                                                                                ggCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAg                                                     <  1:2221712/150‑1 (MQ=255)

                                                                                                 gCACCGGTATCGGCGCCAGCAACCAGTTCGAGCCCGGGGGGTTCATCCGGACTAGTTCGTAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAgg                                                    <  1:1206570/150‑1 (MQ=255)

                                                                                                  cACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc                                                   >  2:1852082/1‑150 (MQ=255)

                                                                                                  cACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc                                                   <  1:1852082/150‑1 (MQ=255)

                                                                                                   aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGAATAATCCGGACTGGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgcg                                                  >  1:1005103/1‑150 (MQ=255)

                                                                                                         aTCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATg                                            <  2:2006933/150‑1 (MQ=255)

                                                                                                          tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                                                                 >  2:956590/1‑96 (MQ=255)

                                                                                                          tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                                                                 <  1:956590/96‑1 (MQ=255)

                                                                                                             gtgcCAGCAACCGGTTCGGGGCCGGGGGGTTCATCCGGGCTCGTCCGGAGTTCAAGTGACCGGACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACGACGGCTCGAATGGCGTGAAGGAACATGGCTTCCAGTCGCACATG‑GGCt                                       <  2:1032988/148‑1 (MQ=255)

                                                                                                             gcgcCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGaa                                                                     >  1:999219/1‑121 (MQ=255)

                                                                                                             gcgcCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCt                                       <  2:972987/150‑1 (MQ=255)

                                                                                                                 cAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCAt                                   >  2:468666/1‑150 (MQ=255)

                                                                                                                 cAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCAt                                   <  2:1006839/150‑1 (MQ=255)

                                                                                                                     aaCCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATgcgc                               >  2:508220/1‑150 (MQ=255)

                                                                                                                       ccAGTTCGAGGCAAGGGGGGTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTc                             >  2:1798323/1‑150 (MQ=255)

                                                                                                                         aGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGc                           >  1:2436820/1‑150 (MQ=255)

                                                                                                                              gAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    >  2:302322/1‑105 (MQ=255)

                                                                                                                              gAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAg                                                                    <  1:302322/105‑1 (MQ=255)

                                                                                                                                     gggggTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGcc               >  2:2324827/1‑150 (MQ=255)

                                                                                                                                          ttCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATc          <  2:2436820/150‑1 (MQ=255)

                                                                                                                                            actCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCa                                    >  1:709643/3‑122 (MQ=255)

                                                                                                                                            actCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCa                                    <  2:709643/120‑1 (MQ=255)

                                                                                                                                            actCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCa        >  2:1767917/3‑150 (MQ=255)

                                                                                                                                             ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       >  1:1277414/2‑150 (MQ=255)

                                                                                                                                             ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGC‑GATGAGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       >  1:1308061/2‑150 (MQ=255)

                                                                                                                                             aTCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       <  2:1861348/150‑1 (MQ=255)

                                                                                                                                             aTCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       <  2:689484/150‑1 (MQ=255)

                                                                                                                                             aTCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg       <  1:18129/150‑1 (MQ=255)

                                                                                                                                              tCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGc                           >  2:1708348/1‑129 (MQ=255)

                                                                                                                                              tCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGc                           <  1:1708348/129‑1 (MQ=255)

                                                                                                                                                 ggACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCa   >  2:1827976/1‑150 (MQ=255)

                                                                                                                                                  gACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCaa  <  1:1798323/150‑1 (MQ=255)

                                                                                                                                                  gACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCaa  <  1:2324827/150‑1 (MQ=255)

                                                                                                                                                   aCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTAc                                                                                                                    >  1:2442880/1‑36 (MQ=255)

                                                                                                                                                   aCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTAc                                                                                                                    <  2:2442880/36‑1 (MQ=255)

CTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG‑GGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAA  >  NC_002947/5776649‑5776944

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