BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,401	G→A	6.0%	D349N (GAC→AAC)	PP_1666 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,401	0	G	A	6.0%	516.4 / 14.9	166	D349N (GAC→AAC)	PP_1666	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (78/78); new base A (4/6); total (82/84)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 166 total aligned reads displayed.

TCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGG  >  NC_002947/1863253‑1863545

  gTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACg                                                                                                                                               <  1:5358341/150‑1 (MQ=255)

      aCAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCTCAGCCAGCACGCAAGCCGTCGAACGGCtt                                                                                                                                           <  1:3830303/150‑1 (MQ=255)

        aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGc                                                                                                                                         >  2:3444769/1‑150 (MQ=255)

        aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCACAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTCGAACGGCTTgg                                                                                                                                         <  1:2386479/150‑2 (MQ=255)

            aCGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt                                                                                                                                     >  1:264679/1‑150 (MQ=255)

              gTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTc                                                                                                                                   <  1:1035542/150‑1 (MQ=255)

                gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGtgctggc                                                                                                                                 <  1:4616228/150‑7 (MQ=255)

                ctgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGc                                                                                                                                 >  2:2699338/4‑150 (MQ=255)

                 cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGTCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt                                                                                                                                >  2:1200048/1‑150 (MQ=255)

                 cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGcg                                                                                                                                >  2:1268401/1‑149 (MQ=255)

                 cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt                                                                                                                                <  1:154107/150‑1 (MQ=255)

                  ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTg                                                                                                                               <  2:3160140/150‑1 (MQ=255)

                   cAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTgg                                                                                                                              >  1:5028533/1‑150 (MQ=255)

                   cAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTgg                                                                                                                              >  1:4436873/1‑150 (MQ=255)

                   cAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGACGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTgg                                                                                                                              >  2:4201076/1‑150 (MQ=255)

                    aGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGt                                                                                                                             >  2:2237194/1‑150 (MQ=255)

                    aGCTTGGCCAGGCCAAGCTTCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGt                                                                                                                             >  1:1445518/1‑150 (MQ=255)

                       ttGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAg                                                                                                                          >  1:5235852/1‑150 (MQ=255)

                        tggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAgg                                                                                                                         >  1:2811189/1‑150 (MQ=255)

                         ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGa                                                                                                                        <  1:651314/150‑1 (MQ=255)

                          gcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAg                                                                                                                       >  2:2661866/1‑150 (MQ=255)

                            ctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAgcg                                                                                                                     >  2:1414337/1‑150 (MQ=255)

                               gccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCtt                                                                                                                  <  2:1030030/150‑1 (MQ=255)

                                ccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTg                                                                                                                 >  2:1342193/1‑150 (MQ=255)

                                 cAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                                                                                                >  2:5053152/1‑150 (MQ=255)

                                 cAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                                                                                                >  2:1374786/1‑150 (MQ=255)

                                  aaGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGc                                                                                                                                             >  2:4557942/1‑120 (MQ=255)

                                  aaGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGc                                                                                                                                             <  1:4557942/120‑1 (MQ=255)

                                    gCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGt                                                                                                             <  2:3198830/150‑1 (MQ=255)

                                     cTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtg                                                                                                            >  2:1063925/1‑150 (MQ=255)

                                      tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGTGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           >  2:865034/1‑150 (MQ=255)

                                      tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           >  2:2599281/1‑150 (MQ=255)

                                      tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           <  1:1951001/150‑1 (MQ=255)

                                       gCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgac                                                                                                          >  1:4012338/1‑150 (MQ=255)

                                         aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactg                                                                                                        <  1:408592/150‑1 (MQ=255)

                                         aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactg                                                                                                        <  1:408539/150‑1 (MQ=255)

                                         aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactg                                                                                                        <  1:3908922/150‑1 (MQ=255)

                                                gCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATg                                                                                                   >  2:3479249/1‑148 (MQ=255)

                                                gCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCg                                                                                                 >  1:971117/1‑150 (MQ=255)

                                                    gTCGAGCAGGCACCGGCGGTGCCGGCCATGCAAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGc                                                                                             >  2:4108109/1‑150 (MQ=255)

                                                          cAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATc                                                                                       >  2:2168205/1‑150 (MQ=255)

                                                             gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                                                                                                                                           <  1:4310739/95‑1 (MQ=255)

                                                             gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                                                                                                                                           >  2:4310739/1‑95 (MQ=255)

                                                               aCCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCGTGCCGTGACTGACATGCGTCGCTGGATctggct                                                                                  >  2:3825822/1‑150 (MQ=255)

                                                                 cggcggTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATctggctgg                                                                                <  2:1830409/150‑1 (MQ=255)

                                                                       tGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTg                                                                          <  2:1837823/150‑1 (MQ=255)

                                                                          cGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCAc                                                                       >  1:1413248/1‑150 (MQ=255)

                                                                             ccATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgc                                                                    <  1:1662200/150‑1 (MQ=255)

                                                                             cactGCCGGTACTGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgc                                                                    >  1:1441190/4‑150 (MQ=255)

                                                                              cATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgct                                                                   >  1:962777/1‑150 (MQ=255)

                                                                              cATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgct                                                                   >  1:342542/1‑150 (MQ=255)

                                                                               aTGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgctg                                                                  >  2:1152818/1‑150 (MQ=255)

                                                                                tGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGa                                                                 >  2:2099214/1‑150 (MQ=255)

                                                                                tGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGATTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGa                                                                 <  1:677785/150‑1 (MQ=255)

                                                                                   cAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCg                                                              <  1:1152818/150‑1 (MQ=255)

                                                                                     gTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActg                                                                     <  1:436540/141‑1 (MQ=255)

                                                                                     gTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActg                                                                     >  2:436540/1‑141 (MQ=255)

                                                                                        ccGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtgtg                                                         >  2:1486258/1‑150 (MQ=255)

                                                                                         cGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgc                                                        >  2:3307981/1‑150 (MQ=255)

                                                                                           gCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgctg                                                      <  1:2730147/150‑1 (MQ=255)

                                                                                               tGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCttt                                                  <  2:5018117/150‑1 (MQ=255)

                                                                                               tGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCttt                                                  <  2:4777521/150‑1 (MQ=255)

                                                                                                 gACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTtat                                                >  1:4519012/1‑150 (MQ=255)

                                                                                                 gACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTtat                                                >  1:3264882/1‑150 (MQ=255)

                                                                                                   cgccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAg                                              >  1:768719/1‑150 (MQ=255)

                                                                                                     tgctggCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGcc                                            >  1:5115019/6‑150 (MQ=255)

                                                                                                        cAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGc                                         >  2:3607802/1‑150 (MQ=255)

                                                                                                        cAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGc                                         <  1:2699338/150‑1 (MQ=255)

                                                                                                        cAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGc                                         >  2:2690422/1‑150 (MQ=255)

                                                                                                             gcgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAat                                    >  1:3362276/1‑150 (MQ=255)

                                                                                                              cgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAata                                   <  2:5225576/150‑1 (MQ=255)

                                                                                                                cAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATAtt                                 <  1:4579661/150‑1 (MQ=255)

                                                                                                                cAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATAtt                                 <  2:3040225/150‑1 (MQ=255)

                                                                                                                cAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATAtt                                 <  2:5415852/150‑1 (MQ=255)

                                                                                                                 aGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTc                                <  1:4201076/150‑1 (MQ=255)

                                                                                                                 aGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTc                                <  1:3326688/150‑1 (MQ=255)

                                                                                                                   cccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTATGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCtt                              >  2:2888307/1‑150 (MQ=255)

                                                                                                                   cccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCtt                              <  1:2570514/150‑1 (MQ=255)

                                                                                                                      gcgcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                          >  2:4100977/1‑135 (MQ=255)

                                                                                                                      gcgcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACg                           >  1:1560388/1‑150 (MQ=255)

                                                                                                                          cTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGt                       >  1:2177055/1‑150 (MQ=255)

                                                                                                                           tGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGtt                      <  2:4012338/150‑1 (MQ=255)

                                                                                                                                   cagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCg              <  2:1560388/150‑1 (MQ=255)

                                                                                                                                   cagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCg              >  1:1569673/1‑150 (MQ=255)

                                                                                                                                   cagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCg              <  2:1560359/150‑1 (MQ=255)

                                                                                                                                   cagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCg              <  2:2811189/150‑1 (MQ=255)

                                                                                                                                         gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgct                                                       >  1:4778779/2‑103 (MQ=255)

                                                                                                                                         gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgct                                                       <  2:4778779/102‑1 (MQ=255)

                                                                                                                                         gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                          >  1:3088214/2‑116 (MQ=255)

                                                                                                                                         gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatca        >  2:1740769/2‑150 (MQ=255)

                                                                                                                                         ccAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatca        >  1:4684970/1‑150 (MQ=255)

                                                                                                                                          cAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatcat       >  1:3371798/1‑150 (MQ=255)

                                                                                                                                           aaGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATg      >  1:2326623/1‑150 (MQ=255)

                                                                                                                                              ccGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTg   <  1:1268401/150‑1 (MQ=255)

                                                                                                                                              ccGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTg   <  1:1200048/150‑1 (MQ=255)

                                                                                                                                               cGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTgg  <  2:5235852/150‑1 (MQ=255)

                                                                                                                                                 ttCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGc                                                             >  1:4816834/1‑89 (MQ=255)

                                                                                                                                                 ttCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGc                                                             <  2:4816834/89‑1 (MQ=255)

                                                                                                                                                 ttCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                          <  2:3698730/108‑1 (MQ=255)

                                                                                                                                                 ttCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                          >  1:3698730/1‑108 (MQ=255)

TCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGG  >  NC_002947/1863253‑1863545

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