BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,615,424	G→A	5.0%	A104A (GCG→GCA)	PP_4084 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,615,424	0	G	A	5.0%	443.4 / 7.2	139	A104A (GCG→GCA)	PP_4084	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (70/62); new base A (5/2); total (75/64)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

ACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGACAAGG  >  NC_002947/4615277‑4615568

aCTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAggcggc                                                                                                                                                >  2:4425150/1‑150 (MQ=255)

  tGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGgcct                                                                                                                                              <  2:1609375/150‑1 (MQ=255)

     aCGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGgcctgcc                                                                                                                                           <  2:2274784/150‑1 (MQ=255)

     aCGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGgcctgcc                                                                                                                                           <  2:4711883/150‑1 (MQ=255)

      cGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGgcctgcct                                                                                                                                          <  1:469234/150‑1 (MQ=255)

            aCGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATAcc                                                                                                                                    >  1:1007568/1‑150 (MQ=255)

            aCGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATAcc                                                                                                                                    >  2:4612992/1‑150 (MQ=255)

             cGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCg                                                                                                                                   >  2:2709325/1‑150 (MQ=255)

             cGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCg                                                                                                                                   >  2:4589511/1‑150 (MQ=255)

             cGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCg                                                                                                                                   <  1:1480738/150‑1 (MQ=255)

             cGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCg                                                                                                                                   <  1:185881/150‑1 (MQ=255)

              ggAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCgg                                                                                                                                  <  1:4197782/150‑1 (MQ=255)

                aagaagTCAAAGGTGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGa                                                                                                                                >  2:3624319/1‑150 (MQ=255)

                aagaagTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGa                                                                                                                                <  1:1041750/150‑1 (MQ=255)

                aagaagTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGa                                                                                                                                <  1:4133057/150‑1 (MQ=255)

                  gaagTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACa                                                                                                                              <  1:3876372/150‑1 (MQ=255)

                      tCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCaa                                                                                                                          >  2:3513663/1‑150 (MQ=255)

                        aaaGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACg                                                                                                                        >  2:4808483/1‑150 (MQ=255)

                        aaaGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACg                                                                                                                        >  2:1667460/1‑150 (MQ=255)

                        aaaGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACg                                                                                                                        >  2:1714927/1‑150 (MQ=255)

                         aaGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACgg                                                                                                                       >  1:2759519/1‑150 (MQ=255)

                             cgcTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGc                                                                                                                   <  2:447483/150‑1 (MQ=255)

                               cTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTg                                                                                                                 <  1:984094/150‑1 (MQ=255)

                               cTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTg                                                                                                                 <  1:4629048/150‑1 (MQ=255)

                                tGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTgg                                                                                                                >  1:4576501/1‑150 (MQ=255)

                                tGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTgg                                                                                                                <  1:1169242/150‑1 (MQ=255)

                                  aCCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGACATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTg                                                                                                              >  2:2118252/1‑150 (MQ=255)

                                  aCCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTg                                                                                                              >  2:2118246/1‑150 (MQ=255)

                                  aCCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTg                                                                                                              >  2:4638043/1‑150 (MQ=255)

                                    cATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCt                                                                                                            >  2:1777682/1‑150 (MQ=255)

                                      ttACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGACATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCa                                                                                                          >  1:1407150/1‑150 (MQ=255)

                                       tACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCaa                                                                                                         <  2:4459236/150‑1 (MQ=255)

                                          gCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTa                                                                                                      >  1:1401173/1‑150 (MQ=255)

                                            ccGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGACATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACa                                                                                                    >  2:2865373/1‑150 (MQ=255)

                                              gcgAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACa                                                                                                    <  1:3452332/148‑1 (MQ=255)

                                              gcgAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACa                                                                                                    >  2:3452332/1‑148 (MQ=255)

                                              gcgAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGc                                                                                                  >  1:2885805/1‑150 (MQ=255)

                                                gAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCaa                                                                                                >  1:3718766/1‑150 (MQ=255)

                                                  aaGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGc                                                                                              >  2:3121199/1‑150 (MQ=255)

                                                  aaGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGc                                                                                              <  2:1401173/150‑1 (MQ=255)

                                                       aCCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTggg                                                                                         >  2:345825/1‑150 (MQ=255)

                                                         ccGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGcc                                                                                       >  1:4637698/1‑150 (MQ=255)

                                                         ccGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGcc                                                                                       >  1:5207668/1‑150 (MQ=255)

                                                         ccGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGcc                                                                                       >  1:762374/1‑150 (MQ=255)

                                                           gCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAg                                                                                                   >  1:1114857/1‑136 (MQ=255)

                                                           gCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAg                                                                                                   <  2:1114857/136‑1 (MQ=255)

                                                            ccTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCttt                                                                                    <  2:2909005/150‑1 (MQ=255)

                                                            ccTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCttt                                                                                    <  1:1777682/150‑1 (MQ=255)

                                                              ttGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                                                  <  2:382556/150‑1 (MQ=255)

                                                              ttGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                                                  <  2:2041776/150‑1 (MQ=255)

                                                              ttGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                                                  <  2:1754123/150‑1 (MQ=255)

                                                                      aaTGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTg                                                                                                                 >  1:1424631/1‑111 (MQ=255)

                                                                      aaTGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTg                                                                                                                 <  2:1424631/111‑1 (MQ=255)

                                                                             aTCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCa                                                                   <  2:2885805/150‑1 (MQ=255)

                                                                              tCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAg                                                                  >  2:2813497/1‑150 (MQ=255)

                                                                              tCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAg                                                                  >  2:2813547/1‑150 (MQ=255)

                                                                              tCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAg                                                                  <  1:4172632/150‑1 (MQ=255)

                                                                               cccGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGc                                                                 >  1:2380858/1‑150 (MQ=255)

                                                                               cccGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGc                                                                 <  2:572088/150‑1 (MQ=255)

                                                                                   tCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACt                                                             >  1:4655590/1‑150 (MQ=255)

                                                                                     gACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGACTGGGCCTTTGAGTGACGACGATATCAGCCACTgg                                                           >  1:1176782/1‑150 (MQ=255)

                                                                                          tCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGgcgcgc                                                      <  2:4176254/150‑1 (MQ=255)

                                                                                          tCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGgcgcgc                                                      >  1:3202664/1‑150 (MQ=255)

                                                                                            gggAAATGGCCATCCAGTTCATGGTCGGCGACAAGTCCCACTACCAGCTTCAGGCAGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGgcgcgcgc                                                    >  2:2164886/1‑150 (MQ=255)

                                                                                             ggAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCt                                                   >  2:2874919/1‑150 (MQ=255)

                                                                                             ggAAATGGCCATCCAGTTCATGGTCGGCGACAAACCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCt                                                   >  2:2117674/1‑150 (MQ=255)

                                                                                                 aTGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTaaa                                               <  1:650606/150‑1 (MQ=255)

                                                                                                    gCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGAc                                            <  1:1597470/150‑1 (MQ=255)

                                                                                                      cATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGa                                          <  1:1437278/150‑1 (MQ=255)

                                                                                                      cATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGa                                          <  1:3198070/150‑1 (MQ=255)

                                                                                                      cATCCAGTTCATGGTCGGCGACAAGCCACACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGa                                          <  2:1191147/150‑1 (MQ=255)

                                                                                                       aTCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAg                                         >  1:2561435/1‑150 (MQ=255)

                                                                                                       aTCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAg                                         >  1:3508696/1‑150 (MQ=255)

                                                                                                         ccAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCt                                       >  2:2378043/1‑150 (MQ=255)

                                                                                                           aGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTcc                                     <  2:2492423/150‑1 (MQ=255)

                                                                                                            gTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGa                                                                                                                                        <  1:1852322/50‑1 (MQ=255)

                                                                                                            gTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGa                                                                                                                                        >  2:1852322/1‑50 (MQ=255)

                                                                                                             ttCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGa                                   <  2:5030870/150‑1 (MQ=255)

                                                                                                             ttCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGa                                   <  1:1683884/150‑1 (MQ=255)

                                                                                                               cATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGAtt                                 <  2:3733328/150‑1 (MQ=255)

                                                                                                                aTGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTGTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGAttt                                >  1:2795987/1‑150 (MQ=255)

                                                                                                                  ggTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGc                              <  2:1194860/150‑1 (MQ=255)

                                                                                                                  ggTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGc                              <  1:3015099/150‑1 (MQ=255)

                                                                                                                  ggTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGc                              >  1:3811820/1‑150 (MQ=255)

                                                                                                                        cGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTg                                                                               >  1:4678913/1‑95 (MQ=255)

                                                                                                                        cGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTg                                                                               <  2:4678913/95‑1 (MQ=255)

                                                                                                                             aGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCtctc                   <  1:2650358/150‑1 (MQ=255)

                                                                                                                                 ccACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTgg                                                                                          >  2:3806754/1‑75 (MQ=255)

                                                                                                                                      aCCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCa          <  2:4183004/150‑1 (MQ=255)

                                                                                                                                       ccggtatCAGGCAGGCCGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAg         >  2:213104/7‑150 (MQ=255)

                                                                                                                                        cggtatCAGGCAGGCCGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTaaa                                               >  2:832600/6‑111 (MQ=255)

                                                                                                                                        cggtatCAGGCAGGCCGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTaaa                                               <  1:832600/106‑1 (MQ=255)

                                                                                                                                        cggtatCAGGCAGGCCGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAgg        >  2:244179/6‑150 (MQ=255)

                                                                                                                                        cAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAgg        >  2:1210479/1‑150 (MQ=255)

                                                                                                                                         aGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGa       <  1:4491970/150‑1 (MQ=255)

                                                                                                                                           cTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGACa     >  1:2880622/1‑150 (MQ=255)

                                                                                                                                            ttCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGTTATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGACaa    <  1:3624319/150‑1 (MQ=255)

                                                                                                                                            ttCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGACaa    >  1:4850666/1‑150 (MQ=255)

                                                                                                                                            atCAGGCAGGCCGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTg                                                            >  1:3695938/2‑94 (MQ=255)

                                                                                                                                              cAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGACAAgg  <  2:3013683/150‑1 (MQ=255)

ACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGCCTCTCCACGCCCCAGGACAAGG  >  NC_002947/4615277‑4615568

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