BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,205,208	T→A	5.5%	N275Y (AAT→TAT)	PP_3690 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,205,208	0	T	A	5.5%	666.9 / 17.1	196	N275Y (AAT→TAT)	PP_3690	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (92/92); new base A (4/7); total (96/100)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

AAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCACTTCTC  >  NC_002947/4205059‑4205357

aaGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTAtt                                                                                                                                                       <  2:3720751/150‑1 (MQ=255)

  gCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGa                                                                                                                                                     <  2:3470179/150‑1 (MQ=255)

   cTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAg                                                                                                                                                    >  1:1691735/1‑150 (MQ=255)

   cTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAg                                                                                                                                                    <  2:4621389/150‑1 (MQ=255)

    tGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGc                                                                                                                                                   >  2:5310830/1‑150 (MQ=255)

       aTCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTg                                                                                                                                                <  2:5483390/150‑1 (MQ=255)

        tCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGt                                                                                                                                               <  2:3220512/150‑1 (MQ=255)

         cACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGtt                                                                                                                                              >  2:1157071/1‑150 (MQ=255)

         cACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggc                                                                                                                                              <  1:3210147/150‑3 (MQ=255)

         cACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggc                                                                                                                                              <  1:2054410/150‑3 (MQ=255)

          aCGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTc                                                                                                                                             >  1:3568510/1‑150 (MQ=255)

          aCGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggca                                                                                                                                             <  2:2192184/150‑4 (MQ=255)

          aCGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggca                                                                                                                                             <  2:701062/150‑4 (MQ=255)

            ggCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCa                                                                                                                                           <  2:3180549/150‑1 (MQ=255)

             gCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAg                                                                                                                                          >  2:5279560/1‑150 (MQ=255)

             gCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAg                                                                                                                                          <  2:1726099/150‑1 (MQ=255)

               tGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCt                                                                                                                                        <  1:1901725/150‑1 (MQ=255)

               tGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCt                                                                                                                                        <  1:1401051/150‑1 (MQ=255)

                ggCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTg                                                                                                                                       >  2:1785107/1‑150 (MQ=255)

                     cGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATa                                                                                                                                  >  1:4276330/1‑150 (MQ=255)

                      gcttgcttGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATaa                                                                                                                                 >  1:777641/1‑150 (MQ=255)

                       cttgcttgGAAATTGCGGCAAAATCGTTGTCTGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATaaa                                                                                                                                >  2:972969/1‑150 (MQ=255)

                        ttgcttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATaaaa                                                                                                                               <  2:4086639/150‑1 (MQ=255)

                        ttgcttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATaaaa                                                                                                                               <  1:2422465/150‑1 (MQ=255)

                            ttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACg                                                                                                                           >  1:5408496/1‑150 (MQ=255)

                            ttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACg                                                                                                                           >  2:4995438/1‑150 (MQ=255)

                            ttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACg                                                                                                                           >  1:4044516/1‑150 (MQ=255)

                               gAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCt                                                                                                                        >  2:1687655/1‑150 (MQ=255)

                                                  tGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATc                                                                                                     <  2:2416831/150‑1 (MQ=255)

                                                    tCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTg                                                                                                   >  1:3154314/1‑150 (MQ=255)

                                                     ccGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGa                                                                                                  <  2:691899/150‑1 (MQ=255)

                                                     ccGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGa                                                                                                  <  2:759532/150‑1 (MQ=255)

                                                      cGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGAc                                                                                                 <  2:3388580/150‑1 (MQ=255)

                                                           cGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTa                                                                                            <  1:5273913/150‑1 (MQ=255)

                                                            gagaTCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAg                                                                                           >  1:5479162/1‑150 (MQ=255)

                                                                  ttCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTaaa                                                                                     >  2:137120/1‑150 (MQ=255)

                                                                  ttCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTaaa                                                                                     >  1:5358712/1‑150 (MQ=255)

                                                                   tCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTaaaa                                                                                    >  1:5147943/1‑150 (MQ=255)

                                                                      aaGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACa                                                                                 <  2:4943501/150‑1 (MQ=255)

                                                                      aaGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACa                                                                                 >  1:2109126/1‑150 (MQ=255)

                                                                        gTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGa                                                                               <  1:4182463/150‑1 (MQ=255)

                                                                        gTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGa                                                                               <  1:1687655/150‑1 (MQ=255)

                                                                        gTACGCCTTCAGAACGCCCTCGATACATACGATAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAgt                                                                               <  2:964711/150‑2 (MQ=255)

                                                                          aCGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGAcc                                                                             <  1:3952480/150‑1 (MQ=255)

                                                                          aCGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGAcc                                                                             <  1:4482989/150‑1 (MQ=255)

                                                                            gCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCg                                                                           >  1:1784268/1‑150 (MQ=255)

                                                                               ttCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaag                                                                        >  1:565874/1‑150 (MQ=255)

                                                                               ttCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaag                                                                        <  1:5495091/150‑1 (MQ=255)

                                                                               ttCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaag                                                                        >  1:3020493/1‑150 (MQ=255)

                                                                                 cAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCgaagaa                                                                      <  2:3469389/150‑1 (MQ=255)

                                                                                   gAACTCCCTCGTTACATACGCTATCGTTCAGTGCTGCCGTATCTATGCCTGCCTTCCAAGCTCTATGGAGCTTGTTCCAGCTGTATTTAAACTCTCCTCCTTCACGCTCAAAAAATCTGACGATTAGGATAAAAACAGACCCGATGAAGt                                                                    <  2:4574995/150‑1 (MQ=255)

                                                                                   gAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGt                                                                    <  2:2029958/150‑1 (MQ=255)

                                                                                    aaCGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTc                                                                   >  2:5510687/1‑150 (MQ=255)

                                                                                    aaCGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTc                                                                   <  1:3404549/150‑1 (MQ=255)

                                                                                        cccTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGaa                                                               <  2:3044446/150‑1 (MQ=255)

                                                                                         ccTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAg                                                              >  1:1868664/1‑150 (MQ=255)

                                                                                         ccTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAg                                                              >  2:631604/1‑150 (MQ=255)

                                                                                           tCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAATATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGca                                                            >  1:3254701/1‑150 (MQ=255)

                                                                                           tCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGca                                                            <  1:4938127/150‑1 (MQ=255)

                                                                                             gATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGcaca                                                          >  1:2924022/1‑150 (MQ=255)

                                                                                                acatacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAAtt                                                       <  1:1660421/150‑1 (MQ=255)

                                                                                                acatacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAAtt                                                       <  1:5533539/150‑1 (MQ=255)

                                                                                                  atacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGa                                                     >  1:817949/1‑150 (MQ=255)

                                                                                                  atacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGa                                                     >  1:4874132/1‑150 (MQ=255)

                                                                                                   tacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGaa                                                    <  2:4076171/150‑1 (MQ=255)

                                                                                                   tacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGaa                                                    <  2:2006638/150‑1 (MQ=255)

                                                                                                   tacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGaa                                                    <  1:4858625/150‑1 (MQ=255)

                                                                                                   tacGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGaa                                                    <  2:2410447/150‑1 (MQ=255)

                                                                                                     cGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAgg                                                  <  1:5021094/150‑1 (MQ=255)

                                                                                                      gCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAggg                                                 >  1:3762817/1‑150 (MQ=255)

                                                                                                       cTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGt                                                >  1:4251044/1‑150 (MQ=255)

                                                                                                        tAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTa                                               >  1:719912/1‑150 (MQ=255)

                                                                                                        tAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTa                                               >  2:1354520/1‑150 (MQ=255)

                                                                                                         aaCGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTAt                                              >  2:617299/1‑150 (MQ=255)

                                                                                                            gTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTAtctc                                           <  2:4172508/150‑1 (MQ=255)

                                                                                                              tCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGa                                         >  2:2058692/1‑150 (MQ=255)

                                                                                                                aGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACa                                       <  1:964060/150‑1 (MQ=255)

                                                                                                                  tgctgcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAgg                                     >  1:1083198/1‑150 (MQ=255)

                                                                                                                      gcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGa                                 >  2:4284476/1‑150 (MQ=255)

                                                                                                                      gcCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGa                                 >  1:2532651/1‑150 (MQ=255)

                                                                                                                       ccGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGaa                                <  2:127049/150‑1 (MQ=255)

                                                                                                                           aTCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAg                            <  2:4241979/150‑1 (MQ=255)

                                                                                                                            tCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGt                           <  2:3568510/150‑1 (MQ=255)

                                                                                                                            tCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGt                           <  2:1691735/150‑1 (MQ=255)

                                                                                                                            tCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGt                           >  2:2693675/1‑150 (MQ=255)

                                                                                                                                  cctgcctgccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGa                     >  1:3123622/1‑150 (MQ=255)

                                                                                                                                   ctgcctgccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGAt                    >  2:3992388/1‑150 (MQ=255)

                                                                                                                                   ctgcctgccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGAt                    >  1:787357/1‑150 (MQ=255)

                                                                                                                                    tgcctgccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGtt                          <  2:1912656/143‑1 (MQ=255)

                                                                                                                                    tgcctgccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGAtt                   >  1:2199687/1‑150 (MQ=255)

                                                                                                                                      tggaacAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTAc                 >  1:3465448/6‑150 (MQ=255)

                                                                                                                                      tggaacAAGCTCAATAGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTAc                 >  2:2819762/6‑150 (MQ=255)

                                                                                                                                      tggaacAAGCTCAATAGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTAc                 >  2:1576674/6‑150 (MQ=255)

                                                                                                                                        gaacAAGCTCAATAGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCAGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGAAAGGCGGACGGAGTTGACGATAACTg               >  1:2896242/4‑150 (MQ=255)

                                                                                                                                          ccAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGACGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGaa             <  2:4044516/150‑1 (MQ=255)

                                                                                                                                               cactAATGAGCTTGTACCAGCTGACATAAAACACCCCTCCCACACGCACAAAAAATCTGACGAGTAGGCTAAAAACACACCCGAAGAAGTCAGAAGCACAATTGAACGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCa        <  2:4276330/148‑1 (MQ=255)

                                                                                                                                                  tATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCACtt     >  1:228021/1‑150 (MQ=255)

                                                                                                                                                   aTTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCACTtc    >  2:4973193/1‑150 (MQ=255)

                                                                                                                                                    ttGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCACTtct   >  1:1479829/1‑150 (MQ=255)

                                                                                                                                                     tGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCACTtctc  <  2:5479130/150‑1 (MQ=255)

AAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAACACGCCTCCGTCACGCTCAAAAAATCTGACGAGTAGGCTAAAAACAGACCCGAAGAAGTCAGAAGCACAATTGAAGGGTATCTCGACAGGCGGAAGGAGTTGACGATTACTGAACTCCACTTCTC  >  NC_002947/4205059‑4205357

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