BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,476	T→C	8.4%	M138T (ATG→ACG) ‡	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,476	0	T	C	8.4%	609.5 / 30.5	203	M138T (ATG→ACG) ‡	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand): ref base T (100/86); new base C (12/5); total (112/91)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 205 total aligned reads displayed.

GAAAGGTGAAGATGTCGAGAAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGC  >  NC_002947/554327‑554625

 aaaGGTGAAGATGTCGAGAAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATg                                                                                                                                                      >  2:1050808/1‑150 (MQ=255)

             gTCGAGAAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGcc                                                                                                                                          >  1:3818077/1‑150 (MQ=255)

             gTCGAGAAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGcc                                                                                                                                          >  1:3336814/1‑150 (MQ=255)

                gagaAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATg                                                                                                                                       >  2:3768189/1‑150 (MQ=255)

                  gaAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGACCGCCGTGTAATGTTCCGTCGCGCCATGaa                                                                                                                                     >  2:4479889/1‑150 (MQ=255)

                  gaAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGaa                                                                                                                                     <  1:5108544/150‑1 (MQ=255)

                         cGTCAGGACCTGTCCAAGCAGTTGTGTGTGCCTGTGCACATCAACTTCGAAGAGATCAGCAAGCCGGATCTCGACGATATGCTGGTTTCGCAGAGCGTAGCTCTGGAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGTAGGGTgag                                                                                                                              <  2:1702790/150‑3 (MQ=255)

                         cGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCg                                                                                                                              >  1:4652729/1‑150 (MQ=255)

                           tCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGt                                                                                                                            >  1:3194735/1‑150 (MQ=255)

                           tCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGt                                                                                                                            >  2:1603679/1‑150 (MQ=255)

                           tCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGt                                                                                                                            >  2:3204718/1‑150 (MQ=255)

                           tCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGt                                                                                                                            >  2:3756197/1‑150 (MQ=255)

                            cAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTa                                                                                                                           <  1:3142760/150‑1 (MQ=255)

                            cAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTa                                                                                                                           <  1:3142656/150‑1 (MQ=255)

                              ggACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACa                                                                                                                         <  1:1177121/150‑1 (MQ=255)

                                aCCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGa                                                                                                                       <  1:1497793/150‑1 (MQ=255)

                                aCCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGa                                                                                                                       <  1:3059511/150‑1 (MQ=255)

                                aCCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGa                                                                                                                       <  2:559190/150‑1 (MQ=255)

                                  cTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAAc                                                                                                                     <  2:144122/150‑1 (MQ=255)

                                   tGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACg                                                                                                                    <  2:103000/150‑1 (MQ=255)

                                      ccAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCa                                                                                                                 <  2:4089985/150‑1 (MQ=255)

                                      ccAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCa                                                                                                                 >  1:1674614/1‑150 (MQ=255)

                                       cAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCACGCCGGAACTCGACGCTATGCTGGTTGCGCAGATCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACTCCAt                                                                                                                <  2:3999632/150‑1 (MQ=255)

                                       cAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCAt                                                                                                                >  1:3557085/1‑150 (MQ=255)

                                       cAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCAt                                                                                                                >  1:3033872/1‑150 (MQ=255)

                                        aaGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATg                                                                                                               >  1:1722938/1‑150 (MQ=255)

                                         aGCAGATGGGTGTGCCTGTGCACATCAACCTCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCATGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGc                                                                                                              <  2:4500398/150‑1 (MQ=255)

                                         aGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGc                                                                                                              >  2:2145815/1‑150 (MQ=255)

                                         aGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGc                                                                                                              <  1:2829368/150‑1 (MQ=255)

                                            aGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTa                                                                                                           <  2:5106266/150‑1 (MQ=255)

                                              aTGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTAtt                                                                                                         >  2:1111423/1‑150 (MQ=255)

                                              aTGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTAtt                                                                                                         >  1:4842002/1‑150 (MQ=255)

                                               tGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTg                                                                                                        <  1:2327704/150‑1 (MQ=255)

                                               tGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTg                                                                                                        <  2:4353263/150‑1 (MQ=255)

                                                gggTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTgg                                                                                                       >  2:4667901/1‑150 (MQ=255)

                                                  gtgtGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTg                                                                                                     >  2:351426/1‑150 (MQ=255)

                                                  gagtgCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTg                                                                                                     >  1:4504044/3‑150 (MQ=255)

                                                     tgCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCa                                                                                                  >  1:1998485/1‑150 (MQ=255)

                                                        cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCTGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                               >  2:2872059/1‑150 (MQ=255)

                                                        cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGTGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                               <  1:2146196/150‑1 (MQ=255)

                                                         tgtgCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAggg                                                                                              >  1:4006177/1‑150 (MQ=255)

                                                           tgCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCa                                                                                            >  2:759/1‑150 (MQ=255)

                                                            gCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCAt                                                                                           <  1:255693/150‑1 (MQ=255)

                                                            gCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCAt                                                                                           <  1:2478206/150‑1 (MQ=255)

                                                             cacaTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATc                                                                                          >  1:1274372/1‑150 (MQ=255)

                                                              acaTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCa                                                                                         >  2:2318747/1‑150 (MQ=255)

                                                                aTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAg                                                                                       >  1:99455/1‑150 (MQ=255)

                                                                aTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAg                                                                                       >  1:564618/1‑150 (MQ=255)

                                                                aTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTAGCGCATAGCGTAGATCAGCAGCAGGAACGCCGTGTAATTTTCCGTCCCGCCATGAAGCGTCCGCTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAg                                                                                       >  1:3678635/1‑150 (MQ=255)

                                                                 tCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGa                                                                                      <  1:3293506/150‑1 (MQ=255)

                                                                 tCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGa                                                                                      >  1:1900506/1‑150 (MQ=255)

                                                                  cAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGAt                                                                                     >  1:2049292/1‑150 (MQ=255)

                                                                      aTCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAg                                                                                 <  1:4782605/150‑1 (MQ=255)

                                                                         gAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTg                                                                              >  1:3944313/1‑150 (MQ=255)

                                                                         gAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTg                                                                              >  2:2085076/1‑150 (MQ=255)

                                                                         gAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTg                                                                              >  1:836752/1‑150 (MQ=255)

                                                                          aaGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGt                                                                                                      <  1:1103677/125‑1 (MQ=255)

                                                                             agaTCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCg                                                                          >  1:713141/1‑150 (MQ=255)

                                                                              gaTCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCgg                                                                         <  1:4070180/150‑1 (MQ=255)

                                                                               aTCCGCAAGCCGGAACTCGACGCTCTGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGACCGCCATGCGGATTGGTGCCAAGGGCAGCAAGATCCAGGTGAGCggg                                                                        >  2:2701599/1‑149 (MQ=255)

                                                                                tCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgtc                                                                       <  1:4547289/150‑1 (MQ=255)

                                                                                  cGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgtcgt                                                                     <  1:4270452/150‑1 (MQ=255)

                                                                                       gCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcgg                                                                <  2:4255044/150‑1 (MQ=255)

                                                                                         cGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcggcg                                                              <  1:2523216/150‑1 (MQ=255)

                                                                                           gAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGt                                                            >  2:3446098/1‑150 (MQ=255)

                                                                                           gAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGt                                                            >  2:2803032/1‑150 (MQ=255)

                                                                                               tCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTg                                                        >  1:2316487/1‑150 (MQ=255)

                                                                                                cGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTga                                                       <  1:4189731/150‑1 (MQ=255)

                                                                                                   cGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGAt                                                    <  2:2937321/150‑1 (MQ=255)

                                                                                                    gCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGAtt                                                   <  2:2384581/150‑1 (MQ=255)

                                                                                                    gCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGAtt                                                   >  2:1912871/1‑150 (MQ=255)

                                                                                                     cTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTg                                                  >  2:1309998/1‑150 (MQ=255)

                                                                                                      tATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGc                                                 <  2:4765579/150‑1 (MQ=255)

                                                                                                            ggTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTAc                                           <  1:3886931/150‑1 (MQ=255)

                                                                                                             gTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTAcc                                          >  1:1641290/1‑150 (MQ=255)

                                                                                                               tGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGa                                        <  1:562852/150‑1 (MQ=255)

                                                                                                               tGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGa                                        <  2:3818852/150‑1 (MQ=255)

                                                                                                                  gcAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTg                                     <  2:2924633/150‑1 (MQ=255)

                                                                                                                        gatcgCTCAGCAGCTCGAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTTCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATcg                               >  2:3519735/5‑150 (MQ=255)

                                                                                                                         gTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATcgc                              <  2:502149/150‑1 (MQ=255)

                                                                                                                          tAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATcgcg                             >  1:2267249/1‑150 (MQ=255)

                                                                                                                          tAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATcgcg                             >  1:2267106/1‑150 (MQ=255)

                                                                                                                            gCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGaa                           <  1:4147372/150‑1 (MQ=255)

                                                                                                                            gCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGaa                           <  1:1003930/150‑1 (MQ=255)

                                                                                                                            gCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGaa                           <  1:3385929/150‑1 (MQ=255)

                                                                                                                            gCTCAGCAGCTCGAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGaa                           >  2:1446298/1‑150 (MQ=255)

                                                                                                                             cTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAg                          <  1:181782/150‑1 (MQ=255)

                                                                                                                             cTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAg                          >  1:1803902/1‑150 (MQ=255)

                                                                                                                              tCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGTGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATAGCGAAgg                         >  2:4911358/1‑150 (MQ=255)

                                                                                                                                  cagcTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCgt                     <  1:3768189/150‑1 (MQ=255)

                                                                                                                                    gcTCGAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCgcgt                   >  2:2350115/1‑150 (MQ=255)

                                                                                                                                      tcgAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCgcgttc                 >  1:3491150/3‑148 (MQ=255)

                                                                                                                                        gAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTccg               >  1:1342420/1‑149 (MQ=255)

                                                                                                                                        gAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCct               >  2:727764/1‑150 (MQ=255)

                                                                                                                                         aaCGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTccgc              <  1:1571251/150‑3 (MQ=255)

                                                                                                                                           cGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTg            <  2:1254350/150‑1 (MQ=255)

                                                                                                                                           cGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTg            <  2:271519/150‑1 (MQ=255)

                                                                                                                                                tgtAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCacac       >  2:4585211/1‑150 (MQ=255)

                                                                                                                                                 gtAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACAcc      <  1:964702/150‑1 (MQ=255)

                                                                                                                                                    aTGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTg   >  1:3586440/1‑150 (MQ=255)

GAAAGGTGAAGATGTCGAGAAGCTGCGTCAGGACCTGACCAAGCAGATGGGTGTGCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGC  >  NC_002947/554327‑554625

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