BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,239,132	T→C	H186R (CAC→CGC)	galS ←	mgl repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,239,132	0	T	C	99.7%	1162.0 / ‑5.3	337	H186R (CAC→CGC)	galS	mgl repressor
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (181/155); minor base A (1/0); total (182/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 350 total aligned reads displayed.

GCGCCAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGTGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAGGCAA  >  NC_000913/2239042‑2239221

 cgcCAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGc                                                                                        <  2:427677/93‑1 (MQ=255)

   ccAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTg                                                                                      >  2:8030122/1‑93 (MQ=255)

    cAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTgg                                                                                     >  1:5507666/1‑93 (MQ=255)

    cAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTgg                                                                                     >  1:1370864/1‑93 (MQ=255)

    cAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTgg                                                                                     >  2:733182/1‑93 (MQ=255)

    cAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTgg                                                                                     >  1:699580/1‑93 (MQ=255)

      aTCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaa                                                                                   >  2:5862114/1‑93 (MQ=255)

       tCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaag                                                                                  >  1:2085310/1‑93 (MQ=255)

       tCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaag                                                                                  >  1:4467414/1‑93 (MQ=255)

          aGCTTTCCGGCTGAATAATATCCTGCTCTTTTAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaa                                                                               <  2:7957505/93‑1 (MQ=255)

           gCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaag                                                                              <  2:8626200/93‑1 (MQ=255)

           gCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaag                                                                              <  1:6188284/93‑1 (MQ=255)

           gCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaag                                                                              >  2:5535362/1‑93 (MQ=255)

             tttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGAt                                                                            >  2:5688826/1‑93 (MQ=255)

             tttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGAt                                                                            >  2:5964018/1‑93 (MQ=255)

             tttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGAGTCATCTTCAATGCCGCGGCTGGAAGAAAGAt                                                                            >  2:6910277/1‑93 (MQ=255)

             cttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGAt                                                                            <  1:2923906/92‑1 (MQ=255)

              ttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGCCATCTTCAATTCCGCGGCTGGAAGAAACATa                                                                           >  1:1184386/1‑93 (MQ=255)

                ccGGCGGAATAATATCCTGCTCTTTCAACGCACTCAGCCAGCATGCGTTACGCAGGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAAc                                                                         <  1:75576/93‑1 (MQ=255)

                ccGGCGGAATAATATCCTGCTCTTTCAAAGCACTCATCCAGCCTGCTTTACGCATCGCGGCCTATTCAATGGCGCGGCTGGCAGAAAGCTAAc                                                                         >  1:3445989/1‑93 (MQ=255)

                   gcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCaa                                                                      <  1:3663401/93‑1 (MQ=255)

                   gcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCaa                                                                      <  2:1706475/93‑1 (MQ=255)

                   gcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCaa                                                                      >  1:6450732/1‑93 (MQ=255)

                    cggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAAt                                                                     <  1:774895/93‑1 (MQ=255)

                    cggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAAt                                                                     >  2:5889569/1‑93 (MQ=255)

                    cggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATCGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAAt                                                                     >  2:8695380/1‑93 (MQ=255)

                     ggAATAATGTCCTGCTCTTTCAACGCGCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATa                                                                    <  2:1762185/93‑1 (MQ=255)

                     ggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATa                                                                    >  1:5211973/1‑93 (MQ=255)

                     ggAATAATATCCTGCTCTTTAAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATa                                                                    <  1:7098632/93‑1 (MQ=255)

                      gAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATAc                                                                   <  1:4641720/93‑1 (MQ=255)

                       aataatATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACg                                                                  <  1:3821742/93‑1 (MQ=255)

                         taataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGtt                                                                <  1:8528904/93‑1 (MQ=255)

                         taataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGtt                                                                >  1:8217640/1‑93 (MQ=255)

                         taataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGtt                                                                <  1:7178611/93‑1 (MQ=255)

                           atatCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtga                                                              <  1:3629533/93‑1 (MQ=255)

                           atatCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtga                                                              <  2:2016862/93‑1 (MQ=255)

                            tatCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATAAGTtgat                                                             >  2:8477272/1‑93 (MQ=255)

                            aatCCTGCTCCTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgat                                                             <  1:7653729/92‑1 (MQ=255)

                             atCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgatg                                                            <  1:1534581/93‑1 (MQ=255)

                             atCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATAAGTtgatg                                                            >  1:6844949/1‑93 (MQ=255)

                              tCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgatga                                                           >  1:3350858/1‑93 (MQ=255)

                              tCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgatga                                                           <  2:8153044/93‑1 (MQ=255)

                                cTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGAcc                                                         <  2:3432155/93‑1 (MQ=255)

                                 tGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCg                                                        >  2:4716411/1‑93 (MQ=255)

                                    tctTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTa                                                     >  2:2546179/1‑93 (MQ=255)

                                    tctTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTa                                                     >  1:5901741/1‑93 (MQ=255)

                                       ttCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTc                                                  <  2:2201323/93‑1 (MQ=255)

                                       ttCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTc                                                  >  2:7724247/1‑93 (MQ=255)

                                        tCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCa                                                 <  1:1919583/93‑1 (MQ=255)

                                         cAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAg                                                <  2:4066775/93‑1 (MQ=255)

                                          aaCGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGc                                               >  1:1000650/1‑93 (MQ=255)

                                           aCGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCa                                              <  2:6568810/93‑1 (MQ=255)

                                              cACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACa                                           <  1:7362635/93‑1 (MQ=255)

                                               aCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACAt                                          >  2:3611700/1‑93 (MQ=255)

                                               aCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACAt                                          <  2:2155226/93‑1 (MQ=255)

                                               aCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACAt                                          >  2:2372230/1‑93 (MQ=255)

                                                cTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATg                                         <  2:6169641/93‑1 (MQ=255)

                                                cTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATg                                         <  1:8073465/93‑1 (MQ=255)

                                                cTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATg                                         <  2:2337708/93‑1 (MQ=255)

                                                 tCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATgc                                        <  2:1417848/93‑1 (MQ=255)

                                                  cATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATgcg                                       <  1:2670428/93‑1 (MQ=255)

                                                  cATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATgcg                                       >  2:1920837/1‑93 (MQ=255)

                                                   aTCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATgcgc                                      >  2:3199863/1‑93 (MQ=255)

                                                     ccAACCTTCTTTTCGGCTGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGt                                    <  1:2295584/93‑1 (MQ=255)

                                                          cTGCTTTACGCATGGCGTCATATTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCa                               >  2:7475135/1‑93 (MQ=255)

                                                           tGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAt                              >  2:8259512/1‑93 (MQ=255)

                                                           tGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAt                              >  2:6413623/1‑93 (MQ=255)

                                                           tGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAt                              >  2:2875568/1‑93 (MQ=255)

                                                           tGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAt                              >  1:6760278/1‑93 (MQ=255)

                                                           tGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAt                              >  2:1542229/1‑93 (MQ=255)

                                                              tttACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCg                           <  1:1260696/93‑1 (MQ=255)

                                                               ttACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCgg                          <  1:6615895/93‑1 (MQ=255)

                                                               ttACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCgg                          >  1:2784382/1‑93 (MQ=255)

                                                                 aCGCAGGGCGTCATCTTCAATGCAGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGc                        <  1:6335143/93‑1 (MQ=255)

                                                                  cGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCa                       >  1:1500954/1‑93 (MQ=255)

                                                                   gCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCAc                      >  2:4838356/1‑93 (MQ=255)

                                                                   gCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCAc                      <  2:6261333/93‑1 (MQ=255)

                                                                   gCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCAc                      <  1:981489/93‑1 (MQ=255)

                                                                     aTGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCg                    >  2:7920038/1‑93 (MQ=255)

                                                                      tGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGc                   >  2:6717510/1‑93 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                 >  2:3744786/1‑93 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                 >  1:8337130/1‑93 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                 <  2:1225239/93‑1 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                 <  2:4983542/93‑1 (MQ=255)

                                                                         cGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTga                <  1:3972810/93‑1 (MQ=255)

                                                                             aTCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGAtt            >  2:6325949/1‑93 (MQ=255)

                                                                             aTCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGAtt            <  1:4961625/93‑1 (MQ=255)

                                                                              tCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTa           <  2:2597784/93‑1 (MQ=255)

                                                                              tCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTa           <  1:4771825/93‑1 (MQ=255)

                                                                                ttCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATc         >  1:4326340/1‑93 (MQ=255)

                                                                                ttCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATc         >  2:3971450/1‑93 (MQ=255)

                                                                                  cAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCa       <  1:3131033/93‑1 (MQ=255)

                                                                                  cAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCa       <  1:5609203/93‑1 (MQ=255)

                                                                                   aaTGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCATCGCTGAGATTATcctg      >  1:7820144/1‑91 (MQ=255)

                                                                                   aaTGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAg      <  1:2445789/93‑1 (MQ=255)

                                                                                   aaTGCCGCGGCTGAAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAg      <  2:5233007/93‑1 (MQ=255)

                                                                                    aTGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAgg     >  1:406687/1‑93 (MQ=255)

                                                                                    aTGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAgg     >  2:4526521/1‑93 (MQ=255)

                                                                                     tGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAGGc    >  2:7361287/1‑93 (MQ=255)

                                                                                       ccGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGCCGCCATTCGGGCACAGCTGAGATTATCCAGGCaa  >  2:279534/1‑93 (MQ=255)

GCGCCAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGTGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAGGCAA  >  NC_000913/2239042‑2239221

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