BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,098	C→T	17.3%	D40D (GAC→GAT)	smrA →	DNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,098	0	C	T	17.3%	408.8 / 65.2	173	D40D (GAC→GAT)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base C (81/61); new base T (19/11); total (101/72)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 181 total aligned reads displayed.

GCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTG  >  NC_000913/1405999‑1406197

gCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                     >  1:2179327/1‑100 (MQ=255)

 cTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                    <  2:285597/100‑1 (MQ=255)

   gTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACg                                                                                                  <  2:3915690/100‑1 (MQ=255)

       cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGc                                                                                              <  2:3808936/100‑1 (MQ=255)

       cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                              <  1:2689834/100‑1 (MQ=255)

       cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATATCCACTGGCATTCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                              >  2:2937403/1‑100 (MQ=255)

         tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                            <  1:2857292/100‑3 (MQ=255)

         tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                            <  2:3377206/100‑3 (MQ=255)

         tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                            <  1:1818799/100‑3 (MQ=255)

          gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                           <  1:3359836/100‑1 (MQ=255)

           ncGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCaaga                                                                                          >  1:1935728/2‑97 (MQ=255)

           aCGCCATGGAAGCTGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                          >  1:2762289/1‑100 (MQ=255)

           aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                          >  1:1941805/1‑100 (MQ=255)

            cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCtt                                                                                         <  1:671675/100‑1 (MQ=255)

             gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGcgca                                                                                        <  2:2881661/100‑4 (MQ=39)

             gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGcgca                                                                                        <  2:1644707/100‑4 (MQ=39)

              ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                       >  1:1038059/1‑100 (MQ=255)

              ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                       >  1:2949974/1‑100 (MQ=255)

                aTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATa                                                                                     <  1:1502876/100‑1 (MQ=255)

                     aGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTc                                                                                >  1:3575687/1‑100 (MQ=255)

                       aTGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCt                                                                              <  1:930982/100‑1 (MQ=255)

                        tGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCACAGCGTATCGACACGCTGCAGCTTGATAATTTCCTc                                                                             >  2:3602287/1‑100 (MQ=255)

                         gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                            >  1:2187126/1‑100 (MQ=255)

                         gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                            <  2:3289715/100‑1 (MQ=255)

                          tCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcac                                                                           <  2:1952519/100‑1 (MQ=255)

                                cGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCgg                                                                     <  1:563941/100‑1 (MQ=255)

                                 gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                    >  2:3046863/1‑100 (MQ=255)

                                 gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                    <  2:3663212/100‑1 (MQ=255)

                                 gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCCACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                    >  2:1451235/1‑100 (MQ=255)

                                           gCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAca                                                          >  1:1469610/1‑100 (MQ=255)

                                            cTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcat                                                         >  1:911692/1‑100 (MQ=255)

                                            cTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcat                                                         >  1:1058265/1‑100 (MQ=255)

                                             tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatc                                                        <  1:1112588/100‑1 (MQ=255)

                                               ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatcat                                                      >  1:3154835/1‑100 (MQ=255)

                                                cGATGTCCACTGGCATCCAACGCGTAACCAATGTGCGACGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACatcatc                                                     <  2:2824924/100‑1 (MQ=255)

                                                 gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcc                                                    >  2:2170179/1‑100 (MQ=255)

                                                 gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcc                                                    <  2:2454982/100‑1 (MQ=255)

                                                  aTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                   <  1:638606/100‑1 (MQ=255)

                                                   tGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATACTTTCCTCACCACCGGATTTCTCGACATCATCCCa                                                  <  1:839868/100‑1 (MQ=255)

                                                   tGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCa                                                  >  2:1025993/1‑100 (MQ=255)

                                                    gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                 <  2:3874431/100‑1 (MQ=255)

                                                    gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                 >  1:3187861/1‑100 (MQ=255)

                                                    gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                 <  1:2343132/100‑1 (MQ=255)

                                                    gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                 <  1:124987/100‑1 (MQ=255)

                                                     tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                                <  1:3955052/100‑1 (MQ=255)

                                                     tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                                >  1:3048277/1‑100 (MQ=255)

                                                     tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                                >  2:82846/1‑100 (MQ=255)

                                                      ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACt                                                >  2:542207/1‑99 (MQ=255)

                                                        aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTaa                                              >  1:4091068/1‑99 (MQ=255)

                                                        aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAg                                             <  2:1424304/100‑1 (MQ=255)

                                                           ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                          >  1:1860951/1‑100 (MQ=255)

                                                           ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                          >  1:2505890/1‑100 (MQ=255)

                                                           ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAGTTCGCGACATCATCCCACTAAGTCa                                          >  2:3992157/1‑100 (MQ=255)

                                                             cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACAGGCTGCAGCTTGATAATTTCCTCACCACCGGATTTGTCGACATCATCCCACTAAGTCAGc                                        >  1:382133/1‑100 (MQ=255)

                                                             cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAg                                         >  2:1685930/1‑99 (MQ=255)

                                                                ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                     >  2:2096006/1‑100 (MQ=255)

                                                                ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                     >  2:1760649/1‑100 (MQ=255)

                                                                     gcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                                <  2:2684494/100‑1 (MQ=255)

                                                                     gcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGATTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                                <  2:1466789/100‑1 (MQ=255)

                                                                      cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGt                               >  1:373343/1‑100 (MQ=255)

                                                                       gTAACCAACGTTCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGtt                              <  2:3592916/100‑1 (MQ=255)

                                                                       gTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCTCCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAg                                <  2:2517630/98‑1 (MQ=255)

                                                                       gTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGtt                              >  2:348337/1‑100 (MQ=255)

                                                                        tAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGttt                             >  1:2627208/1‑100 (MQ=255)

                                                                        tAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGttt                             <  2:3048264/100‑1 (MQ=255)

                                                                           ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                          >  2:87743/1‑100 (MQ=255)

                                                                           ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                          >  2:2311799/1‑100 (MQ=255)

                                                                            cAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggc                         <  2:2966379/100‑1 (MQ=255)

                                                                             aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggcg                        <  2:1126508/100‑1 (MQ=255)

                                                                              aCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggcgg                       <  1:2309658/100‑1 (MQ=255)

                                                                               cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCggg                      <  1:3318066/100‑1 (MQ=255)

                                                                               atcaagctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAACTCAGCCGCTGGAGTTTcggc                         >  2:117208/9‑97 (MQ=38)

                                                                                 caagctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                    >  1:984084/7‑100 (MQ=38)

                                                                                 caagctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTCTCGGCGGGaa                    >  1:3216797/7‑100 (MQ=37)

                                                                                   atctgcAGCGTGTAGATACGCTGCAGCTTGAAAATTTCCTCACCACCGGATTTCTCGACATCATCACACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                  >  1:2822802/5‑100 (MQ=255)

                                                                                   agctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                  >  2:3475478/5‑100 (MQ=255)

                                                                                   agctgcAGCGTGTCGATACGCGGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGCAAgg                  >  2:2063890/5‑100 (MQ=255)

                                                                                    gctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg                 >  1:1184729/4‑100 (MQ=38)

                                                                                    gctgcAGCGTGTCGATACGCGGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTACGTCAGCCGCTGGAGTTTcggcgg                       >  2:538433/4‑94 (MQ=255)

                                                                                     ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                >  1:3737414/3‑100 (MQ=255)

                                                                                     ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                >  1:3686945/1‑100 (MQ=255)

                                                                                      tgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt               >  1:3603562/2‑100 (MQ=255)

                                                                                      tgcAGCGTGTCGATACGCTGCAGCGTGATAATTTCCTCACCCCCGGCTTTCTCGACACCATCCCACTAAGTCAGCCGCTCGCGTTTCGGCAGGa                     >  2:2422260/2‑94 (MQ=255)

                                                                                          gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCaa           >  1:3636602/1‑100 (MQ=255)

                                                                                           cGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAc          >  1:3250409/1‑100 (MQ=255)

                                                                                            gTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACa         <  1:2495990/100‑1 (MQ=255)

                                                                                            gTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACa         <  2:3092930/100‑1 (MQ=255)

                                                                                             tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACAt        >  2:4150104/1‑100 (MQ=255)

                                                                                             tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTAAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACAt        >  2:1864838/1‑100 (MQ=255)

                                                                                               tCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgg      <  1:1025993/100‑1 (MQ=255)

                                                                                                cGACACGCTGCAGCTTGATACTTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg     <  2:3152172/100‑1 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg    >  1:264221/1‑100 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg    <  2:2991152/100‑1 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg    >  1:2130142/1‑100 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg    <  2:575444/100‑1 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg    >  2:802202/1‑100 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg    >  1:1690350/1‑100 (MQ=255)

                                                                                                 gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg     <  1:757400/99‑1 (MQ=255)

                                                                                                   cacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTg  >  1:3886801/1‑100 (MQ=255)

                                                                                                   cacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTg  >  1:1351520/1‑100 (MQ=255)

GCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTG  >  NC_000913/1405999‑1406197

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