BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,292,341	G→A	5.7%	Q24* (CAG→TAG)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,292,341	0	G	A	5.7%	485.7 / 11.8	140	Q24* (CAG→TAG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand): ref base G (80/52); new base A (4/4); total (84/56)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 141 total aligned reads displayed.

ACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAT  >  NC_000913/2292192‑2292489

acCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTg                                                                                                                                                      <  1:1820435/150‑1 (MQ=255)

 ccTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGc                                                                                                                                                     <  1:1634208/150‑1 (MQ=255)

      cTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAacac                                                                                                                                                >  2:1145806/1‑150 (MQ=255)

        gCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACAt                                                                                                                                              <  1:1810477/150‑1 (MQ=255)

        gCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACAt                                                                                                                                              <  1:1065473/150‑1 (MQ=255)

           tCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGa                                                                                                                                           >  2:401920/1‑150 (MQ=255)

           tCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGa                                                                                                                                           >  2:401913/1‑150 (MQ=255)

            cAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGAt                                                                                                                                          <  2:1168979/150‑1 (MQ=255)

              gTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATgg                                                                                                                                        >  1:1364060/1‑150 (MQ=255)

              gTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATgg                                                                                                                                        >  1:1290180/1‑150 (MQ=255)

              gTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATgg                                                                                                                                        >  1:1163352/1‑150 (MQ=255)

               ttCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGt                                                                                                                                       >  2:1126197/1‑150 (MQ=255)

                  aTACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGc                                                                                                                                    >  2:825218/1‑150 (MQ=255)

                   tACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAATCTGCAACACATCGATGGTCGcc                                                                                                                                   >  2:2194254/1‑150 (MQ=255)

                      acTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                                                                                                >  1:2332113/1‑150 (MQ=255)

                       cTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcg                                                                                                                               <  2:230822/150‑1 (MQ=255)

                       cTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcg                                                                                                                               <  2:89591/150‑1 (MQ=255)

                        tttCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcgc                                                                                                                              >  1:2375997/1‑150 (MQ=255)

                        tttCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcgc                                                                                                                              >  1:2375959/1‑150 (MQ=255)

                        tttCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcgc                                                                                                                              >  1:2345453/1‑150 (MQ=255)

                           cTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGATCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGa                                                                                                                           <  2:979039/150‑1 (MQ=255)

                             gACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGc                                                                                                                         <  1:924092/150‑1 (MQ=255)

                               cGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCg                                                                                                                       <  2:1196341/150‑1 (MQ=255)

                                gCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCgg                                                                                                                      <  1:2245587/150‑1 (MQ=255)

                                gCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCgg                                                                                                                      >  1:813039/1‑150 (MQ=255)

                                 cAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGga                                                                                                                     <  2:176434/150‑1 (MQ=255)

                                 cAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGga                                                                                                                     <  1:456193/150‑1 (MQ=255)

                                 cAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGga                                                                                                                     <  1:456120/150‑1 (MQ=255)

                                 cAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGga                                                                                                                     <  1:456119/150‑1 (MQ=255)

                                  aGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgag                                                                                                                    >  2:746650/1‑150 (MQ=255)

                                   ggTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                                                                                                   >  2:1604613/1‑150 (MQ=255)

                                        gTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatcatc                                                                                                              <  1:308188/150‑1 (MQ=255)

                                             aaTCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcat                                                                                                         >  1:1632542/1‑150 (MQ=255)

                                             aaTCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcat                                                                                                         >  2:887043/1‑150 (MQ=255)

                                                catcGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatcag                                                                                                      >  1:1869911/1‑150 (MQ=255)

                                                  tcGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCAtcagcatcagca                                                                                                    <  2:1106551/150‑1 (MQ=255)

                                                    gAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACg                                                                                                  >  2:213499/1‑150 (MQ=255)

                                                    gAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACg                                                                                                  >  2:213381/1‑150 (MQ=255)

                                                     aGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCATTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGc                                                                                                 <  1:1361865/150‑1 (MQ=255)

                                                       ttGGTATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGccc                                                                                               <  2:2035231/150‑1 (MQ=255)

                                                       ttGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGccc                                                                                               <  2:2262258/150‑1 (MQ=255)

                                                       ttGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGccc                                                                                               <  2:495360/150‑1 (MQ=255)

                                                          gAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAAt                                                                                            >  1:1820109/1‑150 (MQ=255)

                                                          gAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAAt                                                                                            >  2:47504/1‑150 (MQ=255)

                                                             tCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                                                                                         >  2:2311775/1‑150 (MQ=255)

                                                             tCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                                                                                         >  1:1514582/1‑150 (MQ=255)

                                                             tCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                                                                                         >  1:2244464/1‑150 (MQ=255)

                                                             tCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                                                                                         >  2:292118/1‑150 (MQ=255)

                                                             tCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                                                                                         >  2:185943/1‑150 (MQ=255)

                                                             tCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                                                                                         >  2:808183/1‑150 (MQ=255)

                                                               ggCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                                                                                       >  1:2260725/1‑150 (MQ=255)

                                                               ggCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                                                                                       >  2:1044338/1‑150 (MQ=255)

                                                               ggCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                                                                                       >  1:2260469/1‑150 (MQ=255)

                                                               ggCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                                                                                       <  1:1586866/150‑1 (MQ=255)

                                                                  aaTGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                                                                                    >  2:611853/1‑150 (MQ=255)

                                                                               tGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcct                                                                       >  2:253328/1‑150 (MQ=255)

                                                                                 gCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctcc                                                                     >  1:314466/1‑150 (MQ=255)

                                                                                    tttCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTa                                                                  <  1:1440785/150‑1 (MQ=255)

                                                                                      tCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACt                                                                <  2:684757/150‑1 (MQ=255)

                                                                                         ggCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTa                                                             >  2:644313/1‑150 (MQ=255)

                                                                                         ggCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTa                                                             <  1:2149801/150‑1 (MQ=255)

                                                                                         ggCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTa                                                             <  1:1178187/150‑1 (MQ=255)

                                                                                          gCAGCGCAGTTCTTCAGTGAGCTGTCGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTAt                                                            >  1:1806301/1‑150 (MQ=255)

                                                                                          gCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTAt                                                            <  2:2215153/150‑1 (MQ=255)

                                                                                               gcAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTtctc                                                       <  1:705439/150‑1 (MQ=255)

                                                                                                       tcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAATTCCTCATTGTGCGGCCTCGTTTCTGTATTTCTCCCACAGCg                                               >  1:259679/1‑150 (MQ=255)

                                                                                                         aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGc                                             <  1:634696/150‑1 (MQ=255)

                                                                                                         aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGc                                             <  1:1126197/150‑1 (MQ=255)

                                                                                                         aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGc                                             <  2:653899/150‑1 (MQ=255)

                                                                                                           tGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCtt                                           <  1:746650/150‑1 (MQ=255)

                                                                                                               cTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGAtc                                       <  2:2126624/150‑1 (MQ=255)

                                                                                                               cTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGAtc                                       <  2:972603/150‑1 (MQ=255)

                                                                                                                          gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCa                            >  2:838241/1‑150 (MQ=255)

                                                                                                                          gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCa                            >  1:1076001/1‑150 (MQ=255)

                                                                                                                          gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCa                            >  1:1674083/1‑150 (MQ=255)

                                                                                                                            tgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGa                          <  1:1883442/150‑1 (MQ=255)

                                                                                                                                 cTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGa                     <  2:341709/150‑1 (MQ=255)

                                                                                                                                      cTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGAtt                <  1:1229535/150‑1 (MQ=255)

                                                                                                                                       ttCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTc               <  2:2092647/150‑1 (MQ=255)

                                                                                                                                        tCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCa              >  1:446061/1‑150 (MQ=255)

                                                                                                                                         cATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAg             >  1:1212353/1‑150 (MQ=255)

                                                                                                                                           tCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGAt           >  2:52765/1‑150 (MQ=255)

                                                                                                                                           tCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGAt           >  2:112868/1‑150 (MQ=255)

                                                                                                                                             tttAAACTACAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATtgtg                                                                             <  2:576785/82‑1 (MQ=255)

                                                                                                                                             tttAAACTACAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATtgtg                                                                             <  2:576946/82‑1 (MQ=255)

                                                                                                                                             tttAAACTACAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATtgtg                                                                             >  1:576785/1‑82 (MQ=255)

                                                                                                                                             tttAAACTACAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATtgtg                                                                             >  1:576946/1‑82 (MQ=255)

                                                                                                                                             tttAAACTACAACACATCGATGGTCGCCAGCGCCGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATtgtg                                                                             <  2:576965/82‑1 (MQ=255)

                                                                                                                                             tttAAACTACAACACATCGATGGTCGCCAGCGCAGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATtgtg                                                                             >  1:576965/1‑82 (MQ=255)

                                                                                                                                               tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGcc       <  1:470433/150‑1 (MQ=255)

                                                                                                                                               tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGcc       <  1:611853/150‑1 (MQ=255)

                                                                                                                                                 aaCTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCg     >  1:2340760/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  2:713558/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  1:987617/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  2:1160259/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  2:1202204/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  2:1238635/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  1:1772718/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  1:177807/1‑150 (MQ=255)

                                                                                                                                                    tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAt  >  2:1309636/1‑150 (MQ=255)

ACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGCAATGCTGTTGTTCTGGCATTTCGGGCAGCGCAGTTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTTACTGTATTTCTCCCACAGCGGCTTGATCTCTTCTTCCCAGACGCGAGGATTCAGATCGCCCGCAT  >  NC_000913/2292192‑2292489

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