BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,116,982	T→A	6.0%	H118Q (CAT→CAA)	PP_5362 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,116,982	0	T	A	6.0%	720.7 / 25.0	218	H118Q (CAT→CAA)	PP_5362	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (98/107); new base A (5/8); total (103/115)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 224 total aligned reads displayed.

TGGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAACCG  >  NC_002947/6116833‑6117131

tgGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAt                                                                                                                                                       >  2:2846368/1‑150 (MQ=255)

tgGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAt                                                                                                                                                       >  1:2259164/1‑150 (MQ=255)

    cacaAACCAACCGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGTCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                                   <  1:3767036/150‑1 (MQ=255)

    cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                                   >  1:5017063/1‑150 (MQ=255)

    cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                                   <  1:2644325/150‑1 (MQ=255)

       aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCt                                                                                                                                                >  1:4177598/1‑150 (MQ=255)

        aaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCtt                                                                                                                                               >  1:1045303/1‑150 (MQ=255)

           cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGAt                                                                                                                                            <  2:2469945/150‑1 (MQ=255)

           cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                            <  2:2013119/150‑2 (MQ=255)

           cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                            <  1:3632372/150‑2 (MQ=255)

            aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATc                                                                                                                                           <  2:2802992/150‑1 (MQ=255)

             aCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCg                                                                                                                                          >  1:5537317/1‑150 (MQ=255)

             aCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaagc                                                                                                                                          <  1:3089503/150‑4 (MQ=255)

               aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                        >  1:821010/1‑150 (MQ=255)

               aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                        <  1:4130942/150‑1 (MQ=255)

               aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                        >  1:3477875/1‑150 (MQ=255)

                gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                       >  2:3346818/1‑150 (MQ=255)

                     tATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGc                                                                                                                                  <  1:1755270/150‑1 (MQ=255)

                         gCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGtt                                                                                                                              >  2:1822025/1‑150 (MQ=255)

                         gCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGtt                                                                                                                              <  1:3565432/150‑1 (MQ=255)

                           acaTTGACTTTTGTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGa                                                                                                                            <  1:26878/150‑1 (MQ=255)

                            caTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGAc                                                                                                                           <  2:5017063/150‑1 (MQ=255)

                            caTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGAc                                                                                                                           <  1:2609910/150‑1 (MQ=255)

                            caTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCGGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGAc                                                                                                                           <  2:1104406/150‑1 (MQ=255)

                              ttGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCATATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGc                                                                                                                         <  2:3337360/150‑1 (MQ=255)

                               tGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcc                                                                                                                        >  2:2999216/1‑150 (MQ=255)

                                 aCTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                                      >  2:4427946/1‑150 (MQ=255)

                                   tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                    >  2:1646292/1‑150 (MQ=255)

                                   tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                    >  2:4481943/1‑150 (MQ=255)

                                   tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                    >  1:5151577/1‑150 (MQ=255)

                                   tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                    >  2:743010/1‑150 (MQ=255)

                                   tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                    <  2:10125/150‑1 (MQ=255)

                                    ttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCtt                                                                                                                   <  2:4297922/150‑1 (MQ=255)

                                    ttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCtt                                                                                                                   >  2:1026581/1‑150 (MQ=255)

                                         tcccGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtt                                                                                                              <  2:5616521/148‑1 (MQ=255)

                                         tACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtt                                                                                                              <  2:825527/150‑1 (MQ=255)

                                         tACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtt                                                                                                              <  1:3687430/150‑1 (MQ=255)

                                            cGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTg                                                                                                           >  2:2277932/1‑150 (MQ=255)

                                            cGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTg                                                                                                           >  2:2233547/1‑150 (MQ=255)

                                              gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGaa                                                                                                         >  1:3354554/1‑150 (MQ=255)

                                                   gCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttt                                                                                                    <  2:356435/150‑1 (MQ=255)

                                                    ccAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttt                                                                                                   <  1:2095649/150‑1 (MQ=255)

                                                        cGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtct                                                                                               >  1:1857749/1‑150 (MQ=255)

                                                         gATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtctc                                                                                              >  1:1962811/1‑150 (MQ=255)

                                                             gTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgct                                                                                          >  2:3261675/1‑150 (MQ=255)

                                                             gTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgct                                                                                          >  1:3423401/1‑150 (MQ=255)

                                                             gTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgct                                                                                          >  2:1263331/1‑150 (MQ=255)

                                                                  gCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTa                                                                                     >  1:431079/1‑150 (MQ=255)

                                                                  gCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTa                                                                                     >  2:4964530/1‑150 (MQ=255)

                                                                    cTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtt                                                                                   <  2:2253359/150‑1 (MQ=255)

                                                                     tgggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAttt                                                                                  <  2:1651693/150‑1 (MQ=255)

                                                                      gggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtttt                                                                                 <  2:4072379/150‑1 (MQ=255)

                                                                      gggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtttt                                                                                 <  1:3346818/150‑1 (MQ=255)

                                                                       ggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTg                                                                                <  1:2390408/150‑1 (MQ=255)

                                                                        gtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                               <  1:2864009/150‑1 (MQ=255)

                                                                        gtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                               <  1:521917/150‑1 (MQ=255)

                                                                          gggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTaa                                                                             <  1:1443860/150‑1 (MQ=255)

                                                                             cTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTc                                                                          >  1:65788/1‑150 (MQ=255)

                                                                             cTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTc                                                                          >  1:65754/1‑150 (MQ=255)

                                                                               aCAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAg                                                                        <  2:525788/150‑1 (MQ=255)

                                                                                   ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTc                                                                    >  1:3084169/1‑150 (MQ=255)

                                                                                    cgccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCg                                                                   <  1:3906821/150‑1 (MQ=255)

                                                                                     gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTTAGTCAGCTTCGa                                                                  <  2:3828182/150‑1 (MQ=255)

                                                                                     gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                                  <  1:1518931/150‑1 (MQ=255)

                                                                                     gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                                  >  1:3187122/1‑150 (MQ=255)

                                                                                     gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                                  >  1:1950050/1‑150 (MQ=255)

                                                                                      ccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAg                                                                 <  2:3932198/150‑1 (MQ=255)

                                                                                          gTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCAcg                                                             <  1:5420572/150‑1 (MQ=255)

                                                                                              gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTGTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCt                                                         >  1:462823/1‑150 (MQ=255)

                                                                                                tGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGa                                                       <  1:931750/150‑1 (MQ=255)

                                                                                                 ggACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaa                                                      <  2:4854001/150‑1 (MQ=255)

                                                                                                    ccaccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAac                                                   <  2:2072295/150‑1 (MQ=255)

                                                                                                     caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAaca                                                  >  1:447811/1‑150 (MQ=255)

                                                                                                      accTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGGTCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAacac                                                 >  2:2225396/1‑150 (MQ=255)

                                                                                                      accTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAacac                                                 >  2:3195560/1‑150 (MQ=255)

                                                                                                        cTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                               <  2:2876018/150‑1 (MQ=255)

                                                                                                        cTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                               >  1:3303352/1‑150 (MQ=255)

                                                                                                        cTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                               >  2:3999001/1‑150 (MQ=255)

                                                                                                            aaCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                                  <  1:1346408/127‑1 (MQ=255)

                                                                                                            aaCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCt                                           >  1:5386373/1‑150 (MQ=255)

                                                                                                               cgcgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGa                                        <  1:4044962/150‑1 (MQ=255)

                                                                                                               cgcgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGa                                        <  2:3277863/150‑1 (MQ=255)

                                                                                                               cgcgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGa                                        >  2:22866/1‑150 (MQ=255)

                                                                                                                    aTCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAc                                   <  1:2996376/150‑1 (MQ=255)

                                                                                                                      ccAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAccc                                 <  1:1023966/150‑1 (MQ=255)

                                                                                                                      ccAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAccc                                 <  1:2643365/150‑1 (MQ=255)

                                                                                                                         gTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCttt                              <  2:236169/150‑1 (MQ=255)

                                                                                                                          ttCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTg                             <  1:4410238/150‑1 (MQ=255)

                                                                                                                          ttCCAGCGGGGAGCGATTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTg                             >  2:250612/1‑150 (MQ=255)

                                                                                                                             cAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCa                          <  1:4253258/150‑1 (MQ=255)

                                                                                                                               gCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGa                        <  2:364392/150‑1 (MQ=255)

                                                                                                                                cGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGAc                       <  2:1962698/150‑1 (MQ=255)

                                                                                                                                 ggggAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                               <  2:3519/93‑1 (MQ=255)

                                                                                                                                      gtgcTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa                 <  1:2441377/148‑1 (MQ=255)

                                                                                                                                          ttCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTa             >  1:4469837/1‑150 (MQ=255)

                                                                                                                                           tCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTAt            <  2:4501208/150‑1 (MQ=255)

                                                                                                                                           tCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTAt            <  2:676311/150‑1 (MQ=255)

                                                                                                                                            cAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATg           >  1:3572947/1‑150 (MQ=255)

                                                                                                                                            cAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATg           >  2:1572938/1‑150 (MQ=255)

                                                                                                                                                  gCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTaa     >  2:5544894/1‑150 (MQ=255)

TGGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAACCG  >  NC_002947/6116833‑6117131

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