BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,093	A→G	16.5%	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	16.5%	236.4 / 39.6	109	I39V (ATC→GTC)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base A (46/45); new base G (8/10); total (54/55)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 111 total aligned reads displayed.

AAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATG  >  NC_000913/1405994‑1406192

aaaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                     >  1:2812275/1‑100 (MQ=255)

  aTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATc                                                                                                   >  1:4081711/1‑100 (MQ=255)

  aTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATc                                                                                                   <  2:2154462/100‑1 (MQ=255)

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

      cTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                <  2:3923129/99‑1 (MQ=255)

        gTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACg                                                                                             >  2:818752/1‑100 (MQ=255)

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

          ttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCt                                                                                           >  2:108883/1‑100 (MQ=255)

           tCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTg                                                                                          >  2:1299120/1‑100 (MQ=255)

            cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                         >  2:3047614/1‑100 (MQ=255)

             ttGACGCCATGGAAGATGTCCAGCCGCTGAACCGTGCTACCGATGTCCAGTGGCATCCAACGCGTAACCAAGGTGCGCCGCAGCGTATCGACACGCTg                                                                                          >  2:3553381/1‑98 (MQ=255)

             ttGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCa                                                                                        <  1:2660770/100‑1 (MQ=255)

              tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACGAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                       <  1:3918588/100‑3 (MQ=255)

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

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

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

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

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

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

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

                 cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                      <  2:1184694/98‑1 (MQ=255)

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

                  gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCTGCAGCGTGTCGATACGCTGcggcgca                                                                                   <  1:1741864/100‑7 (MQ=38)

                   ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTgn                                                                                  <  2:3372177/100‑2 (MQ=255)

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

                    cATGGAAGATGTCCAGCCGCTGAAACGGGCTACCAATGTGCACTGGCATCCGAAGCGTAACAAACGAGCGACGCAGCGTATCGACACGCTGCAGATTg                                                                                   >  2:3442643/1‑98 (MQ=255)

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

                       ggAAGCTGTCCAGCCGCTGAAACGTGCTCCCGATGTCCACTGGCATCCAACGCGTCACCAACGTGAGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                              <  2:398463/100‑1 (MQ=255)

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

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

                         aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAttt                                                                            >  1:1896235/1‑100 (MQ=255)

                           gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                          >  2:2824132/1‑100 (MQ=255)

                             tGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTc                                                                        <  2:2501026/100‑1 (MQ=255)

                             tGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTc                                                                        >  1:4032298/1‑100 (MQ=255)

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

                                ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTcacc                                                                     <  1:3088761/100‑1 (MQ=255)

                                   gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCaccacc                                                                  >  2:697325/1‑100 (MQ=255)

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

                                       cTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAt                                                              >  2:2835221/1‑100 (MQ=255)

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

                                           aaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtct                                                          <  1:7105/100‑1 (MQ=255)

                                             cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                        <  2:574667/100‑1 (MQ=255)

                                             cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                        <  2:3567848/100‑1 (MQ=255)

                                              gTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGa                                                       <  2:1442326/100‑1 (MQ=255)

                                               tGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAc                                                      <  1:3582125/100‑1 (MQ=255)

                                                  tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatc                                                   >  2:4056364/1‑100 (MQ=255)

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

                                                   aCCGATGTCCACTGGCATCCAACGCGTAACCAACGTTCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                        <  2:2037590/94‑1 (MQ=255)

                                                    ccGCTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGGATCGACCCGCTGCAGCTTGCTAATTTCCTCGCCACCGGATTTCTCGAcatcat                                                 >  2:949716/1‑100 (MQ=255)

                                                    ccGATGTCCACTGGCATCCAACGCGACACCCCCCTGCGCAGCCGCGGATACACACACTGAAGCTTGATAATTTCCTCACCACCGGAGTtctc                                                         >  2:949703/1‑92 (MQ=255)

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

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

                                                          tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCACACt                                           >  2:3875786/1‑100 (MQ=255)

                                                          tCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCCACACGCTGAAGCTGGATAATTTCCTCACAAGCGGAATTCTCGACAGCATCCAAc                                            >  2:3726913/1‑99 (MQ=255)

                                                           ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGTTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                          <  2:3562622/100‑1 (MQ=255)

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

                                                            cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTaa                                         >  1:1138325/1‑100 (MQ=255)

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

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

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

                                                              cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                       >  2:3143886/1‑100 (MQ=255)

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

                                                              cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                       >  2:3143289/1‑100 (MQ=255)

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

                                                                     ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                <  2:3289043/100‑1 (MQ=255)

                                                                      cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCt                               <  1:2625267/100‑1 (MQ=255)

                                                                      cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCt                               >  1:1382779/1‑100 (MQ=255)

                                                                         cgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGa                            <  2:3767524/100‑1 (MQ=255)

                                                                         cgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGa                            >  1:3586658/1‑100 (MQ=255)

                                                                               aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcg                      >  1:1404072/1‑100 (MQ=255)

                                                                                ccAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAGGTCAGCCGCTGGAGTTTcgg                     <  2:290796/100‑1 (MQ=255)

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

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

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

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

                                                                                  aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcggc                    >  2:1582363/1‑99 (MQ=255)

                                                                                  aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGGCAGCCGCTGGAGTTTcggcg                   <  2:2523974/100‑1 (MQ=255)

                                                                                     tcaagctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCggg                 >  1:3814193/8‑99 (MQ=38)

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

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

                                                                                      tGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAGTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTgcggcg                   <  2:594161/96‑1 (MQ=255)

                                                                                         tccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg            <  1:2018382/99‑1 (MQ=255)

                                                                                         gctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg            >  2:2950791/4‑100 (MQ=255)

                                                                                         gccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg            >  1:3689618/1‑100 (MQ=255)

                                                                                         gatgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCCATAAGCCAACCGCTGAAGCTTCGGCGGGAAgg             >  2:3661825/4‑99 (MQ=255)

                                                                                          ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg            >  2:1002879/3‑99 (MQ=255)

                                                                                          ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg            >  2:3566309/3‑99 (MQ=255)

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

                                                                                          ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt           <  2:2950745/100‑1 (MQ=255)

                                                                                          ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt           <  2:1321530/100‑1 (MQ=255)

                                                                                          ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAACGGt           >  2:3668817/1‑100 (MQ=255)

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

                                                                                            gcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACAGCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt           <  2:1604994/98‑1 (MQ=255)

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

                                                                                               gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTCAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCaa      >  2:418595/1‑100 (MQ=255)

                                                                                               gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTTCaa      >  2:416571/1‑100 (MQ=255)

                                                                                               gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCaa      <  2:2779132/100‑1 (MQ=255)

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

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

                                                                                                   aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg  <  1:591235/100‑1 (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` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`