BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,971,708	G→A	10.2%	L522L (CTG→TTG)	PP_5236 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,971,708	0	G	A	10.2%	465.3 / 32.2	166	L522L (CTG→TTG)	PP_5236	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (72/77); new base A (7/10); total (79/87)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

ACCAAACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAG  >  NC_002947/5971559‑5971855

 ccAAACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAgg                                                                                                                                                    >  1:2693444/1‑150 (MQ=255)

    aaCCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTg                                                                                                                                                 <  2:83554/150‑1 (MQ=255)

     cccGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGt                                                                                                                                                <  1:880155/149‑1 (MQ=255)

       cGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTAc                                                                                                                                              >  2:3683535/1‑150 (MQ=255)

         gCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACtt                                                                                                                                            >  2:1342139/1‑150 (MQ=255)

         gCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACtt                                                                                                                                            <  1:3896706/150‑1 (MQ=255)

           gTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                          >  1:4103881/1‑150 (MQ=255)

           gTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                          <  1:2767608/150‑1 (MQ=255)

            tAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgc                                                                                                                                         <  2:3699488/150‑1 (MQ=255)

            tAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGtacaac                                                                                                                                              >  1:3663252/1‑140 (MQ=255)

             aaCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcc                                                                                                                                        >  1:4526770/1‑150 (MQ=255)

             aaCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcc                                                                                                                                        <  1:2500132/150‑1 (MQ=255)

              aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                       <  1:1270129/150‑1 (MQ=255)

              aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                       <  1:2198535/150‑1 (MQ=255)

              aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCGGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                       <  2:5339554/150‑1 (MQ=255)

               cAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAgg                                                                                                                                                    >  1:3373669/1‑136 (MQ=255)

               cAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctg                                                                                                                                      <  1:576728/150‑1 (MQ=255)

               cAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctg                                                                                                                                      >  1:2969579/1‑150 (MQ=255)

               cAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctg                                                                                                                                      >  2:2648455/1‑150 (MQ=255)

                aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                     <  1:5256918/150‑1 (MQ=255)

                aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                     >  1:3338656/1‑150 (MQ=255)

                aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                     <  2:246012/150‑1 (MQ=255)

                aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                     <  2:245913/150‑1 (MQ=255)

                aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcctgc                                                                                                                                     <  2:5289309/150‑1 (MQ=255)

                 aaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcc                                                                                                                                    <  1:769907/150‑1 (MQ=255)

                  atgtcccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcct                                                                                                                                   <  1:158660/146‑1 (MQ=255)

                    tCCCAAGGTGCAACCGCGCACAGGCAGTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACAGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCCGGCCAGGTTGTACTTGGCctgcctgc                                                                                                                                 >  2:4992512/1‑150 (MQ=255)

                       ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCtgctgc                                                                                                                              >  1:1664965/1‑150 (MQ=255)

                       ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCtgctgc                                                                                                                              >  1:1229074/1‑150 (MQ=255)

                       ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTGCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGGTGTACTTGGCCTGCCtgctgc                                                                                                                              >  2:4469812/1‑150 (MQ=255)

                           gTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACACGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGc                                                                                                                          <  2:2391405/150‑1 (MQ=255)

                               aaCCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCtgc                                                                                                                      >  2:5298487/1‑148 (MQ=255)

                               aaCCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTcc                                                                                                                      >  1:1270972/1‑150 (MQ=255)

                                  cgcgCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCAcc                                                                                                                   <  2:2997173/150‑1 (MQ=255)

                                       acacGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACg                                                                                                              <  2:1433466/150‑1 (MQ=255)

                                       acacGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACg                                                                                                              >  2:2116876/1‑150 (MQ=255)

                                         acGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTa                                                                                                                                               <  2:209679/115‑1 (MQ=255)

                                         acGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTa                                                                                                                                               >  1:209679/1‑115 (MQ=255)

                                         acGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCt                                                                                                            >  1:1024081/1‑150 (MQ=255)

                                           gCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGc                                                                                                          >  1:3932787/1‑150 (MQ=255)

                                           gCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGc                                                                                                          >  1:5129586/1‑150 (MQ=255)

                                             aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACTGCACGCTGCAg                                                                                                        >  2:1615286/1‑150 (MQ=255)

                                                ggTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                     <  1:5399063/150‑1 (MQ=255)

                                                 gTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGAGCAGGCAGGCCAAGTACAACCTGgcctgc                                                                                                                                     >  1:5115156/1‑117 (MQ=38)

                                                 gTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGAGCAGGCAGGCCAAGTACAACCTGgcctgc                                                                                                                                     <  2:5115156/117‑1 (MQ=38)

                                                    tGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAc                                                                                                 >  2:3410240/1‑150 (MQ=255)

                                                     gCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAcc                                                                                                >  2:4097169/1‑150 (MQ=255)

                                                     gCCCGGTGCCCAGCGCGTCGACACGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAcc                                                                                                <  2:3217230/150‑1 (MQ=255)

                                                              ccAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCt                                                                                       >  1:2844861/1‑150 (MQ=255)

                                                               cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCTTGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                      <  1:3013681/150‑1 (MQ=255)

                                                                aGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCg                                                                                     <  2:4343842/150‑1 (MQ=255)

                                                                     gTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATc                                                                                <  1:4887604/150‑1 (MQ=255)

                                                                        gACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAg                                                                             >  1:4295136/1‑150 (MQ=255)

                                                                         aCAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAgg                                                                            <  2:2724132/150‑1 (MQ=255)

                                                                         aCAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAgg                                                                            <  1:4067264/150‑1 (MQ=255)

                                                                           aaGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTa                                                                          >  1:5352852/1‑150 (MQ=255)

                                                                             gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                        >  1:5368213/1‑150 (MQ=255)

                                                                             gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                        >  1:4913804/1‑150 (MQ=255)

                                                                                    ccTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgg                                                                 <  2:1270972/150‑1 (MQ=255)

                                                                                     ctcaCCGCCAGCATGGCACTGCGCACGTTCAGATCCTCCACGCCGTCCAGGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                <  1:208873/147‑1 (MQ=255)

                                                                                     cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                <  2:2302813/150‑1 (MQ=255)

                                                                                     cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                <  1:4964167/150‑1 (MQ=255)

                                                                                           gCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTg                                                          <  1:1103434/150‑1 (MQ=255)

                                                                                                 aTGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAgg                                                    <  1:878673/150‑1 (MQ=255)

                                                                                                      aCTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgcc                                               >  2:3108626/1‑150 (MQ=255)

                                                                                                        tGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCAt                                             >  2:3756883/1‑150 (MQ=255)

                                                                                                        tGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCAt                                             <  1:81657/150‑1 (MQ=255)

                                                                                                        tGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCAt                                             >  2:2214085/1‑150 (MQ=255)

                                                                                                        tGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCAt                                             >  2:2214270/1‑150 (MQ=255)

                                                                                                          cgcTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATgg                                           <  1:4269338/150‑1 (MQ=255)

                                                                                                            cTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTg                                         >  2:4106946/1‑150 (MQ=255)

                                                                                                              cGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCa                                       <  1:1840444/150‑1 (MQ=255)

                                                                                                               gTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAg                                      <  2:2995576/150‑1 (MQ=255)

                                                                                                               gTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAg                                      >  1:4640976/1‑150 (MQ=255)

                                                                                                               gTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAg                                      <  2:3795607/150‑1 (MQ=255)

                                                                                                                ttCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGc                                     >  1:3694666/1‑150 (MQ=255)

                                                                                                                   aGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCa                                  <  1:4727427/150‑1 (MQ=255)

                                                                                                                        cTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACaa                             <  2:1080732/150‑1 (MQ=255)

                                                                                                                         tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATcc                                                                               >  2:1007244/1‑99 (MQ=255)

                                                                                                                         tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATcc                                                                               <  1:1007244/99‑1 (MQ=255)

                                                                                                                         tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                            >  2:4016779/1‑150 (MQ=255)

                                                                                                                         tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                            <  2:4644621/150‑1 (MQ=255)

                                                                                                                         tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                            >  1:4644621/1‑150 (MQ=255)

                                                                                                                          cATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAgg                           >  1:15983/1‑150 (MQ=255)

                                                                                                                           aTCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGc                          <  2:599629/150‑1 (MQ=255)

                                                                                                                                 gTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgc                    >  1:4770459/1‑150 (MQ=255)

                                                                                                                                 gTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgc                    <  2:5352852/150‑1 (MQ=255)

                                                                                                                                  gccATGCGCAGGCAGGCCAGGTGGGACGTGGCCTGCCTGCTGCCCGCCTCCCCCGCACGCTGCAGCCCGGACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgca                   <  2:4837962/149‑1 (MQ=255)

                                                                                                                                     aTGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCCTACAAGGCCAGGgcagca                >  2:4010535/1‑150 (MQ=255)

                                                                                                                                       gcgcAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCAcc              <  2:1229074/150‑1 (MQ=255)

                                                                                                                                        cgcAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCAccc             >  2:4813494/1‑150 (MQ=255)

                                                                                                                                          caggcaggCCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCgg           >  2:3348290/1‑150 (MQ=255)

                                                                                                                                          caggcaggCCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCgg           >  1:2982395/1‑150 (MQ=255)

                                                                                                                                            ggcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCa         <  2:1664965/150‑1 (MQ=255)

                                                                                                                                             gcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGGACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg        <  1:3683535/150‑1 (MQ=255)

                                                                                                                                             gcaggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCa                                  >  2:856930/1‑124 (MQ=255)

                                                                                                                                             gcaggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg        >  2:1845911/1‑150 (MQ=255)

                                                                                                                                              caggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAgg       >  1:2368631/1‑150 (MQ=255)

                                                                                                                                                 gcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTc    <  1:2260771/150‑1 (MQ=255)

                                                                                                                                                  ccAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCa   >  1:4784939/1‑150 (MQ=255)

ACCAAACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAG  >  NC_002947/5971559‑5971855

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