BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	962,276:1	+A	5.3%	coding (142/1035 nt)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	962,276	1	.	A	5.3%	463.6 / 18.9	113	coding (142/1035 nt)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand): ref base . (52/55); new base A (2/4); total (54/59)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 116 total aligned reads displayed.

CGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGCAC  >  NC_002947/962132‑962423

cgccCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGGCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATg                                                                                                                                                  >  2:260127/1‑150 (MQ=255)

  cccATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCg                                                                                                                                                >  2:746914/1‑150 (MQ=255)

  cccATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCG‑CGCGACGATGcc                                                                                                                                                <  1:3005109/150‑2 (MQ=255)

   ccATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCG‑CGCGACGATGcag                                                                                                                                               <  1:1283951/150‑3 (MQ=255)

         cATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGt                                                                                                                                               <  2:2407909/144‑1 (MQ=255)

          aTAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGAc                                                                                                                                        <  1:1974573/150‑1 (MQ=255)

          aTAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGAc                                                                                                                                        >  1:1824566/1‑150 (MQ=255)

           tAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCTAcc                                                                                                                                       <  1:42866/150‑1 (MQ=255)

               tGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTTCCGACCcaca                                                                                                                                   >  2:1246750/1‑150 (MQ=255)

                     tGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCg                                                                                                                             <  1:3000042/150‑1 (MQ=255)

                      gATGGCGTTGGTGGCGGCCATGCGGCCCCGGCCGTCCGGGTGGATGCCGTCATGCATGAAGAAGGCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            <  2:674851/150‑1 (MQ=255)

                      gATGGCGTTGGTGGCGGCAATGCGGCGGAAGTAGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            <  2:2489823/150‑1 (MQ=255)

                      gATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            <  1:2737924/150‑1 (MQ=255)

                       aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCggg                                                                                                                           <  2:2847540/150‑1 (MQ=255)

                       aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCggg                                                                                                                           <  2:2912584/150‑1 (MQ=255)

                            gTTGGTGGCGGCAATGCGGCCGAAGTCATAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCG‑CGCGACGATGca                                                                                                                                                <  1:1070528/124‑2 (MQ=255)

                             ttGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                                                     <  1:260127/150‑1 (MQ=255)

                             ttGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                                                     <  1:214701/150‑1 (MQ=255)

                             ttGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                                                     >  2:1685021/1‑150 (MQ=255)

                                gtggCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGcc                                                                                                                  >  1:1667536/1‑150 (MQ=255)

                                  ggcggcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTAGGGTCCAGGGCCTg                                                                                                                <  2:559559/150‑1 (MQ=255)

                                     ggcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGt                                                                                                             <  2:1824566/150‑1 (MQ=255)

                                      gcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTg                                                                                                            >  1:588231/1‑150 (MQ=255)

                                        aaTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGcc                                                                                                          <  2:1080925/150‑1 (MQ=255)

                                         aTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCg                                                                                                         <  1:542384/150‑1 (MQ=255)

                                                 gAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGACCTgg                                                                                                               <  1:2327571/136‑1 (MQ=255)

                                                 gAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGACCTgg                                                                                                               >  2:2327571/1‑136 (MQ=255)

                                                  aaGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCa                                                                                                >  1:1292380/1‑150 (MQ=255)

                                                    gTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAg                                                                                              <  1:2166954/150‑1 (MQ=255)

                                                     tCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAgg                                                                                             >  1:2576620/1‑150 (MQ=255)

                                                      cGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGa                                                                                            <  2:1303183/150‑1 (MQ=255)

                                                       gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTTCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGAt                                                                                           >  2:914671/1‑150 (MQ=255)

                                                       gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGAt                                                                                           >  1:2897251/1‑150 (MQ=255)

                                                          aGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                                                        <  2:1667536/150‑1 (MQ=255)

                                                          aGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCG‑CGCGACGATGcagca                                                                                                                                             >  2:2253515/1‑93 (MQ=255)

                                                           gggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACg                                                                                       >  2:1133757/1‑150 (MQ=255)

                                                            ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            <  2:435807/112‑1 (MQ=255)

                                                            ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            >  1:435807/1‑112 (MQ=255)

                                                            ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                      >  1:236766/1‑150 (MQ=255)

                                                            ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                      <  2:1206500/150‑1 (MQ=255)

                                                              tCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTcc                                                                                    >  2:1162145/1‑150 (MQ=255)

                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                                                               >  2:247897/1‑118 (MQ=255)

                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                                                               <  1:247897/118‑1 (MQ=255)

                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                                                        <  2:2759920/141‑1 (MQ=255)

                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                      <  2:2240294/143‑1 (MQ=255)

                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGt                                                                                      >  1:2240294/1‑143 (MQ=255)

                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCAc                                                                               >  1:776677/1‑150 (MQ=255)

                                                                      gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCt                                                                                                                 <  1:2377645/113‑1 (MQ=255)

                                                                      gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCt                                                                                                                 >  2:2377645/1‑113 (MQ=255)

                                                                      gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCg                                                                                                   <  1:2619989/127‑1 (MQ=255)

                                                                      gaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCtt                                                                            >  2:868404/1‑150 (MQ=255)

                                                                       aTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                                                           <  1:868404/150‑1 (MQ=255)

                                                                       aTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                                                           <  2:90795/150‑1 (MQ=255)

                                                                        tCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑g                                                                         <  2:222841/150‑1 (MQ=255)

                                                                          ggCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑ggg                                                                       <  2:1454789/150‑1 (MQ=255)

                                                                           gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑gggg                                                                      >  2:1564489/1‑150 (MQ=255)

                                                                              tGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑g                                                                         <  2:11608/144‑1 (MQ=255)

                                                                              tGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑g                                                                         >  1:11608/1‑144 (MQ=255)

                                                                               gaagaaGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GAt                                                                 <  2:1285602/150‑1 (MQ=255)

                                                                                aagaagTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATc                                                                <  2:2897251/150‑1 (MQ=255)

                                                                                  gaagTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCtt                                                              <  1:1564489/150‑1 (MQ=255)

                                                                                     gTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            >  2:2822891/1‑87 (MQ=255)

                                                                                     gTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCgg                                                                                                                            <  1:2822891/87‑1 (MQ=255)

                                                                                     gTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCCTG‑GGGGA‑GATCTTGCa                                                           >  1:1668912/1‑150 (MQ=255)

                                                                                      tCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGAATCTTG‑GGGGA‑GATCTTGCaa                                                          <  2:2576620/150‑1 (MQ=255)

                                                                                       cGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCAACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAAc                                                         <  2:513204/150‑1 (MQ=255)

                                                                                             tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTggg                                                                                                              >  2:1043658/1‑93 (MQ=255)

                                                                                             tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTggg                                                                                                              <  1:1043658/93‑1 (MQ=255)

                                                                                             tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGc                                                   >  1:2639262/1‑150 (MQ=255)

                                                                                                     aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCa                                           <  2:236766/150‑1 (MQ=255)

                                                                                                        ccgcGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATgg                                        <  1:1526437/150‑1 (MQ=255)

                                                                                                         cgcgGTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCT‑CGATGCTTTGCCGATCCACTGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATGACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGc                                       <  1:2082584/150‑1 (MQ=255)

                                                                                                            gcTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGc                                                   >  1:2344271/1‑135 (MQ=255)

                                                                                                                gtcgtcCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑ggg                                                                       >  1:71987/1‑112 (MQ=255)

                                                                                                                gtcgtcCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg‑ggg                                                                       <  2:71987/112‑1 (MQ=255)

                                                                                                                 tcgtcCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTgg                               >  1:983745/1‑150 (MQ=255)

                                                                                                                 tcgtcCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTgg                               <  2:1635518/150‑1 (MQ=255)

                                                                                                                      aaGCGCCTACACTGCTTCATAGTCGCG‑CGATGCGTTGCCGACCCACGTTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAg                          <  1:1427655/149‑1 (MQ=255)

                                                                                                                        gcgcGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGc                        <  2:3068981/150‑1 (MQ=255)

                                                                                                                              acacTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCg                  <  1:1162145/150‑1 (MQ=255)

                                                                                                                              acacTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCg                  <  2:3038234/150‑1 (MQ=255)

                                                                                                                              acacTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCg                  >  1:2997869/1‑150 (MQ=255)

                                                                                                                                acTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgc                >  1:2322343/1‑150 (MQ=255)

                                                                                                                                  tgctgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCaa                                                          <  2:2753932/106‑1 (MQ=255)

                                                                                                                                  tgctgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCaa                                                          >  1:2753932/1‑106 (MQ=255)

                                                                                                                                  tgctgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGc                        <  1:1518958/140‑1 (MQ=255)

                                                                                                                                  tgctgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctc              >  1:511018/1‑150 (MQ=255)

                                                                                                                                  tgctgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctc              >  2:1768594/1‑150 (MQ=255)

                                                                                                                                  tgctgcATCGTCGCG‑AGCTGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGc                        >  2:1518958/1‑140 (MQ=255)

                                                                                                                                    ctgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCCTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctcgc            <  2:1668912/150‑1 (MQ=255)

                                                                                                                                     tgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTg                                                                                                                <  1:2262999/51‑1 (MQ=255)

                                                                                                                                     tgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTg                                                                                                                >  2:2262999/1‑51 (MQ=255)

                                                                                                                                     tgcATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGgctcgct           >  2:2329382/1‑150 (MQ=255)

                                                                                                                                       cATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCa         <  1:914671/150‑1 (MQ=255)

                                                                                                                                        aTCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGT‑‑AACGCTTTGCAGGCATTAACTGGCAACTACAGCCGGCGCCCCGGCTGTCTTAGGTCAGGCGAATCGGCTCGCTCAt        >  1:1595585/1‑150 (MQ=255)

                                                                                                                                           gtcgCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACg     >  1:425632/1‑150 (MQ=255)

                                                                                                                                            tcgCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGc    >  1:2823314/1‑150 (MQ=255)

                                                                                                                                            tcgCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGc    >  2:1733106/1‑150 (MQ=255)

                                                                                                                                              gcg‑cgATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGGCGCCTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGCCGGTACTGGGTCAGGCGAATCGGCTCACTCATACGCAc  <  2:1143044/150‑1 (MQ=255)

                                                                                                                                              gcg‑cgATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTcc                                                                                    <  2:304782/70‑1 (MQ=255)

CGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCG‑CGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTG‑GGGGA‑GATCTTGCAACCACAGCCGGCGCCATGGCTGTACTGGGTCAGGCGAATCGGCTCGCTCATACGCAC  >  NC_002947/962132‑962423

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