BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,816	A→G	9.5%	K58K (AAA→AAG)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,816	0	A	G	9.5%	1056.6 / 69.0	360	K58K (AAA→AAG)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base A (145/180); new base G (20/14); total (165/195)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 372 total aligned reads displayed.

GGTGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGC  >  NC_002947/1549667‑1549965

  tGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGa                                                                                                                                                     >  1:991219/1‑150 (MQ=255)

   gATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGaa                                                                                                                                                    <  1:5647980/150‑1 (MQ=255)

   gATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGaa                                                                                                                                                    >  2:2422152/1‑150 (MQ=255)

     ttCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAAt                                                                                                                                                  >  1:4795340/1‑150 (MQ=255)

     ttCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAAt                                                                                                                                                  >  2:1071663/1‑150 (MQ=255)

      tCCGCTCGAAACAAAATGCTGGTTGGTGTCAACGATCTGGCTGACGCGGAAAAAGCGACCCTGCGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGCACGGCGTTTCCGTCGCCAAAGAAATc                                                                                                                                                 <  2:7554467/150‑1 (MQ=255)

        cGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTCCCGTCGCCAAAGAAATCGa                                                                                                                                               <  2:6229894/150‑1 (MQ=255)

              tAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCATGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGat                                                                                                                                         >  1:8099473/1‑149 (MQ=255)

               aaGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTATAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCTAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGaaa                                                                                                                                        <  2:3750814/150‑1 (MQ=255)

               aaGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGaaa                                                                                                                                        >  2:7823822/1‑150 (MQ=255)

               aaGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGaaa                                                                                                                                        <  2:6644914/150‑1 (MQ=255)

               aaGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGaaa                                                                                                                                        >  2:7786658/1‑150 (MQ=255)

                aGAAAATGCTGGTTGGTGTCACCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAg                                                                                                                                       <  2:1627767/150‑1 (MQ=255)

                aGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAg                                                                                                                                       >  1:8184803/1‑150 (MQ=255)

                 gAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGa                                                                                                                                      >  1:3224800/1‑150 (MQ=255)

                    aaTGCTGGTTGGTGTCAAGGTTCTGGCTTACGCTGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGc                                                                                                                                   <  1:2203685/150‑1 (MQ=255)

                      tGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCtt                                                                                                                                 >  1:2276774/1‑150 (MQ=255)

                         tggttggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCg                                                                                                                              >  2:5795038/1‑150 (MQ=255)

                           gttggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGaa                                                                                                                            <  2:4273077/150‑1 (MQ=255)

                           gttggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGaa                                                                                                                            >  2:3264616/1‑150 (MQ=255)

                            ttggtGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGaaa                                                                                                                           >  1:5965038/1‑150 (MQ=255)

                               gtgtCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACa                                                                                                                        >  2:1787049/1‑150 (MQ=255)

                               gtgtCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACa                                                                                                                        >  2:6832878/1‑150 (MQ=255)

                               gtgtCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACa                                                                                                                        >  2:7277128/1‑150 (MQ=255)

                                 gtCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATg                                                                                                                      <  1:2817960/150‑1 (MQ=255)

                                  tCACCGTTCTGGCTGACGCGTTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCACCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                     <  1:5279693/150‑1 (MQ=255)

                                  tCAACTTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                     <  2:940018/150‑1 (MQ=255)

                                  tCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                     <  2:5721226/150‑1 (MQ=255)

                                  tCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                     <  1:2417150/150‑1 (MQ=255)

                                  tCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATgg                                                                                                                     <  1:7613369/150‑1 (MQ=255)

                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGgc                                                                                                                   >  2:6544103/1‑150 (MQ=255)

                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGgc                                                                                                                   <  1:2055613/150‑1 (MQ=255)

                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGgc                                                                                                                   >  1:4483883/1‑150 (MQ=255)

                                    aaCGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCGGTCGAAAACATGGgc                                                                                                                   >  1:718387/1‑150 (MQ=255)

                                       gTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGcc                                                                                                                >  2:4663998/1‑150 (MQ=255)

                                            ggCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCt                                                                                                           >  1:7416238/1‑150 (MQ=255)

                                               tGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTCCTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTGTCGAAAACATGGGCGCCCAGCTGGt                                                                                                        >  2:1435608/1‑150 (MQ=255)

                                                gACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTc                                                                                                       >  1:1679498/1‑150 (MQ=255)

                                                  cgcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCaa                                                                                                     <  1:5911193/150‑1 (MQ=255)

                                                  cgcgGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCaa                                                                                                     <  2:2748822/150‑1 (MQ=255)

                                                     ggTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGa                                                                                                  <  1:1068850/150‑1 (MQ=255)

                                                     ggTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGa                                                                                                  <  1:7093295/150‑1 (MQ=255)

                                                      gTAATAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGaa                                                                                                 <  1:8365601/150‑1 (MQ=255)

                                                      gTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGaa                                                                                                 <  2:8723361/150‑1 (MQ=255)

                                                       tAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAg                                                                                                <  1:6489119/150‑1 (MQ=255)

                                                         aaaGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGtt                                                                                              >  2:266328/1‑150 (MQ=255)

                                                           aGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGc                                                                                            >  1:1921103/1‑150 (MQ=255)

                                                            gCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCt                                                                                           >  2:7092095/1‑150 (MQ=255)

                                                               aCCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTcc                                                                                        <  1:439938/150‑1 (MQ=255)

                                                                 ccTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCCCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCaa                                                                                      >  1:8751356/1‑150 (MQ=255)

                                                                      gCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCa                                                                                 >  2:7880973/1‑150 (MQ=255)

                                                                      gCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCa                                                                                 >  1:2007896/1‑150 (MQ=255)

                                                                          gAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAAGATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacga                                                                             <  2:8184803/150‑1 (MQ=255)

                                                                           aaaGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgac                                                                            >  2:5467243/1‑150 (MQ=255)

                                                                                 cGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAacgacg                                                                           >  1:8707930/1‑145 (MQ=255)

                                                                                      cGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGa                                                                 <  2:8498258/150‑1 (MQ=255)

                                                                                         gttACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACgg                                                              <  1:5570938/148‑1 (MQ=255)

                                                                                            aCTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcac                                                           >  1:4046825/1‑150 (MQ=255)

                                                                                            aCTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcac                                                           <  2:1575597/150‑1 (MQ=255)

                                                                                            aCTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcac                                                           <  1:6850464/150‑1 (MQ=255)

                                                                                                gCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcaccacc                                                       >  2:2574393/1‑150 (MQ=255)

                                                                                                gCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGcaccacc                                                       >  1:8676202/1‑150 (MQ=255)

                                                                                                   aaGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCaccaccacc                                                    >  2:2290210/1‑150 (MQ=255)

                                                                                                    agagCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCg                                                   >  1:5735667/1‑150 (MQ=255)

                                                                                                      tccTTCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGcg                                                 >  2:8657728/3‑149 (MQ=255)

                                                                                                         ttCGGCGCTCCACGCATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGAcc                                              <  2:4897482/150‑1 (MQ=255)

                                                                                                         ttCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAcc                                              >  2:8896725/1‑150 (MQ=255)

                                                                                                             gcgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTc                                          >  1:8511707/1‑150 (MQ=255)

                                                                                                             gcgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTc                                          <  2:7322791/150‑1 (MQ=255)

                                                                                                             gcgcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTc                                          <  1:1350752/150‑1 (MQ=255)

                                                                                                               gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTg                                        <  1:391720/150‑1 (MQ=255)

                                                                                                               gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTg                                        <  1:1440568/150‑1 (MQ=255)

                                                                                                               gcgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTg                                        <  1:3004318/150‑1 (MQ=255)

                                                                                                                cgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                       <  2:5281161/150‑1 (MQ=255)

                                                                                                                cgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                       <  1:4219214/150‑1 (MQ=255)

                                                                                                                cgcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTgg                                       <  1:4168424/150‑1 (MQ=255)

                                                                                                                      gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                 <  2:3824054/149‑1 (MQ=255)

                                                                                                                       cATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGc                                >  1:7878930/1‑150 (MQ=255)

                                                                                                                       cATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGc                                >  2:638454/1‑150 (MQ=255)

                                                                                                                       cATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGc                                >  1:8414004/1‑150 (MQ=255)

                                                                                                                           aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                            >  1:8919293/1‑150 (MQ=255)

                                                                                                                           aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                            >  1:6338255/1‑150 (MQ=255)

                                                                                                                             cAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGt                          >  1:4561629/1‑150 (MQ=255)

                                                                                                                             cAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                          >  2:3442859/1‑150 (MQ=255)

                                                                                                                             cAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                          >  2:2929410/1‑150 (MQ=255)

                                                                                                                                    ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                   <  2:4087666/150‑1 (MQ=255)

                                                                                                                                    ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                   <  1:3538322/150‑1 (MQ=255)

                                                                                                                                         ttCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt              <  2:1184586/150‑1 (MQ=255)

                                                                                                                                           ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTATCGAAAACATGCGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGCCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAACGCCTGa            <  2:8035503/150‑1 (MQ=255)

                                                                                                                                             gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa          <  2:4483883/150‑1 (MQ=255)

                                                                                                                                                gCCAAAGAAATCGCGCTGAAAGACGCTTTCGAAAACCTGGGCGCCCCGCTGGTCAAGGACGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGcc       <  2:3570140/150‑1 (MQ=255)

                                                                                                                                                 ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTg                                                                                             >  2:4349809/1‑63 (MQ=255)

                                                                                                                                                  cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt     >  2:279513/1‑150 (MQ=255)

                                                                                                                                                  cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt     >  1:4970929/1‑150 (MQ=255)

                                                                                                                                                  cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt     <  1:663989/150‑1 (MQ=255)

                                                                                                                                                   aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc    >  2:2169093/1‑150 (MQ=255)

                                                                                                                                                   aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc    >  2:3676190/1‑150 (MQ=255)

                                                                                                                                                    aaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGa                    <  2:8363192/133‑1 (MQ=255)

                                                                                                                                                     aGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCCGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCCTCGTCAACGAAGGCCTGAGAGCCGTCGc  >  1:2691113/1‑150 (MQ=255)

                                                                                                                                                     aGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGc  >  1:5379934/1‑150 (MQ=255)

GGTGATTCCGCTCGTAAGAAAATGCTGGTTGGTGTCAACGTTCTGGCTGACGCGGTAAAAGCGACCCTGGGCCCGAAAGGCCGTAACGTGGTACTGGCCAAGAGCTTCGGCGCGCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGC  >  NC_002947/1549667‑1549965

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