BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	733,650	A→T	6.2%	intergenic (+21/‑26)	PP_t09 → / → PP_t10	tRNA‑Lys/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	733,650	0	A	T	6.2%	446.7 / 10.4	149	intergenic (+21/‑26)	PP_t09/PP_t10	tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand): ref base A (77/62); new base T (4/6); total (81/68)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

TTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAA  >  NC_002947/733503‑733799

tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                     <  1:2379638/150‑1 (MQ=255)

  tGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAgaga                                                                                                                                                   <  1:762867/150‑1 (MQ=255)

  tGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAgaga                                                                                                                                                   <  1:762980/150‑1 (MQ=255)

      aGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                               <  1:3295492/150‑1 (MQ=255)

      aGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                               <  1:3295428/150‑1 (MQ=255)

      aGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGcc                                                                                                                                               <  2:1704218/150‑1 (MQ=255)

         aGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCtt                                                                                                                                            <  2:1831442/150‑1 (MQ=255)

           gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttt                                                                                                                                          <  1:1772872/150‑1 (MQ=255)

            aTGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTATCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATTTTGAAAAGGGCATCTCTTTGAGATGCCCttttc                                                                                                                                         <  1:4675381/150‑2 (MQ=17)

            aTGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCttttc                                                                                                                                         <  1:4080020/150‑2 (MQ=255)

               tCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttttt                                                                                                                                      <  2:1891401/150‑1 (MQ=255)

                cAGGCGTGCGAAGCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCCGTGGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGACAAGGGCATCTCAAAGAGATGCCCTTTTTTTTg                                                                                                                                     <  1:1029240/150‑1 (MQ=17)

                 aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGa                                                                                                                                    <  1:1327858/150‑1 (MQ=255)

                   gAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaaa                                                                                                                                  <  2:2352035/150‑1 (MQ=255)

                    aGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAac                                                                                                                                 <  1:4705694/150‑1 (MQ=255)

                    aGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAac                                                                                                                                 <  1:4705567/150‑1 (MQ=255)

                     gTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAAATCTTAGTAGGTTCTAATAAAACACGACCCACAGTTTAGTAAAGCGCAACTCAATGTCATGCCCTTGTTGTTGAAaca                                                                                                                                >  2:3516632/1‑150 (MQ=14)

                      ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttt                                                                                                                                          >  2:1919135/1‑139 (MQ=255)

                        cGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTc                                                                                                                             <  1:1236347/150‑1 (MQ=255)

                         gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                            >  2:4285538/1‑150 (MQ=255)

                         gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                            <  1:1276952/150‑1 (MQ=255)

                          caTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                           <  2:2533131/149‑1 (MQ=255)

                          aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                           <  1:5117666/150‑1 (MQ=255)

                          aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                           <  1:381726/150‑1 (MQ=255)

                           aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTTGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                          <  2:2561272/150‑1 (MQ=255)

                              cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtg                                                                                                                       >  1:205402/1‑150 (MQ=255)

                               ccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtgt                                                                                                                      >  1:2219109/1‑150 (MQ=255)

                                cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                     >  2:1857315/1‑150 (MQ=255)

                                cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGGTGGCTTTTAACCAATTGGTCGTAGGTGCGAATCCTACACGACCCACCAGGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGGGTa                                                                                                                     >  2:3965356/1‑150 (MQ=255)

                                 tccccgagACCATTTTCCGGGTCGTTAGCTCCGTTGGTAGAGAAGTTTGCTTTTATCCAATTTGTAGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCCAAGAGTTGAACTTTTTTTTGAAACACTCGGGGTGTAg                                                                                                                    <  1:1745599/143‑1 (MQ=21)

                                 tCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTAAACGACCCACCAAGTAGAAAAGGGCATCTCAAACAGATGCCCTTTTTTTTTAAACACTCGGGGTGTAg                                                                                                                    >  1:4605756/1‑150 (MQ=255)

                                   aCCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAATGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAgcg                                                                                                                  >  1:5077164/1‑150 (MQ=255)

                                   aCCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAgcg                                                                                                                  <  1:4285538/150‑1 (MQ=255)

                                     cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                                <  1:4324123/150‑1 (MQ=255)

                                        gACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAATAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                             >  1:684397/1‑150 (MQ=255)

                                         aCCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAAGGCCATAGCAAAGAGATGCCCTTAGTTATGAAACAATCGGGATGTAGCGCAGTcg                                                                                                            >  2:2350443/1‑149 (MQ=17)

                                           cATTTTCCGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGGCTTTCAACCCATTGGTCGTAGGTTCGACACGCACCCGAGCCACCCTTTAGAAAAGGGCATCTCCAACTGTTGCACTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCgg                                                                                                          <  2:1721366/150‑1 (MQ=21)

                                                        gTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCt                                                                                             >  2:428770/1‑150 (MQ=255)

                                                           aGCCCAGCTGGTCGAGCAGTTGGCTTTTCCCCAATTGGCCGTAGGATCGAATCCTACACGACCCACCATGTACAACAGGGCACCTCAAAGAGATGCCCTTATCTTTGAAACCCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                          <  2:1041379/150‑1 (MQ=21)

                                                            gCTCCGTTGGTAGAGCAGTTGGCATTTAACCATTTGGTCGTAGGTTCGAATCCTAGACTACCGAGCATGTAGAAAAGGTCATCTCAAATTTTTGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                         <  1:4712282/150‑1 (MQ=25)

                                                            gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                         <  1:2178660/150‑1 (MQ=255)

                                                            gCTCAGTTGGTAGAGCAGTTGGCTTTGAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                         <  2:2733940/150‑1 (MQ=255)

                                                               cAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAg                                                                                      >  1:106665/1‑150 (MQ=255)

                                                                 gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTATTTTTGAAACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                    <  1:2775066/150‑1 (MQ=255)

                                                                  ttGCTAGAGCAGATGGCATTAAACCAATTGGTCGTAGCTTCGAATCCTACAAGACACACCATGTAGAAAAGGTCATCACAAAGACATGACCATTTTTATGAAAAACTCGGTGTGTAGAGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                   <  2:1113391/150‑1 (MQ=25)

                                                                        gagCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTc                                                                             >  2:2331408/1‑150 (MQ=255)

                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGtt                                                                      >  1:2824631/1‑150 (MQ=255)

                                                                                 gCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATAACAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGAAGCGCAGTCCGGTAGCGCGCCAGCTTTGGGAGCAGGATGTCAGGAGTTCg                                                                    >  2:831819/1‑150 (MQ=255)

                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTCGCGCGCCAGCTTTGGGAGCAGGAGGTCAGGCGTTCGaa                                                                  >  2:2875279/1‑150 (MQ=255)

                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                  >  1:3644553/1‑150 (MQ=255)

                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTACAAAAGGGCATCACAAAGAGATGCCCATTTATTTGAAACACACGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATccc                                                              >  2:3893473/1‑150 (MQ=37)

                                                                                        aCCAATTGGTCGTAGGTGCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTTAAACACTAGGGGTGTAGCCCAGTCCGGTAGCGCGCCTGCTTGGGGAGCAGGATGTCAGGAGTTCGAATcccc                                                             >  1:352919/1‑150 (MQ=37)

                                                                                          cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAATGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTGGGGAGCAGGATGACAGGAGTTCTAATccccac                                                           >  2:5569559/1‑148 (MQ=255)

                                                                                          cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                           >  2:1269831/1‑150 (MQ=255)

                                                                                           aaTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCa                                                          >  1:10483/1‑150 (MQ=255)

                                                                                            aTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAGGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAc                                                         >  1:9657/1‑150 (MQ=255)

                                                                                            aTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAc                                                         >  1:206131/1‑150 (MQ=255)

                                                                                            aTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAc                                                         >  1:3021255/1‑150 (MQ=255)

                                                                                            aTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAc                                                         >  1:4496326/1‑150 (MQ=255)

                                                                                            aTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAc                                                         >  1:9598/1‑150 (MQ=255)

                                                                                                 tCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGa                                                    >  2:1595052/1‑150 (MQ=255)

                                                                                                   gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                       <  2:4516750/49‑1 (MQ=255)

                                                                                                   gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                       >  1:4516750/1‑49 (MQ=255)

                                                                                                      ggTTCGAATCCTACACGACCCACCATGTAGAAATGGCCATCACAAAGAGATGCCATTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCGCTCGCCCCGGccggt                                               >  2:690198/1‑147 (MQ=25)

                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTa                                                                                                        <  2:1196245/91‑1 (MQ=255)

                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTa                                                                                                        >  1:1196245/1‑91 (MQ=255)

                                                                                                            aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTggg                                         >  1:1077150/1‑150 (MQ=255)

                                                                                                            aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTggg                                         >  1:633229/1‑150 (MQ=255)

                                                                                                                cTACACGACCCACCATGTAGACAAGGGCATCTCACAGAGATGCCCTTTTTTTTGAACCACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                     <  1:3407305/150‑1 (MQ=255)

                                                                                                                cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                     >  2:3223456/1‑150 (MQ=255)

                                                                                                                 tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAGTTTTGGGTCGtt                                    >  2:5397588/1‑150 (MQ=255)

                                                                                                                 gacacGCCCCTCCCTGCAGACTTGGGCATCTCCCAGAGATGCCCTTATTTTAGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGtt                                    <  1:1641495/149‑1 (MQ=17)

                                                                                                                  acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                             <  1:4865990/76‑1 (MQ=255)

                                                                                                                  acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                             >  2:4865990/1‑76 (MQ=255)

                                                                                                                  acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTa                                   >  2:4021393/1‑150 (MQ=255)

                                                                                                                  acacGACCCACCATGTAGAAAAGGCCAACACAAAGAGATGCCCTTTTTATTTAAACACTCGGGGTGAAGCGCAGTCCGGAACCGCGCCTGCTTACGGAGCAGGATGTCAGGAGTACGAATCCCCTCACCCCGACCATGTATGGGTCGTTa                                   >  2:4021373/1‑150 (MQ=21)

                                                                                                                    acGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGc                                 <  1:1869856/150‑1 (MQ=255)

                                                                                                                    acGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGc                                 <  1:1482000/150‑1 (MQ=255)

                                                                                                                        cccACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAg                             <  2:205402/150‑1 (MQ=255)

                                                                                                                         ccaccaTGTAGCAAAGGGCATCTCAAAGAGATGCCCAATTTTTTGAATCTCTCGGGGTGTAGCGCAGTGCGGTAGCGCTCCTGCTTTGTGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGt                            <  2:692993/150‑1 (MQ=18)

                                                                                                                         ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGt                            >  2:3478524/1‑150 (MQ=255)

                                                                                                                          caccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGtt                           <  2:496127/150‑1 (MQ=255)

                                                                                                                            ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                         <  1:4082718/150‑1 (MQ=255)

                                                                                                                             caTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                        <  1:4498306/150‑1 (MQ=255)

                                                                                                                             caTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                        >  1:907371/1‑150 (MQ=255)

                                                                                                                                 tAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTagag                    <  2:106665/150‑1 (MQ=255)

                                                                                                                                 aaaaaaaGGGCATCTCAAAGAGATGCCCTTTTTTTTTAAAAAAGCGGGGTGTAGCGCAGTCCGGTAGCGCGCCAGCTTTGGGAGCAGGATGTCAGGAGTTCGAATACCCTCACCCCGACCAATTTTTGGTCGTTATCAAAGTTGGTagag                    >  2:140352/4‑150 (MQ=18)

                                                                                                                                    aaaaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAg                 >  2:2112653/1‑150 (MQ=255)

                                                                                                                                       aGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg              >  1:1241245/1‑150 (MQ=255)

                                                                                                                                       aGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg              >  2:5237892/1‑150 (MQ=255)

                                                                                                                                       aGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATGTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg              >  2:1282808/1‑150 (MQ=255)

                                                                                                                                         ggCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGc            >  1:5000857/1‑150 (MQ=255)

                                                                                                                                         ggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTGGGTAGAGCAGTTGGc            >  1:5199559/1‑150 (MQ=255)

                                                                                                                                          gCATCTCAAAGAGATGCCCTTTTTTTAGAAACACTCGCGGAGTAGCGCAGTCCGGTAGCGCGCCTGCGTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCACACCCCGACAATTTATGGGTCGGTAGCTCAGGTGGTAGAGCAGTAGGCt           >  2:5358140/1‑150 (MQ=18)

                                                                                                                                           cATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt          >  1:3507226/1‑150 (MQ=255)

                                                                                                                                               tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa      >  1:1139969/1‑150 (MQ=255)

                                                                                                                                               tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa      >  1:3019487/1‑150 (MQ=255)

                                                                                                                                               tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa      >  2:516589/1‑150 (MQ=255)

                                                                                                                                                cAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAc     >  2:546052/1‑150 (MQ=255)

                                                                                                                                                   agagACGCCATTTATTTTTTTCCACTCTGTGTGTAGCGCTGTCCGGTCGAGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCaa  <  1:797904/150‑1 (MQ=11)

TTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAA  >  NC_002947/733503‑733799

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