BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,098	C→T	18.1%	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	18.1%	230.9 / 40.0	105	D40D (GAC→GAT)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base C (37/49); new base T (5/14); total (42/63)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-02

Only 100 of 110 total aligned reads displayed.

CTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGT  >  NC_000913/1406000‑1406196

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

 tGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAAAGTGCGCCACAGCGTATCGacac                                                                                                  >  1:512377/1‑100 (MQ=255)

   tttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGc                                                                                                <  1:2499280/100‑1 (MQ=255)

      cTTGACGCTATGGAAGATGGGCAGCCGCTGAAACGTGCTACCGGTGTGCATTGTCATCCAATTTTTTTCTATTGTTCGTCGTAGCGTATCGACACGCTgg                                                                                             <  2:614288/100‑2 (MQ=255)

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

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

       ttGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCATAACCAACGTGCGCCGCAGCGTATCGACACGCTGCa                                                                                            >  1:1422857/1‑100 (MQ=255)

        tGACGCCATGGAAGATGTCCAGCTGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                           <  1:617911/100‑3 (MQ=255)

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

         gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTATCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                          <  1:2814543/100‑1 (MQ=255)

         gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                          <  1:1704726/100‑1 (MQ=255)

         gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                          <  1:2814730/100‑1 (MQ=255)

         gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                          <  2:2703934/100‑1 (MQ=255)

         gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                          <  1:1974206/100‑1 (MQ=255)

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

         gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                          >  2:593100/1‑100 (MQ=255)

          aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGcg                                                                                         <  1:2696222/100‑2 (MQ=255)

          aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                         <  1:2691615/100‑1 (MQ=255)

          aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                         <  1:2551525/100‑1 (MQ=255)

           cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCCACGTGCGCCGCAGCGTATCGACACGCTGCAGCtt                                                                                        >  1:2213178/1‑100 (MQ=255)

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

            gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcggcgca                                                                                       <  2:31383/100‑7 (MQ=39)

            gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcggcgca                                                                                       <  1:2545421/100‑7 (MQ=39)

             ccATGGAATATGTCCAGCCGCTGAAACGTCCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGCCGATACGCTGCGGc                                                                                          <  2:2264679/96‑1 (MQ=255)

             ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcggcgcac                                                                                      <  1:1333197/100‑8 (MQ=38)

              cATGGAAGATGTCCAGCCGCTAAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGAt                                                                                     <  2:1941220/100‑1 (MQ=255)

               aTGGAAGATGTCCAGCCGCTGAAACGTGTTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGAt                                                                                     <  1:2916457/99‑1 (MQ=255)

                 ggAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                                  >  1:2646927/1‑100 (MQ=255)

                   aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCATCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAttt                                                                                <  2:1646675/100‑1 (MQ=255)

                   aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAttt                                                                                <  1:1216581/100‑1 (MQ=255)

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

                     gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                              <  1:2223614/100‑1 (MQ=255)

                     gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                              <  1:2271584/100‑1 (MQ=255)

                     gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                              <  1:465943/100‑1 (MQ=255)

                     gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                              <  1:615081/100‑1 (MQ=255)

                     gATGTCCAGCCGCTGAAACGTGCTAACGATGTCAACTGGCATACAACGCGTAACCAACGTGCGCCGCAGCGTATCGACAGGCTGCAGCTGGATAATTTc                                                                               >  2:1273284/1‑99 (MQ=255)

                         tgcAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcac                                                                          >  1:2059177/3‑100 (MQ=255)

                          ccAGCTGCTGAAACGTGCTACCGATGTCCACTTGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcacc                                                                         <  1:1359909/100‑1 (MQ=255)

                           cAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcacca                                                                        <  1:61325/100‑1 (MQ=255)

                             gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCaccacc                                                                      <  1:767249/100‑1 (MQ=255)

                                 cTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAt                                                                  <  1:593103/100‑1 (MQ=255)

                                       cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                            >  2:1421330/1‑100 (MQ=255)

                                        gTGCTACCGATGTCCACTGGCATCCAACGTGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGa                                                           <  1:848542/100‑1 (MQ=255)

                                         tGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAc                                                          >  2:1398661/1‑100 (MQ=255)

                                          gCTACCGATGTCCACTGGCATCCAACGCGTAATCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGa                                                           <  2:2019887/98‑1 (MQ=255)

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

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

                                               cGATGTCCACTGGCATCCAACGCGTAACCATCGTGCCCCGCAGCGTCTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatcat                                                     <  2:2051031/99‑1 (MQ=255)

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

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

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

                                                 aTGTCCACTGGCATCCAACGTGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                  <  1:1288139/100‑1 (MQ=255)

                                                 aTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGTGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                  <  1:2689698/100‑1 (MQ=255)

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

                                                     ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                              >  2:1783921/1‑100 (MQ=255)

                                                     ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                              >  1:334956/1‑100 (MQ=255)

                                                     ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                              <  1:1889811/100‑1 (MQ=255)

                                                        cTGGCATCCAACGCGTAACCAACGTGCGCTGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                           <  2:242978/100‑1 (MQ=255)

                                                        cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                           <  1:1201444/100‑1 (MQ=255)

                                                        cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                           <  1:2018712/100‑1 (MQ=255)

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

                                                            cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGc                                       >  1:1145582/1‑100 (MQ=255)

                                                            cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGc                                       >  1:1555307/1‑100 (MQ=255)

                                                             aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                      >  2:1774682/1‑100 (MQ=255)

                                                             aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                      <  2:459260/100‑1 (MQ=255)

                                                              tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAg                                        <  2:3025011/97‑1 (MQ=255)

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

                                                                cAACGCGTAACCAACGTGCGCCGCAGCGAATCCACACGCTGCAGCTTGATACTTTCCTCACCACCCGATTTCTCGACATCATCGCACTAAGTCAGCCGc                                    >  2:1343730/1‑99 (MQ=255)

                                                                 aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                  <  2:2421018/100‑1 (MQ=255)

                                                                 aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                  <  1:1735141/100‑1 (MQ=255)

                                                                     cgTAACCAACGTGCGCTACAGCATATCGACACGATGAAGCTTGATAATTTCCTCACTAACAGATTTATCGACATTATCCAACTATGTCAGTCGCTTGAGt                              >  2:1436833/1‑100 (MQ=255)

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

                                                                         aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcg                          <  1:2309853/100‑1 (MQ=255)

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

                                                                          aattatcaagccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                         <  1:2654741/91‑1 (MQ=255)

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

                                                                                tGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                   >  2:2103048/1‑100 (MQ=255)

                                                                                  cgccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                 <  2:1957145/100‑1 (MQ=255)

                                                                                   gctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg                >  1:377187/4‑100 (MQ=255)

                                                                                   gctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg                >  1:50968/4‑100 (MQ=255)

                                                                                   gctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCAACATCATACAACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg                >  1:820094/4‑100 (MQ=255)

                                                                                    ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt               >  2:2545096/3‑100 (MQ=255)

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

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

                                                                                     cgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt              >  2:1719354/1‑100 (MQ=255)

                                                                                       cAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTg             <  2:2451566/99‑1 (MQ=255)

                                                                                        aGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCa           <  2:654245/100‑1 (MQ=255)

                                                                                        aGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCa           <  1:1967161/100‑1 (MQ=255)

                                                                                          cGTATCGACACGCTGCAGCTTGATAATTTCCTCACCGCCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGATTTTCGGTGGGAAGGGTTGCa           <  2:1879932/98‑1 (MQ=255)

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

                                                                                             aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg      <  1:1731375/100‑1 (MQ=255)

                                                                                              tcaacacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgg     <  1:2917528/97‑1 (MQ=255)

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

                                                                                               cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg    >  2:2969683/1‑100 (MQ=255)

                                                                                               cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg    >  2:420084/1‑100 (MQ=255)

                                                                                                gACACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg   <  2:1198202/100‑1 (MQ=255)

                                                                                                gACACGCTGCAGCTTGATAATTTCCTCACCACCGTTTTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg   <  1:1804693/100‑1 (MQ=255)

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

                                                                                                 acacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt  >  1:462726/1‑100 (MQ=255)

                                                                                                 acacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAATCAACCGCTGGAGTTTAAGCGGAAAGGGTTGCAAAATGGGGt  >  1:463053/1‑100 (MQ=255)

CTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGT  >  NC_000913/1406000‑1406196

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