BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,009,719	C→T	A230T (GCG→ACG)	metR ←	MetR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,719	0	C	T	92.8%	1255.8 / 42.6	400	A230T (GCG→ACG)	metR	MetR
Reads supporting (aligned to +/- strand): ref base C (19/10); new base T (186/185); total (205/195)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

Only 100 of 413 total aligned reads displayed.

CAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCA  >  NC_000913/4009627‑4009811

caAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCcgcgc                                                                                              >  1:6581789/1‑93 (MQ=255)

  aGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTgg                                                                                            >  1:780812/1‑93 (MQ=255)

   gCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGc                                                                                           <  2:4791779/93‑1 (MQ=255)

      ttCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAAc                                                                                        <  2:4276229/93‑1 (MQ=255)

        cGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCa                                                                                      <  1:386507/93‑1 (MQ=255)

           ccAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCt                                                                                   >  1:8901529/1‑93 (MQ=255)

           ccAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCt                                                                                   >  2:3850151/1‑93 (MQ=255)

            cAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTg                                                                                  >  1:7436713/1‑93 (MQ=255)

              gggTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGaa                                                                                <  1:7626888/93‑1 (MQ=255)

              gggTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGaa                                                                                <  2:4336795/93‑1 (MQ=255)

              gggTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGaa                                                                                <  2:5131141/93‑1 (MQ=255)

                gTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATc                                                                              <  1:7016852/93‑1 (MQ=255)

                 ttttGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATca                                                                             >  1:9757868/1‑93 (MQ=255)

                 ttttGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATca                                                                             >  1:8755354/1‑93 (MQ=255)

                 ttttGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATca                                                                             >  2:1176288/1‑93 (MQ=255)

                 gtttGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATca                                                                             <  1:1690167/92‑1 (MQ=255)

                   ttGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATcaac                                                                           >  1:2113257/1‑93 (MQ=255)

                    tGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATcaaca                                                                          <  2:1799536/93‑1 (MQ=255)

                     ggTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATcaacaa                                                                         <  1:3415826/93‑1 (MQ=255)

                     ggTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATcaacaa                                                                         >  2:2993647/1‑93 (MQ=255)

                       tCACCACCAGACCCTGGGGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATa                                                                       <  2:7760347/93‑1 (MQ=255)

                       tCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATa                                                                       <  1:5649241/93‑1 (MQ=255)

                          ccaccaGACCCTGGCGCTCAAAGCTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAgg                                                                    <  1:8242603/93‑1 (MQ=255)

                          ccaccaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAgg                                                                    <  1:3675441/93‑1 (MQ=255)

                           caccaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGgt                                                                   >  2:6600242/1‑93 (MQ=255)

                            accaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGgtg                                                                  >  1:9693037/1‑93 (MQ=255)

                            accaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGgtg                                                                  >  1:411458/1‑93 (MQ=255)

                            accaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGgtg                                                                  >  1:3119348/1‑93 (MQ=255)

                            accaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGgtg                                                                  >  1:10004611/1‑93 (MQ=255)

                             ccaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGgtgt                                                                 <  2:4001595/93‑1 (MQ=255)

                               aGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTa                                                               >  1:5899699/1‑93 (MQ=255)

                                 aCCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATc                                                             >  2:2022890/1‑93 (MQ=255)

                                 aCCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATc                                                             <  1:5549577/93‑1 (MQ=255)

                                 aCCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATc                                                             <  2:2726942/93‑1 (MQ=255)

                                   ccTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGa                                                           <  1:6781074/93‑1 (MQ=255)

                                   ccTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGa                                                           <  1:2299981/93‑1 (MQ=255)

                                   ccTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGa                                                           <  1:7183022/93‑1 (MQ=255)

                                    cTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGAc                                                          >  1:9208293/1‑93 (MQ=255)

                                    cTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGAc                                                          <  2:6942114/93‑1 (MQ=255)

                                    cTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGAc                                                          >  2:6348582/1‑93 (MQ=255)

                                     tGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACg                                                         >  2:715786/1‑93 (MQ=255)

                                     tGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACg                                                         >  2:9293201/1‑93 (MQ=255)

                                         gcTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCttt                                                     <  2:1920193/93‑1 (MQ=255)

                                         gcTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCttt                                                     >  1:4560120/1‑93 (MQ=255)

                                           tCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTc                                                   <  2:9822042/93‑1 (MQ=255)

                                           tCAAAACTCTCTACTACCCAATGCGGTAGAGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTc                                                   >  1:8064205/1‑93 (MQ=255)

                                           ccAAAACTCTCTCCTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTc                                                   <  2:2343626/92‑1 (MQ=255)

                                            cAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCa                                                  <  2:5472029/93‑1 (MQ=255)

                                            cAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCa                                                  >  1:5202578/1‑93 (MQ=255)

                                             aaaaCTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAg                                                 >  2:8066296/1‑93 (MQ=255)

                                              aaaCTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGt                                                >  1:6016794/1‑93 (MQ=255)

                                               aaCTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTg                                               >  1:3989260/1‑93 (MQ=255)

                                                  tctctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACg                                            >  2:9537634/1‑93 (MQ=255)

                                                  tctctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACg                                            >  1:3231918/1‑93 (MQ=255)

                                                   ctctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACgg                                           <  1:2678711/93‑1 (MQ=255)

                                                    tctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACggg                                          <  2:8400667/93‑1 (MQ=255)

                                                    tctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACggg                                          >  2:5946811/1‑93 (MQ=255)

                                                     ctactaCCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGc                                         >  1:9619108/1‑93 (MQ=255)

                                                      tactacCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCt                                        <  2:7288009/93‑1 (MQ=255)

                                                      tactacCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCt                                        >  1:5281312/1‑93 (MQ=255)

                                                       actacCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTg                                       >  2:3374301/1‑93 (MQ=255)

                                                        ctacCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGa                                      <  1:358147/93‑1 (MQ=255)

                                                        ctacCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGa                                      <  2:9109714/93‑1 (MQ=255)

                                                         tacCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGAc                                     >  2:6410289/1‑93 (MQ=255)

                                                         tacCCAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGAc                                     >  1:1602165/1‑93 (MQ=255)

                                                            ccAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGcc                                  >  2:8108570/1‑93 (MQ=255)

                                                            ccAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGcc                                  <  1:1934808/93‑1 (MQ=255)

                                                             cAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCt                                 >  2:2112953/1‑93 (MQ=255)

                                                             cAATGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCt                                 <  1:6212466/93‑1 (MQ=255)

                                                              aaTGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGTCtt                                >  1:4970037/1‑92 (MQ=255)

                                                               aTGCGGTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGc                               >  2:9207966/1‑93 (MQ=255)

                                                                    gTAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCt                          >  2:1455924/1‑93 (MQ=255)

                                                                     tAGCGCGGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTg                         <  2:5732974/93‑1 (MQ=255)

                                                                     tAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTg                         <  2:10011489/93‑1 (MQ=255)

                                                                          cgGCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTgaagaa                    <  1:6937006/93‑1 (MQ=255)

                                                                            gCAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGaaaa                  >  2:8183866/1‑93 (MQ=255)

                                                                            gCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGaaaa                  <  2:8755353/93‑1 (MQ=255)

                                                                             cAATACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAAt                 <  1:9343217/93‑1 (MQ=255)

                                                                              aaTACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATg                >  2:4501988/1‑93 (MQ=255)

                                                                              aaTACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATg                <  2:9374358/93‑1 (MQ=255)

                                                                              aaTACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATg                <  1:9871602/93‑1 (MQ=255)

                                                                               acaCCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAACATgc               >  1:7357812/3‑93 (MQ=255)

                                                                               aTACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATgc               <  1:8477094/93‑1 (MQ=255)

                                                                               aTACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATgc               >  2:1717910/1‑93 (MQ=255)

                                                                                tACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATgcc              <  1:629316/93‑1 (MQ=255)

                                                                                tACCCATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATgcc              <  2:2592456/93‑1 (MQ=255)

                                                                                  cccATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATgccgc            >  1:4819708/1‑93 (MQ=255)

                                                                                   ccATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATgccgac           >  2:9017526/1‑91 (MQ=255)

                                                                                    cATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCa          <  1:1085046/93‑1 (MQ=255)

                                                                                    cATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCa          >  2:7856400/1‑93 (MQ=255)

                                                                                    cATCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCa          <  2:3188459/93‑1 (MQ=255)

                                                                                     aTCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAg         >  2:5332855/1‑93 (MQ=255)

                                                                                     aTCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAg         >  1:539508/1‑93 (MQ=255)

                                                                                     aTCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAg         <  1:2421198/93‑1 (MQ=255)

                                                                                      tCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGa        >  2:2678458/1‑93 (MQ=255)

                                                                                      tCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGa        >  2:8315909/1‑93 (MQ=255)

                                                                                      tCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGa        <  1:4353978/93‑1 (MQ=255)

                                                                                      tCCGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGa        >  1:3313752/1‑93 (MQ=255)

                                                                                       ccGCGTGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGAc       <  1:65094/93‑1 (MQ=255)

                                                                                            cgGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCa  <  2:347489/93‑1 (MQ=255)

CAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCA  >  NC_000913/4009627‑4009811

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`