BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,098	C→T	15.9%	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	15.9%	271.0 / 35.5	113	D40D (GAC→GAT)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base C (49/46); new base T (8/10); total (57/56)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 115 total aligned reads displayed.

CTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGG  >  NC_000913/1406000‑1406195

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

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

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

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

          aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                        <  1:617045/100‑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)

                tGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGAGCCGCAGCGTATCGACACGCTGCAGATTGATaa                                                                                  >  2:1129019/1‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                 aTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                 <  1:2449208/100‑1 (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)

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

                                                      cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTaa                                            >  1:1138325/1‑100 (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)

                                                             aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCACACTAAGTCAGcc                                     <  2:2544199/100‑1 (MQ=255)

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

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

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

                                                                  aCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTgg                                >  1:1791236/1‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                tGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                  >  2:1190889/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)

                                                                                    ctgcAGCGTGTCGATACGCTGGAGCTTGATAATTTCCACACCACCGGATTTCTCGACATCATCCCACTAATTCAGCCGCTGGAGTTTCGGCGGGAAgg                >  2:4033892/3‑98 (MQ=255)

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

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

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

                                                                                    ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggga              >  2:2018369/1‑99 (MQ=255)

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

                                                                                    ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt              <  2:2950745/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)

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

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

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

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

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

                                                                                               cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg   <  1:4032091/100‑1 (MQ=255)

                                                                                                gccacGCTGCAGCTTGATAATTTCCGCACCACGGGATGTCGCGACATCAGCCCACTAAGTCAGCCGGGGGAGTTTAGGAGGGAAGGGTTGCGAGATgggg  >  2:3554331/3‑100 (MQ=255)

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

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

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

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

CTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGG  >  NC_000913/1406000‑1406195

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