BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,093	A→G	18.4%	I39V (ATC→GTC)	smrA →	DNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,093	0	A	G	18.4%	207.5 / 44.1	104	I39V (ATC→GTC)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base A (39/45); new base G (5/14); total (44/60)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-02

Only 100 of 107 total aligned reads displayed.

AAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATG  >  NC_000913/1405994‑1406192

aaaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                     <  2:2621147/100‑1 (MQ=255)

aaaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                     <  2:2618219/100‑1 (MQ=255)

 aaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTAt                                                                                                    >  1:2262906/1‑100 (MQ=255)

 aaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTAt                                                                                                    >  1:2855708/1‑100 (MQ=255)

 aaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGACAACTGGAATCCAACGCGTAACCAACGTGCGCCGCAGCGTAt                                                                                                    >  1:1001479/1‑100 (MQ=255)

   tCGCTGTTTCTTGACGCCATGGAAGATGTTCAGCCGTTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCg                                                                                                  <  1:2383141/100‑1 (MQ=255)

   tCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCg                                                                                                  <  1:646967/100‑1 (MQ=255)

   tCGCTGTTTCTCGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCg                                                                                                  <  2:1887930/100‑1 (MQ=255)

      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:325548/100‑1 (MQ=255)

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

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

              tGACGCCATGGAAGATGTCCAGTCGTTGAAACGTGCTATCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  1:2475656/100‑3 (MQ=255)

              tGACGCCATGGAAGATGTCCAGTCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  1:1486622/100‑3 (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:1974206/100‑1 (MQ=255)

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

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

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

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

               gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                      <  2:532878/100‑1 (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)

                     aTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATa                                                                                >  2:148548/1‑100 (MQ=255)

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

                        gAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAtt                                                                             >  1:2002190/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:615081/100‑1 (MQ=255)

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

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

                           gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                          <  1:2271584/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)

                                   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:1433528/1‑100 (MQ=255)

                                                      gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATcc                                               >  1:2451595/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                                          <  1:1889811/100‑1 (MQ=255)

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

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

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

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

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

                                                                ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                     >  1:383820/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)

                                                                       aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                              >  2:2230324/1‑100 (MQ=255)

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

                                                                             tAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGttt                        >  1:1917041/1‑100 (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           >  1:3013046/3‑100 (MQ=255)

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

                                                                                          ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt           >  1:731879/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)

AAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATG  >  NC_000913/1405994‑1406192

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