BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	966,598	A→T	7.3%	L41* (TTG→TAG) ‡	PP_0828 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	966,598	0	A	T	7.3%	514.0 / 21.7	159	L41* (TTG→TAG) ‡	PP_0828	membrane protein
Reads supporting (aligned to +/- strand): ref base A (78/69); new base T (7/5); total (85/74)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

AGCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATG  >  NC_002947/966449‑966747

aGCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCa                                                                                                                                                       <  1:883818/150‑1 (MQ=255)

        gAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGc                                                                                                                                               >  2:1120223/1‑150 (MQ=255)

          aCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             <  2:2723759/150‑2 (MQ=255)

          aCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             <  1:3626415/150‑2 (MQ=255)

          aCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             <  2:3804458/150‑2 (MQ=255)

              ttCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGaaa                                                                                                                                         <  2:3953444/150‑1 (MQ=255)

                gaggctcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATc                                                                                                                                       >  2:268444/7‑150 (MQ=255)

                cACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATc                                                                                                                                       >  1:1480829/1‑150 (MQ=255)

                 aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      >  2:3933887/1‑150 (MQ=255)

                 aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      >  1:2219997/1‑150 (MQ=255)

                 aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      >  1:4112168/1‑150 (MQ=255)

                 aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      <  1:4587375/150‑1 (MQ=255)

                   cAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    <  2:4599606/150‑1 (MQ=255)

                   cAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    >  2:3707639/1‑150 (MQ=255)

                    ctcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATg                                                                                                                                   <  1:268444/148‑1 (MQ=255)

                      cgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGt                                                                                                                                 <  2:4895583/150‑1 (MQ=255)

                         gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                              >  1:1227780/1‑150 (MQ=255)

                         gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                              <  1:4294817/150‑1 (MQ=255)

                           cccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCa                                                                                                                            >  1:3635703/1‑150 (MQ=255)

                           cccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCa                                                                                                                            >  1:3473499/1‑150 (MQ=255)

                             cacaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGc                                                                                                                          >  1:2832338/1‑150 (MQ=255)

                              acaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGca                                                                                                                         <  1:3306555/150‑1 (MQ=255)

                                     gCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCa                                                                                                                  <  1:548757/150‑1 (MQ=255)

                                       gAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGa                                                                                                                <  1:4539158/150‑1 (MQ=255)

                                          ccTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                             <  2:3137299/150‑1 (MQ=255)

                                           cTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcag                                                                                                            >  2:2350625/1‑150 (MQ=255)

                                               caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  2:2960507/1‑150 (MQ=255)

                                               caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  1:1533513/1‑150 (MQ=255)

                                               caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  1:1533711/1‑150 (MQ=255)

                                                 ggtatagCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCa                                                                                                      <  1:4024424/144‑1 (MQ=255)

                                                  caGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAg                                                                                                     >  2:599484/1‑150 (MQ=255)

                                                      ggCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                                 <  1:2313129/150‑1 (MQ=255)

                                                       gCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAg                                                                                                <  2:2045033/150‑1 (MQ=255)

                                                        cTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               <  1:3049841/150‑1 (MQ=255)

                                                        cTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               <  1:1338062/150‑1 (MQ=255)

                                                           tAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAGGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAc                                                                                            <  1:1550290/150‑1 (MQ=255)

                                                            aGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAcc                                                                                           <  1:4427461/150‑1 (MQ=255)

                                                            aGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAcc                                                                                           <  1:4427478/150‑1 (MQ=255)

                                                               cgcTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATTAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGc                                                                                        <  2:3658340/150‑1 (MQ=255)

                                                               cgcTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGc                                                                                        <  2:515731/150‑1 (MQ=255)

                                                                  cgctgTAGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGGGACCAGCACCTCCCAGCTCAGCAGCCAGAACAGCGACCAGCCGa                                                                                     <  2:1086893/149‑1 (MQ=255)

                                                                     tgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGatat                                                                                  >  2:4144843/1‑150 (MQ=255)

                                                                           gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            <  2:396093/150‑1 (MQ=255)

                                                                           gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            <  2:396283/150‑1 (MQ=255)

                                                                           gAAAACTCACCAGCACCACCAGCGCAGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            >  1:5059965/1‑150 (MQ=255)

                                                                            aaaaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGAt                                                                           <  1:751982/150‑1 (MQ=255)

                                                                             aaaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCGGATATAGCCGATg                                                                          <  2:1946094/150‑1 (MQ=255)

                                                                              aaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                                 <  2:4233621/126‑1 (MQ=255)

                                                                              aaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                                 >  1:4233621/1‑126 (MQ=255)

                                                                              aaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATgg                                                                         <  1:2477050/150‑1 (MQ=255)

                                                                               aCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGt                                                                        >  1:1771031/1‑150 (MQ=255)

                                                                                cTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTc                                                                       <  2:3240417/150‑1 (MQ=255)

                                                                                 tCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGCTGGTGACCAGCACCTCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCt                                                                      >  2:2585016/1‑150 (MQ=255)

                                                                                 tCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCt                                                                      >  1:3247647/1‑150 (MQ=255)

                                                                                  cACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCtt                                                                     <  2:2407355/150‑1 (MQ=255)

                                                                                   aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGCTGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGTGACCAGCCGATATAGCCGATGGTCTTc                                                                    <  1:1734941/150‑1 (MQ=255)

                                                                                   aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTc                                                                    >  1:2319312/1‑150 (MQ=255)

                                                                                    ccAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                                   <  1:2634259/150‑1 (MQ=255)

                                                                                     cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcac                                                                  <  2:3999595/150‑1 (MQ=255)

                                                                                     cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcac                                                                  >  1:530109/1‑150 (MQ=255)

                                                                                       gCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcacca                                                                >  1:1847907/1‑150 (MQ=255)

                                                                                           caccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAg                                                            >  2:559797/1‑150 (MQ=255)

                                                                                           caccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAg                                                            >  2:3895488/1‑150 (MQ=255)

                                                                                           caccaGCGCTGAG‑CCAGCAGGGTCAGCGGGATCAGTGGCAGGGTTAACTTGCTGTTCAAGAACAGCATGAAATAGATGGTGACCAGCACATACCAGAACAGCAGCCAGAACAGCGACCAGCCGAAATAGCCGATGGTCATCACCAACAgg                                                           >  1:3754825/1‑150 (MQ=255)

                                                                                              caGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCg                                                         >  2:4528888/1‑150 (MQ=255)

                                                                                                gcgcTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTc                                                                       >  2:2263947/1‑134 (MQ=255)

                                                                                                    tgtgCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACtt                                                   >  1:1932836/4‑150 (MQ=255)

                                                                                                    tGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagca                                                                                                          <  1:425465/95‑1 (MQ=255)

                                                                                                    tGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagca                                                                                                          >  2:425465/1‑95 (MQ=255)

                                                                                                        cccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg                                               >  2:3836733/1‑150 (MQ=255)

                                                                                                        cccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg                                               <  1:1710323/150‑1 (MQ=255)

                                                                                                         ccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCCGCAGCCAGAACAGCGACCCGCCGATATAGCCGATGGTCGTCACCACCAGGCGGGAcgggcggg                                              >  1:4711727/1‑150 (MQ=255)

                                                                                                         ccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCggg                                              >  2:4064617/1‑150 (MQ=255)

                                                                                                           agcagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGATGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCg                                            <  2:955699/150‑1 (MQ=255)

                                                                                                            gcagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGa                                           <  1:3933887/150‑1 (MQ=255)

                                                                                                             cagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGACATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          <  2:3928153/150‑1 (MQ=255)

                                                                                                             cagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          >  2:258074/1‑150 (MQ=255)

                                                                                                             cagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          >  1:876906/1‑150 (MQ=255)

                                                                                                                ggTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGAc                                       <  2:1431921/150‑1 (MQ=255)

                                                                                                                    aGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCt                                   >  2:4474652/1‑150 (MQ=255)

                                                                                                                    aGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTAGCGGGCGATGACGTCt                                   >  2:331780/1‑150 (MQ=255)

                                                                                                                         gATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAg                                                                                         <  2:713803/91‑1 (MQ=255)

                                                                                                                         gATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAg                                                                                         >  1:712709/1‑91 (MQ=255)

                                                                                                                            cAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCAcgcg                           >  1:1560544/1‑150 (MQ=255)

                                                                                                                                 gCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGt                      >  2:84790/1‑150 (MQ=255)

                                                                                                                                   aGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGACCATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTAc                    >  2:4351436/1‑150 (MQ=255)

                                                                                                                                     ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt                  <  2:181191/150‑1 (MQ=255)

                                                                                                                                     ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt                  >  2:1200727/1‑150 (MQ=255)

                                                                                                                                        tttcatGCTGTTCTTGAACAGCAAGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGa               >  1:3499884/6‑150 (MQ=255)

                                                                                                                                            ttGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcg           >  1:266161/1‑150 (MQ=255)

                                                                                                                                            atGCTGTTCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcg           >  1:5008176/2‑150 (MQ=255)

                                                                                                                                            atGCTGTTCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcg           >  2:254576/2‑150 (MQ=255)

                                                                                                                                            atGCTGTTCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcg           >  1:4584951/2‑150 (MQ=255)

                                                                                                                                                tGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       <  2:4042084/150‑1 (MQ=255)

                                                                                                                                                  ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggggg     >  1:2101790/1‑150 (MQ=255)

                                                                                                                                                    cAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGAt   <  2:266161/150‑1 (MQ=255)

                                                                                                                                                     aaGAACATCAAGAAATCGATGGTGACCAGAACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGTGGGATg  <  2:3635703/150‑1 (MQ=255)

                                                                                                                                                     aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATg  >  2:84951/1‑150 (MQ=255)

                                                                                                                                                     aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATg  >  1:4990156/1‑150 (MQ=255)

                                                                                                                                                     aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGACGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATg  >  2:4445925/1‑150 (MQ=255)

AGCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATG  >  NC_002947/966449‑966747

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