BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,131,657	T→G	5.3%	D950A (GAC→GCC)	gcvP‑I ←	glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,131,657	0	T	G	5.3%	743.8 / 17.9	219	D950A (GAC→GCC)	gcvP‑I	glycine dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (110/96); new base G (7/6); total (117/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

TGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCC  >  NC_002947/1131511‑1131804

tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTAGCCGCGAAAGCGGAGGTATATGCACCAGCGCATTCGCGGGTAAAACCCCCCCAAAACAGAGCCCCCCACCCCAGCCAGGTGGCGGGGGGGGCGAGAAGGCGGccg                                                                                                                                                   >  1:5054913/1‑149 (MQ=255)

tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTAACCCCGAATGAGGAGGTATTGCCACCACCGCATTAGCGGGTAAACCCGCAACGAAACAGAGCCCGCCCCCCAAGCCCGGTGGCGGGGGGGGCTATCAGGCGTccg                                                                                                                                                   >  2:722390/1‑149 (MQ=255)

tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTAACCCCCCAAGCGGCGGGATAGCCACCCTCGCATTACCGGGTAAAACCCCACCAAAACAGCACCCCCCACCCAACCCACGTGGCCGGGCGGGCGAAAAGGCGGccg                                                                                                                                                   >  2:415310/1‑149 (MQ=255)

  cctggTTGTTGGGTGGACATTGAGCCCGGTCGGCGG‑GGTTTTCCCCGGGAAGGCGGTGGTATCCCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                                                                                 <  1:3686212/147‑1 (MQ=255)

  cccggTTGTGGGGGGGACTTTGCGCCCTGTAGGCGC‑GGGTTTCCCCGCGACTGCGGAGGTATCCCCACCATCGCATTCGCGGGAAAACCCGCTCCTACACCGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                                                                                 <  1:3593359/147‑1 (MQ=255)

  ccAGGTTGTTGGGTGGACATTGAGCCCTGTCGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCTTTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                                                                                 <  1:2100492/150‑1 (MQ=255)

  ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTGACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGTCGTCCTGa                                                                                                                                                 <  2:4209346/150‑1 (MQ=255)

    aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGTTTTACCCGCGATTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa                                                                                                                                               <  1:556475/150‑1 (MQ=255)

     ggTTGTGGGGTGGAAATTGCGCCGTGTCGGCGC‑GGGTTTACCCGCTAATGCGGTGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                                                                                              <  1:6359273/150‑1 (MQ=38)

       ttgttgGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGcc                                                                                                                                            <  1:6277358/150‑2 (MQ=255)

        tgttgGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                                                                                           <  1:6706409/150‑1 (MQ=255)

        tgttgGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                                                                                           <  2:4262454/150‑1 (MQ=255)

        tgttgGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                                                                                           <  2:4781627/150‑1 (MQ=255)

            gggTGGACATTGAGCCCTGTAGGCGC‑GGGTTTAACCCCGAATGCGGTGGTAAAACCACCAAAGCAATGGCGGGTAAAACCGCTCCAACACAGAGCCCCCCACCCAAGCCACGTGGCGGGGCGGGCGAACAGGCGTACAGATAGCTCacga                                                                                                                                       >  2:1912372/1‑150 (MQ=255)

                    aTTGAGCCCTGTAGGCGCGGGGTTTACCC‑CGCAAGCGGGGGTTTAGCCCCCCACGCAGTAGCGGGGAAACCCGCGCCAACACACGGCCCCCCCCCCCACCCCGGGGGGGGGGGGGGCGATCCGGCGGCCGGCTAGCTCGCGATCCACggg                                                                                                                               >  2:1834892/1‑150 (MQ=255)

                       gcgCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGTGTTATAGCCACCATCGCATTCGCGGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAg                                                                                                                            <  1:4272147/148‑1 (MQ=255)

                        cgCCCTGTAGGCGC‑GGGTTGACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGTGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                                                                                           <  2:4739185/149‑1 (MQ=255)

                        aGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                                                                                           <  2:2916160/150‑1 (MQ=255)

                         gCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                                                                                          <  2:2608642/150‑1 (MQ=255)

                         gCCCTGTAGGCGC‑GGGTTTACCCGCGAAGGCGGAGTTATCGCCACCATCGCATGCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                                                                                          <  1:4302699/150‑1 (MQ=255)

                              gTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGcaca                                                                                                                     <  2:5815790/150‑1 (MQ=255)

                               tAGGCGC‑GGGTTTACCCGCGATTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac                                                                                                                    <  1:5597855/150‑1 (MQ=255)

                               cattcgc‑gGGTAAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac                                                                                                                    >  1:4641329/5‑150 (MQ=255)

                                         gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                                                                                           <  1:2231617/150‑1 (MQ=255)

                                          tttAACCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCACGCCACCCAAGCCAGGTGGCGGGGAGGGCTCTCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCgg                                                                                                          >  2:2426406/1‑150 (MQ=255)

                                            tACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                                                                                        <  2:1113574/150‑1 (MQ=255)

                                               ccGCGAATGGGGAGGTACAGCCACCAGCGCATTGGGGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                                                                                     <  1:1716703/150‑1 (MQ=255)

                                               ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                                                                                          >  1:2436333/1‑113 (MQ=255)

                                               ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                                                                                     <  2:4869964/150‑1 (MQ=255)

                                               ccGCGAATGCGGAGGTATAGCCACCATCGCAGTCGCGGGTAAACCCGCTACTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                                                                                     >  1:6624430/1‑150 (MQ=255)

                                                cgcgACTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGATGGGCAGGCGCACACCAGGTTGCGGTCGCCg                                                                                                    <  1:1373349/150‑1 (MQ=255)

                                                cgcgAATGCGGAGGTATCGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGATGGGCAGGCGCACACCAGGTTGCGGTCGCCg                                                                                                    <  1:4818805/150‑1 (MQ=255)

                                                cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                                                                                                    >  1:1122975/1‑150 (MQ=255)

                                                 gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                                                                                                  >  2:2292561/1‑135 (MQ=255)

                                                 gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa                                                                                                   >  2:2298152/1‑150 (MQ=255)

                                                 gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa                                                                                                   >  1:302875/1‑150 (MQ=255)

                                                 gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa                                                                                                   >  1:2731076/1‑150 (MQ=255)

                                                 gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCACTCGATCGACGGGCAGGCGCACACCAGGTAGCGGTCGCCGa                                                                                                   <  1:2257789/150‑1 (MQ=255)

                                                   gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac                                                                                                 >  1:755121/1‑150 (MQ=255)

                                                       gCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTTCGGTCGCCGAACACGt                                                                                             >  2:4899879/1‑150 (MQ=255)

                                                         ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                                                                                           >  2:2794830/1‑150 (MQ=255)

                                                           aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                                                                                         <  2:1844952/150‑1 (MQ=255)

                                                             gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGa                                                                                       >  1:3760981/1‑150 (MQ=255)

                                                               ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                                                                                     >  1:3476381/1‑150 (MQ=255)

                                                                     accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGTTAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg                                                                               >  1:5991109/1‑150 (MQ=255)

                                                                     accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg                                                                               >  1:5393381/1‑150 (MQ=255)

                                                                       caccGCATTCGCGGGTAAACCCGCTCCTACACCGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc                                                                             <  1:828558/147‑1 (MQ=255)

                                                                         tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg                                                                           <  1:2349740/150‑1 (MQ=255)

                                                                         tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg                                                                           >  1:4724407/1‑150 (MQ=255)

                                                                         tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg                                                                           >  2:6776777/1‑150 (MQ=255)

                                                                           gCATTCGCGGGTAAACCCGCTCCTACACAGCGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc                                                                         <  1:4143090/150‑1 (MQ=255)

                                                                            cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg                                                                        >  1:5425624/1‑150 (MQ=255)

                                                                                cgcgGGTAACCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa                                                                    <  1:1316776/150‑1 (MQ=255)

                                                                                cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa                                                                    <  2:5737694/150‑1 (MQ=255)

                                                                                 gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg                                                                   >  2:4461169/1‑150 (MQ=255)

                                                                                 gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg                                                                   >  2:5905811/1‑150 (MQ=255)

                                                                                        aaCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGc                                                            <  2:3056717/150‑1 (MQ=255)

                                                                                         aCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGcc                                                           >  1:6329221/1‑150 (MQ=255)

                                                                                           ccGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCtt                                                         >  2:3210000/1‑150 (MQ=255)

                                                                                               tCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTccac                                                     >  2:5317054/1‑150 (MQ=255)

                                                                                               tCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTccac                                                     >  2:2211509/1‑150 (MQ=255)

                                                                                               tCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTccac                                                     >  1:1677286/1‑150 (MQ=255)

                                                                                                  tACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAg                                                  >  1:389586/1‑150 (MQ=255)

                                                                                                  tACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAg                                                  >  2:1610301/1‑150 (MQ=255)

                                                                                                    cacaGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGc                                                >  1:1787197/1‑150 (MQ=255)

                                                                                                     acaGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCt                                               >  1:321734/1‑150 (MQ=255)

                                                                                                        tcaatgtccACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGc                                            >  1:3668200/8‑150 (MQ=255)

                                                                                                         tccccGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac                                                                                                 >  2:1772391/3‑96 (MQ=255)

                                                                                                           cccGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑g                                        >  1:1168900/1‑150 (MQ=255)

                                                                                                           cccGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑g                                        >  1:1862563/1‑150 (MQ=255)

                                                                                                            ccGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑gg                                       <  2:332494/150‑1 (MQ=255)

                                                                                                             cGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑ggg                                      <  1:4534236/150‑1 (MQ=255)

                                                                                                             cGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑ggg                                      >  1:3106360/1‑150 (MQ=255)

                                                                                                               ccACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑ggggg                                    >  1:14391/1‑150 (MQ=255)

                                                                                                                 aCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGta                                  >  1:2765852/1‑150 (MQ=255)

                                                                                                                   ccAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGtata                                >  2:373638/1‑150 (MQ=255)

                                                                                                                    cAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATAc                               >  2:3128886/1‑150 (MQ=255)

                                                                                                                     aaGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATAcc                              >  2:2225341/1‑150 (MQ=255)

                                                                                                                     aaGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATAcc                              >  2:4058650/1‑150 (MQ=255)

                                                                                                                         cAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATCCCGCCt                          <  1:716635/150‑1 (MQ=255)

                                                                                                                         cAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCt                          >  2:3355010/1‑150 (MQ=255)

                                                                                                                         cAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCt                          <  1:604818/150‑1 (MQ=255)

                                                                                                                             tGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTc                      <  1:5913160/150‑1 (MQ=255)

                                                                                                                               gCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTcgc                    >  1:5494906/1‑150 (MQ=255)

                                                                                                                                 ggggAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCgg                  >  2:6645646/1‑150 (MQ=255)

                                                                                                                                  gggAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGc                 >  2:2801801/1‑150 (MQ=255)

                                                                                                                                     agagCTATCAGGACGCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGt              >  1:2190901/4‑150 (MQ=255)

                                                                                                                                      gggCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGtc             >  2:13578/1‑149 (MQ=255)

                                                                                                                                      gggCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTa             <  1:2794830/150‑1 (MQ=255)

                                                                                                                                      gggCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTa             <  2:6864468/150‑1 (MQ=255)

                                                                                                                                       ggCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAg            >  2:2089639/1‑150 (MQ=255)

                                                                                                                                           aTCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACCCGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCAt        >  2:3707587/1‑150 (MQ=255)

                                                                                                                                           aTCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCAt        <  1:524655/150‑1 (MQ=255)

                                                                                                                                           aTCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCAt        <  1:6228584/150‑1 (MQ=255)

                                                                                                                                            tCAGGCGTCCTGATAGCTCTCGATCGACGGGCCGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCCGGCGGCCAGTACTTGCCGTCCACCAGGCTCGCCAGGGGGGTATACCGCCTGCTCGCGGCTGTAGCCAt        <  2:6763958/150‑1 (MQ=255)

                                                                                                                                            tCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATg       >  2:6326897/1‑150 (MQ=255)

                                                                                                                                             cAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATgg      <  1:1565845/150‑1 (MQ=255)

                                                                                                                                              aGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATggg     <  1:6768351/150‑1 (MQ=255)

                                                                                                                                                gCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTc   <  1:6183253/150‑1 (MQ=255)

                                                                                                                                                 cGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTcc  <  1:2879829/150‑1 (MQ=255)

TGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGC‑GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCCTTCCACCAGGCTCGCCA‑GGGGGTATACCGCCTGCTCGCGGCTGTAGCCATGGGTCC  >  NC_002947/1131511‑1131804

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