BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,971,711	T→A	9.5%	N521Y (AAC→TAC)	PP_5236 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,971,711	0	T	A	9.5%	490.7 / 32.7	165	N521Y (AAC→TAC)	PP_5236	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (71/77); new base A (7/9); total (78/87)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 167 total aligned reads displayed.

AACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGG  >  NC_002947/5971563‑5971860

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

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

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

      cGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTg                                                                                                                                                <  1:535994/150‑1 (MQ=255)

        tAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgc                                                                                                                                              <  2:3699488/150‑1 (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)

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

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

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

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

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

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

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

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

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

             aaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcc                                                                                                                                         <  1:770090/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)

                   ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGAAGGCCAAGTACAACCTGgcctgc                                                                                                                                          >  1:3839628/1‑143 (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                                                                                                                                                    >  1:209679/1‑115 (MQ=255)

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

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

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

                                        cATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCa                                                                                                              <  1:5182547/150‑1 (MQ=255)

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

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

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

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

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

                                                    cGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCgg                                                                                                  <  1:2973756/150‑1 (MQ=255)

                                                     gTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGc                                                                                                 >  2:4269338/1‑150 (MQ=255)

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

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

                                                           cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCTTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                           <  2:468549/150‑1 (MQ=255)

                                                           cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCTGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                           <  1:4097169/150‑1 (MQ=255)

                                                           cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                           >  1:4945505/1‑150 (MQ=255)

                                                           cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                           <  1:889584/150‑1 (MQ=255)

                                                             gcgcGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGt                                                                                         <  1:3356062/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)

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

                                                                       aaGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcct                                                                                                                                        >  1:2736853/1‑93 (MQ=255)

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

                                                                                ccTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgg                                                                      <  2:1270972/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)

                                                                                                    tGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCAt                                                  <  2:2378278/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:4694747/1‑150 (MQ=255)

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

                                                                                                            ttCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGc                                          <  2:2610732/150‑1 (MQ=255)

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

                                                                                                                  tCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATAc                                    >  1:4845842/1‑150 (MQ=255)

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

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

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

                                                                                                                     tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTg                                              >  1:2057099/1‑137 (MQ=255)

                                                                                                                     tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTg                                              <  2:2057099/137‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                        tCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCCGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGcc                              >  2:216608/1‑150 (MQ=255)

                                                                                                                             gTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgc                         >  1:4770221/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)

                                                                                                                                      caggcaggCCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCgg                >  2:3348290/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)

                                                                                                                                            ggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGtt          <  1:4393614/150‑1 (MQ=255)

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

                                                                                                                                                   ttGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCg   <  2:4913804/150‑1 (MQ=255)

AACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGG  >  NC_002947/5971563‑5971860

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