BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,251,181	C→G	5.4%	S29S (TCG→TCC)	alr ←	alanine racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,251,181	0	C	G	5.4%	1316.9 / 44.2	387	S29S (TCG→TCC)	alr	alanine racemase
Reads supporting (aligned to +/- strand): ref base C (166/200); new base G (18/3); total (184/203)
Fisher's exact test for biased strand distribution p-value = 4.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 398 total aligned reads displayed.

TAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGA  >  NC_002947/4251032‑4251330

tAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCATGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                                       <  2:8138739/150‑1 (MQ=255)

tAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                                       <  2:5018526/150‑1 (MQ=255)

 aGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCg                                                                                                                                                      >  2:1511023/1‑150 (MQ=255)

 aGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCAACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCg                                                                                                                                                      <  1:2235197/150‑1 (MQ=255)

  ggCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGa                                                                                                                                                     <  2:7496096/150‑1 (MQ=255)

  ggCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGa                                                                                                                                                     <  2:4866721/150‑1 (MQ=255)

   gCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGAc                                                                                                                                                    <  1:656433/150‑1 (MQ=255)

   gCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGAc                                                                                                                                                    >  1:7415710/1‑150 (MQ=255)

     aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACaa                                                                                                                                                  >  1:2124509/1‑150 (MQ=255)

     aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACaa                                                                                                                                                  >  2:3456208/1‑150 (MQ=255)

     aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACaa                                                                                                                                                  >  1:7211935/1‑150 (MQ=255)

      tCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAAc                                                                                                                                                 >  2:4070210/1‑150 (MQ=255)

         gCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggt                                                                                                                                              >  2:2768765/1‑150 (MQ=255)

         gCCTTGAGCACGGCGCACAGCTTGGACTGGCCGGCCAGCTCGGCCTGCTGCGTGCGGATGTTGTGCTGCAGGGCTCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggt                                                                                                                                              <  1:7719235/150‑1 (MQ=255)

          ccTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCAGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggtg                                                                                                                                             >  2:8673402/1‑150 (MQ=255)

           cTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggtgg                                                                                                                                            >  2:6244370/1‑150 (MQ=255)

           cTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggtgg                                                                                                                                            >  2:4441709/1‑150 (MQ=255)

            ttGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACggtggt                                                                                                                                           >  1:5786249/1‑150 (MQ=255)

               aGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGcc                                                                                                                                        <  2:5694913/150‑1 (MQ=255)

                gCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCg                                                                                                                                       >  2:1976676/1‑150 (MQ=255)

                 cACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGc                                                                                                                                      <  2:2297548/150‑1 (MQ=255)

                    ggCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATa                                                                                                                                   >  1:1830190/1‑150 (MQ=255)

                     gcgcACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATAc                                                                                                                                  >  2:430256/1‑150 (MQ=255)

                         acaGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                              <  1:5914857/150‑1 (MQ=255)

                           aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            <  1:5670910/150‑1 (MQ=255)

                           aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            >  1:3033735/1‑150 (MQ=255)

                            gCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAgggggg                                                                                                                           <  1:6749368/150‑1 (MQ=255)

                            gCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAgggggg                                                                                                                           <  1:6838654/150‑1 (MQ=255)

                             cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                          <  1:6260191/150‑1 (MQ=255)

                                  aCTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTc                                                                                                                     <  2:2572228/150‑1 (MQ=255)

                                    ttGCCGGCCAGCTCGGCCTGCAGCGTGGGCATGTAGTGATGCAGGGCGCTGGAGCTGACTTAGACCCAGGCATTGCTGTTTTGCACGTTCATGGTGGTGGTGCCGATGTCCATCGACAACGGTGGTTCCGCATACAGGGGGGCCTGTCCg                                                                                                                   >  1:7971619/1‑150 (MQ=255)

                                    ttGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCTGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCg                                                                                                                   <  1:1203544/150‑1 (MQ=255)

                                      gCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGt                                                                                                                 <  1:4632233/150‑1 (MQ=255)

                                      gCCGGCCAACTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGt                                                                                                                 <  1:7034607/150‑1 (MQ=255)

                                       ccGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTg                                                                                                                <  2:6600263/150‑1 (MQ=255)

                                         ggCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGAt                                                                                                              >  2:6315803/1‑150 (MQ=255)

                                         ggCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGAt                                                                                                              <  2:1136570/150‑1 (MQ=255)

                                            cAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAg                                                                                                           >  2:554938/1‑150 (MQ=255)

                                             aGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCAGCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGa                                                                                                          >  2:7872766/1‑150 (MQ=255)

                                              gCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGaa                                                                                                         <  1:8595577/150‑1 (MQ=255)

                                                tCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGt                                                                                                       >  2:4658695/1‑150 (MQ=255)

                                                        cAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGa                                                                                               <  1:3825783/150‑1 (MQ=255)

                                                               cGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGcc                                                                                        >  2:1771985/1‑150 (MQ=255)

                                                                ggATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCa                                                                                       >  2:8770000/1‑150 (MQ=255)

                                                                 gATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAg                                                                                      >  1:4981435/1‑150 (MQ=255)

                                                                   tgttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                                    <  1:5437037/150‑1 (MQ=255)

                                                                   tgttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                                    <  2:7662550/150‑1 (MQ=255)

                                                                   tgttgtGCTGCAGGGCGCTGGCGCTGACCTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                                    <  1:1386414/150‑1 (MQ=255)

                                                                      tgtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAggg                                                                                 >  2:1307530/1‑150 (MQ=255)

                                                                           tgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACgg                                                                            >  2:4154952/1‑150 (MQ=255)

                                                                             cAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCg                                                                          >  1:6172455/1‑150 (MQ=255)

                                                                                ggcgctggcgctGACTTCGACCCCGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGaaa                                                                       <  1:144679/150‑1 (MQ=255)

                                                                                ggcgctggcgctGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGaaa                                                                       <  2:4733487/150‑1 (MQ=255)

                                                                                ggcgctggcgctGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGaaa                                                                       >  2:2123430/1‑150 (MQ=255)

                                                                                   gctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAggg                                                                    >  2:6903539/1‑150 (MQ=255)

                                                                                      ggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCAt                                                                 <  2:2485793/150‑1 (MQ=255)

                                                                                      ggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGATACGGCGAAAGGGCAt                                                                 <  2:5186361/150‑1 (MQ=255)

                                                                                        cgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATgg                                                               <  2:5275790/150‑1 (MQ=255)

                                                                                        cgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATgg                                                               <  2:1904573/150‑1 (MQ=255)

                                                                                               ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTAc                                                        <  2:4939983/150‑1 (MQ=255)

                                                                                               ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTAc                                                        >  2:7803239/1‑150 (MQ=255)

                                                                                               ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTAc                                                        >  2:7394028/1‑150 (MQ=255)

                                                                                                 cGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                                      <  1:8203286/150‑1 (MQ=255)

                                                                                                 cGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                                      <  1:948917/150‑1 (MQ=255)

                                                                                                     ccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgtgt                                                  >  1:8030508/1‑150 (MQ=255)

                                                                                                      cAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGtt                                                 >  1:4459941/1‑150 (MQ=255)

                                                                                                       aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTc                                                >  1:3395301/1‑150 (MQ=255)

                                                                                                       aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTc                                                <  1:4896394/150‑1 (MQ=255)

                                                                                                       aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTc                                                <  2:7002352/150‑1 (MQ=255)

                                                                                                         gCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCtt                                              <  2:7210391/150‑1 (MQ=255)

                                                                                                            ttGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTc                                           <  1:8767109/150‑1 (MQ=255)

                                                                                                            ttGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTc                                           >  2:2662002/1‑150 (MQ=255)

                                                                                                               cTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                        <  1:5130954/150‑1 (MQ=255)

                                                                                                               cTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                        <  2:5767763/150‑1 (MQ=255)

                                                                                                                tGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAg                                       >  1:2325911/1‑150 (MQ=255)

                                                                                                                tGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAg                                       >  1:6957020/1‑150 (MQ=255)

                                                                                                                     tGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACACCGGTGGTGCCGCATACCGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTTCGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTgcgc                                  <  1:347087/150‑1 (MQ=255)

                                                                                                                            tCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCTGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                           <  1:6193126/150‑1 (MQ=255)

                                                                                                                            tCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                           >  2:8254938/1‑150 (MQ=255)

                                                                                                                            tCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                           <  2:2209417/150‑1 (MQ=255)

                                                                                                                            tCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                           >  2:1661964/1‑150 (MQ=255)

                                                                                                                            tCAGGGGGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTCCGGCGAAAGGGCATGGGGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTaa                           >  2:2431703/1‑150 (MQ=255)

                                                                                                                             cAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                        >  2:3011148/1‑136 (MQ=255)

                                                                                                                             cAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                        <  1:3011148/136‑1 (MQ=255)

                                                                                                                                ggTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGa                       >  1:8408003/1‑150 (MQ=255)

                                                                                                                                    gcaccaccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCa                   >  2:3147219/7‑150 (MQ=255)

                                                                                                                                     tGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg                  <  2:4510857/150‑1 (MQ=255)

                                                                                                                                      ggTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAgg                 >  1:1978603/1‑150 (MQ=255)

                                                                                                                                       gTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt                <  2:5942978/150‑1 (MQ=255)

                                                                                                                                        caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg               >  1:2483255/3‑150 (MQ=255)

                                                                                                                                        caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg               >  1:5587453/3‑150 (MQ=255)

                                                                                                                                        caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg               >  1:1077020/3‑150 (MQ=255)

                                                                                                                                        caccGTTGTCGATGGACAACGGCGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg               >  1:3222926/3‑150 (MQ=255)

                                                                                                                                         gCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtgt              >  1:5299395/1‑150 (MQ=255)

                                                                                                                                         accGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtgt              >  2:2862855/2‑150 (MQ=255)

                                                                                                                                             ttGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGttt          >  2:3859542/1‑150 (MQ=255)

                                                                                                                                               gTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGa        >  2:2533836/1‑150 (MQ=255)

                                                                                                                                                 ccATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGt      <  1:6472937/150‑1 (MQ=255)

                                                                                                                                                  cATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGtt     >  2:852025/1‑150 (MQ=255)

                                                                                                                                                     cGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGa  <  1:7616006/150‑1 (MQ=255)

TAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGA  >  NC_002947/4251032‑4251330

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