BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,109,690	C→T	5.8%	E62E (GAG→GAA)	PP_1879 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,109,690	0	C	T	5.8%	980.3 / 32.4	307	E62E (GAG→GAA)	PP_1879	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (153/136); new base T (13/5); total (166/141)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 311 total aligned reads displayed.

TCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTT  >  NC_002947/2109541‑2109839

tCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAAc                                                                                                                                                       <  2:5740825/150‑1 (MQ=255)

    tAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCa                                                                                                                                                   <  2:4229417/150‑1 (MQ=255)

      aGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAcc                                                                                                                                                 <  2:6973524/150‑1 (MQ=255)

       gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAccc                                                                                                                                                >  1:3622892/1‑150 (MQ=255)

        cTTGTCGACCGGTACCAATATGGATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCg                                                                                                                                               >  1:5042507/1‑150 (MQ=255)

           gTCGACCGGTACCAAGATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAg                                                                                                                                            <  1:5808019/150‑1 (MQ=255)

               aCCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                        >  1:4391670/1‑150 (MQ=255)

                 cGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGa                                                                                                                                      <  2:6310790/150‑1 (MQ=255)

                   gTACCAATATGAATACCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAAt                                                                                                                                    <  1:3251577/150‑1 (MQ=255)

                    tACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATa                                                                                                                                   <  2:3909108/150‑1 (MQ=255)

                     aCCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAg                                                                                                                                  <  2:1085966/150‑1 (MQ=255)

                      ccAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGACTAgg                                                                                                                                 >  1:3565395/1‑150 (MQ=255)

                      ccAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAgg                                                                                                                                 >  1:3122628/1‑150 (MQ=255)

                         atatGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTg                                                                                                                              >  1:1489606/1‑150 (MQ=255)

                                tAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGaa                                                                                                                       <  2:4844005/150‑1 (MQ=255)

                                   ctctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTc                                                                                                                    <  1:2783313/150‑1 (MQ=255)

                                      ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGAc                                                                                                                 >  1:1777005/1‑150 (MQ=255)

                                           gTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCgtgt                                                                                                            >  1:1366120/1‑150 (MQ=255)

                                               gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTggg                                                                                                        >  2:3466817/1‑150 (MQ=255)

                                               gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTggg                                                                                                        >  2:5622349/1‑150 (MQ=255)

                                               gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTggg                                                                                                        >  2:2575009/1‑150 (MQ=255)

                                                       tGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTccccc                                                                                                <  2:2487673/150‑1 (MQ=255)

                                                        gCTGAAGGTCTTCTGGTCGCCCAGTACGTCAGTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTGTcccccc                                                                                               <  2:1584683/150‑1 (MQ=255)

                                                                tcttctTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCgt                                                                                                                                              >  1:6280618/1‑94 (MQ=255)

                                                                tcttctTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCgt                                                                                                                                              <  2:6280618/95‑2 (MQ=255)

                                                                 cttcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATc                                                                                      >  2:6119613/1‑150 (MQ=255)

                                                                  ttcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                                     >  2:2778746/1‑150 (MQ=255)

                                                                    cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTc                                                                                   <  2:3279586/150‑1 (MQ=255)

                                                                       gTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGGGTTGGGTTTCCCCCCTTCAACATCGTccac                                                                                >  1:4497556/1‑150 (MQ=255)

                                                                        tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCa                                                                                                                                             <  1:4961369/88‑1 (MQ=255)

                                                                           aCCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACg                                                                            >  1:2904981/1‑150 (MQ=255)

                                                                             cAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCt                                                                          >  1:4449865/1‑150 (MQ=255)

                                                                             cAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCt                                                                          >  2:3186320/1‑150 (MQ=255)

                                                                                tACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTgg                                                                       >  1:1642792/1‑150 (MQ=255)

                                                                                tACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTgg                                                                       >  2:6152166/1‑150 (MQ=255)

                                                                                tACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTgg                                                                       >  2:2900778/1‑150 (MQ=255)

                                                                                     cATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACgg                                                                  <  1:6930725/150‑1 (MQ=255)

                                                                                     cATTCCAGTAGAGG‑AAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGt                                                                 <  2:4120870/150‑1 (MQ=255)

                                                                                        tCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTg                                                               >  1:190450/1‑150 (MQ=255)

                                                                                        tCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTg                                                               >  1:2527932/1‑150 (MQ=255)

                                                                                         ccAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgg                                                              <  1:744941/150‑1 (MQ=255)

                                                                                         ccAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgg                                                              <  1:6036313/150‑1 (MQ=255)

                                                                                          cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCg                                                                                                                   >  1:5127801/1‑96 (MQ=255)

                                                                                          cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                             <  2:7400432/150‑1 (MQ=255)

                                                                                          cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                             <  2:7276395/150‑1 (MQ=255)

                                                                                          cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                             <  1:796195/150‑1 (MQ=255)

                                                                                          cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                             <  2:3588897/150‑1 (MQ=255)

                                                                                           aGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgggg                                                            >  2:4026008/1‑150 (MQ=255)

                                                                                           aGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgggg                                                            >  1:500464/1‑150 (MQ=255)

                                                                                            gTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                           >  1:6256304/1‑150 (MQ=255)

                                                                                            gTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAAGAGGTAGCCAGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                           >  1:1253978/1‑150 (MQ=255)

                                                                                             tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAc                                                                                                                                                  <  1:1066578/62‑1 (MQ=255)

                                                                                               gagGAAAAGGTATTCCTCATCTATTTCCGGGTTCTAATTTGCCCACTGGTTGAACTACACCCGCAGGGGGGAATAGGATGCCCGGAACTAGACCGTGATGGGTTTCCCCCCTTCAACATCGTCCACCACTCTTGGTACGGAAGGGGGGAc                                                        >  1:5518301/1‑150 (MQ=255)

                                                                                               gagGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGAc                                                        >  1:6878658/1‑150 (MQ=255)

                                                                                                agGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACa                                                       >  2:2329291/1‑150 (MQ=255)

                                                                                                  gAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCgt                                                                                                              >  2:5979179/1‑93 (MQ=255)

                                                                                                    aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagag                                                   >  1:3024901/1‑150 (MQ=255)

                                                                                                     aaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagaga                                                  <  2:5121035/150‑1 (MQ=255)

                                                                                                      aGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagagag                                                 <  2:1172931/150‑1 (MQ=255)

                                                                                                      aGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagagag                                                 <  1:2170252/150‑1 (MQ=255)

                                                                                                      aGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagagag                                                 <  2:4973562/150‑1 (MQ=255)

                                                                                                         tATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGTcc                                              <  2:1301941/150‑1 (MQ=255)

                                                                                                           ttCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGa                                            <  1:1174461/150‑1 (MQ=255)

                                                                                                           ttCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGa                                            <  2:6831683/150‑1 (MQ=255)

                                                                                                           ttCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGa                                            <  1:4557436/150‑1 (MQ=255)

                                                                                                            tCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAg                                           <  2:4391670/150‑1 (MQ=255)

                                                                                                              cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGGTGGGGGGACAAGAGAGGCCGAGTc                                         >  2:6221931/1‑150 (MQ=255)

                                                                                                               tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTcc                                        >  1:2206419/1‑150 (MQ=255)

                                                                                                               tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTcc                                        >  2:6615256/1‑150 (MQ=255)

                                                                                                               tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTGGGGGGGACAAGAGAGGCCGAGTcc                                        >  2:7160063/1‑150 (MQ=255)

                                                                                                               tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTAGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGTGGGGACAAGAGAGGCCGAGTcc                                        >  2:674414/1‑150 (MQ=255)

                                                                                                               tCATCTATTACCGGGTGCGAATTTGACCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAAAATCGTCCAACACGCTTGGTACGGATGGGGGGACAAGAGAGGCCGAGTcc                                        >  2:7210315/1‑150 (MQ=255)

                                                                                                                cATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCa                                       >  2:1468989/1‑150 (MQ=255)

                                                                                                                cATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCa                                       <  1:1517500/150‑1 (MQ=255)

                                                                                                                    tATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGa                                                         >  2:5286835/1‑128 (MQ=255)

                                                                                                                    tATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGa                                                         <  1:5286835/128‑1 (MQ=255)

                                                                                                                      tttCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                                 >  1:4311564/1‑150 (MQ=255)

                                                                                                                           gggTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGtt                            <  2:5445683/150‑1 (MQ=255)

                                                                                                                           gggTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGtt                            <  2:4923395/150‑1 (MQ=255)

                                                                                                                            ggTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                           <  1:6738473/150‑1 (MQ=255)

                                                                                                                              tGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGa                         <  2:1816388/150‑1 (MQ=255)

                                                                                                                              tGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGa                         <  1:2976651/150‑1 (MQ=255)

                                                                                                                              tGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGa                         >  1:6643557/1‑150 (MQ=255)

                                                                                                                                 gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGtt                      <  1:3994184/150‑1 (MQ=255)

                                                                                                                                   aTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCt                    <  2:5042507/150‑1 (MQ=255)

                                                                                                                                   aTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCt                    <  1:5397526/150‑1 (MQ=255)

                                                                                                                                    tttGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTg                   >  1:1994421/1‑150 (MQ=255)

                                                                                                                                     ttGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGc                  <  2:5197024/150‑1 (MQ=255)

                                                                                                                                      tGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg                 <  1:5668336/150‑1 (MQ=255)

                                                                                                                                      cccctgcGGGTGGAGTTCCACCCGTGGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg                 >  2:4725610/7‑150 (MQ=255)

                                                                                                                                      cccctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTCCg                 >  2:828013/7‑150 (MQ=255)

                                                                                                                                         ctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGtt              >  2:2016759/4‑150 (MQ=255)

                                                                                                                                          cACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTa             <  1:5242057/150‑1 (MQ=255)

                                                                                                                                             gggTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCg          <  1:3492569/150‑1 (MQ=255)

                                                                                                                                               gTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGa        <  1:7236020/150‑1 (MQ=255)

                                                                                                                                                tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCGAGTTCAGGGGTTGGAGTTCTGCGGTTAGCGGAg       >  2:2446032/1‑150 (MQ=255)

                                                                                                                                                tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAg       >  1:4771175/1‑150 (MQ=255)

                                                                                                                                                tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAg       <  2:3025009/150‑1 (MQ=255)

                                                                                                                                                    aCTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCt   >  2:1031767/1‑150 (MQ=255)

                                                                                                                                                     cTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCtt  <  2:864945/150‑1 (MQ=255)

TCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTT  >  NC_002947/2109541‑2109839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`