BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,974,739	G→C	7.0%	Y520* (TAC→TAG) ‡	PP_5236 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,974,739	0	G	C	7.0%	995.4 / 47.8	302	Y520* (TAC→TAG) ‡	PP_5236	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (134/146); new base C (5/16); total (139/163)
Fisher's exact test for biased strand distribution p-value = 4.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 307 total aligned reads displayed.

AACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGC  >  NC_002947/5974590‑5974888

aaCCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTg                                                                                                                                                       >  2:1087514/1‑150 (MQ=255)

 aCCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGt                                                                                                                                                      <  2:1274780/150‑1 (MQ=255)

  ccGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTTGTa                                                                                                                                                     <  2:35227/150‑1 (MQ=255)

    ggCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACt                                                                                                                                                   >  1:8323531/1‑150 (MQ=255)

    ggCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACt                                                                                                                                                   >  2:7172805/1‑150 (MQ=255)

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

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

          aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcct                                                                                                                                             <  2:2719713/150‑1 (MQ=255)

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

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

            aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGAGCAGGCAGGCCAAGTACAACCTGgcctgc                                                                                                                                           <  1:4045826/150‑1 (MQ=255)

            aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCAGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                           >  2:3152883/1‑150 (MQ=255)

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

             aaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcc                                                                                                                                          >  2:5475254/1‑150 (MQ=255)

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

             aaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcc                                                                                                                                          >  2:563609/1‑150 (MQ=255)

              aaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcct                                                                                                                                         <  2:2670720/150‑1 (MQ=255)

              aaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcctgcct                                                                                                                                         <  2:5503864/150‑1 (MQ=255)

                 ccccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCctgct                                                                                                                                      <  1:5927227/150‑1 (MQ=255)

                  cccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCctgctg                                                                                                                                     <  2:173226/150‑1 (MQ=255)

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

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

                   ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGGCctgcctgcgca                                                                                                                                    <  1:6972290/150‑4 (MQ=39)

                     aGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGccc                                                                                                                                  <  2:6232150/150‑1 (MQ=255)

                     aGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGccc                                                                                                                                  <  1:684589/150‑1 (MQ=255)

                         gCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCt                                                                                                                              <  2:2274995/150‑1 (MQ=255)

                         gCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCt                                                                                                                              >  1:6182375/1‑150 (MQ=255)

                          cAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTc                                                                                                                             >  2:5519438/1‑150 (MQ=255)

                          cAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTc                                                                                                                             >  2:2213064/1‑150 (MQ=255)

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

                                       gCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCTTGGCACTGCGCTGGTTCAGATCCTCATCGCCGTCCATGCGCTGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGc                                                                                                                <  1:8439691/150‑1 (MQ=255)

                                          ttGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGc                                                                                                             <  1:6559711/150‑1 (MQ=255)

                                            ggTCTGCCCGTTGCCCAGCGCGTCGACACGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                           <  1:5157159/150‑1 (MQ=255)

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

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

                                                 gCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAcc                                                                                                      <  1:6658067/150‑1 (MQ=255)

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

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

                                                          ccAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCt                                                                                             <  2:5299291/150‑1 (MQ=255)

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

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

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

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

                                                           cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCATCATCGCCGTCCATGCGCAGGAAGGCAAGTTTGTAATTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                            <  1:2385678/150‑1 (MQ=255)

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

                                                             gcgcGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGt                                                                                          <  1:1304626/150‑1 (MQ=255)

                                                                cgTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGtgcg                                                                                       <  1:178167/150‑4 (MQ=255)

                                                                    gACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAg                                                                                   <  1:4918759/150‑1 (MQ=255)

                                                                     aCAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAgg                                                                                  >  1:5500541/1‑150 (MQ=255)

                                                                     aCAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGACTGCCTCGTCATCCAgg                                                                                  >  2:7378826/1‑150 (MQ=255)

                                                                       aaGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCTGGTTGTACTTGGCGTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTa                                                                                <  1:4603228/150‑1 (MQ=255)

                                                                        cgCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTaa                                                                               <  1:2682539/149‑1 (MQ=255)

                                                                        aGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTaa                                                                               <  1:3801772/150‑1 (MQ=255)

                                                                           aGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATgc                                                                            >  1:2000524/1‑150 (MQ=255)

                                                                                ccTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTCCTGCCCACCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgg                                                                       <  2:1949399/150‑1 (MQ=255)

                                                                                 cTGACCGCCAGCTTGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCTGGCAGGCCAGTTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                      <  1:7375215/150‑1 (MQ=255)

                                                                                 cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                      >  2:6993798/1‑150 (MQ=255)

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

                                                                                 cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                      >  1:1376636/1‑150 (MQ=255)

                                                                                  tGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgggg                                                                     >  1:4004647/1‑150 (MQ=255)

                                                                                  tGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgggg                                                                     >  2:6259013/1‑150 (MQ=255)

                                                                                  tGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgggg                                                                     >  2:3951502/1‑150 (MQ=255)

                                                                                  tGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgggg                                                                     >  2:851162/1‑150 (MQ=255)

                                                                                     ccGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGgtgt                                                                  <  2:4699997/150‑1 (MQ=255)

                                                                                      cGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGtt                                                                 >  2:4166247/1‑150 (MQ=255)

                                                                                      cGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGtt                                                                 >  1:1317105/1‑150 (MQ=255)

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

                                                                                       gCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTg                                                                >  2:2825019/1‑150 (MQ=255)

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

                                                                                       gCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGCGTTg                                                                <  1:4281150/150‑1 (MQ=255)

                                                                                         cAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGc                                                              >  1:6380374/1‑150 (MQ=255)

                                                                                           gCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCa                                                            <  1:7981921/150‑1 (MQ=255)

                                                                                            cATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAg                                                           <  2:1370820/150‑1 (MQ=255)

                                                                                              tGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgc                                                         <  1:2125978/150‑1 (MQ=255)

                                                                                                 cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                      <  2:259870/150‑1 (MQ=255)

                                                                                                 cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                      <  1:5029206/150‑1 (MQ=255)

                                                                                                   cTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCa                                                    <  1:3032704/150‑1 (MQ=255)

                                                                                                     gcgcTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATg                                                  <  2:7560840/150‑1 (MQ=255)

                                                                                                         tCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGc                                              <  2:1462940/150‑1 (MQ=255)

                                                                                                         tCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGc                                              >  1:8134413/1‑150 (MQ=255)

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

                                                                                                             tCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCt                                          <  2:4024427/150‑1 (MQ=255)

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

                                                                                                                      cATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAgg                                 <  1:1617712/150‑1 (MQ=255)

                                                                                                                         cGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCa                              <  2:3298121/150‑1 (MQ=255)

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

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

                                                                                                                              gcgggcagcaggcagGCCAAGTACAACCTGGCCTGCCTGC‑GCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgcag                        >  2:8561129/8‑150 (MQ=38)

                                                                                                                                 aTGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgcagca                      >  2:8069101/1‑150 (MQ=255)

                                                                                                                                      caggcaggCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCgg                 <  1:7838740/150‑1 (MQ=255)

                                                                                                                                       aggcaggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGc                >  1:7704827/1‑150 (MQ=255)

                                                                                                                                         gcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg              >  1:7515132/1‑150 (MQ=255)

                                                                                                                                           aggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGt            <  1:3718552/150‑1 (MQ=255)

                                                                                                                                             gcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTc          >  2:5607372/1‑150 (MQ=255)

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

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

                                                                                                                                               cAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAg        >  2:5880728/1‑150 (MQ=255)

                                                                                                                                                  gTTGTACTTGGCCTGGCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCTCCCGGCAGGTTCAGTTc     >  2:3741782/1‑150 (MQ=255)

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

                                                                                                                                                     gTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGgc  <  2:5645258/150‑1 (MQ=255)

AACCGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGC  >  NC_002947/5974590‑5974888

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