BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,503,532	T→C	T43A (ACG→GCG)	yhfS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,503,532	0	T	C	85.7%	657.7 / 101.0	280	T43A (ACG→GCG)	yhfS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (21/19); new base C (93/147); total (114/166)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 285 total aligned reads displayed.

CCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAA  >  NC_000913/3503441‑3503624

 ccTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTa                                                                                            >  1:4013961/1‑93 (MQ=255)

  cTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTaa                                                                                           <  2:40965/93‑1 (MQ=255)

   tGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAg                                                                                          >  2:7348494/1‑93 (MQ=255)

      aCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCt                                                                                       >  2:5750700/1‑93 (MQ=255)

       ccAGCGCCGCAGCCAGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTa                                                                                      >  2:7489694/1‑93 (MQ=255)

        cAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTaa                                                                                     <  1:2396020/93‑1 (MQ=255)

         aGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTaaa                                                                                    <  2:2398218/93‑1 (MQ=255)

         aGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTaaa                                                                                    <  1:3133918/93‑1 (MQ=255)

                  gCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCa                                                                           >  1:1382811/1‑93 (MQ=255)

                     tgtgCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTa                                                                        >  2:3591161/1‑93 (MQ=255)

                     tgtgCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTa                                                                        >  2:6493827/1‑93 (MQ=255)

                      gtgCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTaa                                                                       <  1:5796710/93‑1 (MQ=255)

                      gtgCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTaa                                                                       >  1:464399/1‑93 (MQ=255)

                        gCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGa                                                                     <  1:6855058/93‑1 (MQ=255)

                         cGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGaa                                                                    <  1:7109378/93‑1 (MQ=255)

                           tGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAAc                                                                  <  2:4955023/93‑1 (MQ=255)

                            gAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACt                                                                 <  1:4719174/93‑1 (MQ=255)

                             aaaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTc                                                                <  2:2909683/93‑1 (MQ=255)

                             aaaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTc                                                                <  2:2799314/93‑1 (MQ=255)

                             aaaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTc                                                                >  1:5221551/1‑93 (MQ=255)

                             aaaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTc                                                                <  2:2599352/93‑1 (MQ=255)

                             aaaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTc                                                                <  2:4604931/93‑1 (MQ=255)

                              aaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCg                                                               <  1:4318658/93‑1 (MQ=255)

                               aTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGc                                                              <  1:3460102/93‑1 (MQ=255)

                               aTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCGGTAAGAAACTCGc                                                              >  2:3797218/1‑93 (MQ=255)

                               aTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGc                                                              >  1:4040668/1‑93 (MQ=255)

                               aTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGc                                                              <  1:779531/93‑1 (MQ=255)

                                tGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCt                                                             <  2:4425206/93‑1 (MQ=255)

                                  cATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGc                                                           <  2:5677389/93‑1 (MQ=255)

                                  cATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGc                                                           <  1:7240203/93‑1 (MQ=255)

                                   aTCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGcc                                                          >  1:2717362/1‑93 (MQ=255)

                                    tCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCg                                                         <  2:2686336/93‑1 (MQ=255)

                                    tCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCg                                                         >  2:2038402/1‑93 (MQ=255)

                                       gccagcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCgggg                                                      >  1:528024/1‑93 (MQ=255)

                                        ccagcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGa                                                     >  2:1355281/1‑93 (MQ=255)

                                          agcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaaa                                                   <  1:1972038/93‑1 (MQ=255)

                                          agcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaaa                                                   <  2:162521/93‑1 (MQ=255)

                                          agcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaaa                                                   <  2:37181/93‑1 (MQ=255)

                                           gcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAAt                                                  >  2:6075166/1‑93 (MQ=255)

                                           gcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAAt                                                  <  1:2997960/93‑1 (MQ=255)

                                           gcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAAt                                                  <  1:3866021/93‑1 (MQ=255)

                                            cACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATg                                                 <  1:7456772/93‑1 (MQ=255)

                                            cACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATg                                                 <  2:1815343/93‑1 (MQ=255)

                                            cACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATg                                                 <  2:5530452/93‑1 (MQ=255)

                                             aCCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATgg                                                <  1:561844/93‑1 (MQ=255)

                                              ccTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGc                                               >  1:1515440/1‑93 (MQ=255)

                                              ccTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGc                                               <  2:4553303/93‑1 (MQ=255)

                                                tGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGa                                             >  1:2107853/1‑93 (MQ=255)

                                                tGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGa                                             <  2:6686316/93‑1 (MQ=255)

                                                 gTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGAc                                            <  1:5511772/93‑1 (MQ=255)

                                                 gTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGAc                                            <  1:5022055/93‑1 (MQ=255)

                                                  ttCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACa                                           >  2:2283551/1‑93 (MQ=255)

                                                   tCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaa                                          <  1:5147376/93‑1 (MQ=255)

                                                    ccACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaaa                                         >  1:3660358/1‑93 (MQ=255)

                                                    ccACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaaa                                         >  2:3535222/1‑93 (MQ=255)

                                                     cacaCGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAAt                                        <  2:4692921/93‑1 (MQ=255)

                                                       cacGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGc                                      <  1:3862465/93‑1 (MQ=255)

                                                        acGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCt                                     <  2:7055554/93‑1 (MQ=255)

                                                        acGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCt                                     <  1:2814110/93‑1 (MQ=255)

                                                         cGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTa                                    <  1:6340291/93‑1 (MQ=255)

                                                           cTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATc                                  <  2:4023795/93‑1 (MQ=255)

                                                           cTGGGTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATc                                  <  1:101290/93‑1 (MQ=255)

                                                             gggTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATcca                                <  2:6214271/93‑1 (MQ=255)

                                                               gTAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccacc                              <  2:1926234/93‑1 (MQ=255)

                                                                tAACGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccacca                             >  2:7258465/1‑93 (MQ=255)

                                                                 aaCGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAg                            >  2:6399242/1‑93 (MQ=255)

                                                                 aaCGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAg                            <  2:2556390/93‑1 (MQ=255)

                                                                 aaCGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAg                            <  1:6781367/93‑1 (MQ=255)

                                                                  aCGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgc                           <  2:5891489/93‑1 (MQ=255)

                                                                  aCGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgc                           <  1:1073065/93‑1 (MQ=255)

                                                                  aCGCGCGGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgc                           <  2:5935361/93‑1 (MQ=255)

                                                                    gcgcgGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgcgc                         >  2:5128711/1‑93 (MQ=255)

                                                                      gcgGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCaa                       <  2:173729/93‑1 (MQ=255)

                                                                       cgGTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCaaa                      <  1:935856/93‑1 (MQ=255)

                                                                        ggTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAAc                     <  1:243195/93‑1 (MQ=255)

                                                                         gTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACt                    >  2:3133006/1‑93 (MQ=255)

                                                                         gTTGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACt                    >  1:2574260/1‑93 (MQ=255)

                                                                          ttGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTg                   <  2:5357114/93‑1 (MQ=255)

                                                                           tGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGc                  <  1:3706756/93‑1 (MQ=255)

                                                                           tGATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGc                  <  1:2133782/93‑1 (MQ=255)

                                                                            gATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCt                 <  2:4999101/93‑1 (MQ=255)

                                                                            gATTCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCt                 <  2:2255109/93‑1 (MQ=255)

                                                                             aTTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtt                <  2:7161894/93‑1 (MQ=255)

                                                                              ttCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCttt               <  2:4858817/93‑1 (MQ=255)

                                                                              ttCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCttt               >  1:169501/1‑93 (MQ=255)

                                                                               tCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtttt              <  2:2139554/93‑1 (MQ=255)

                                                                               tCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtttt              <  2:6167232/93‑1 (MQ=255)

                                                                               tCAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtttt              >  1:5705929/1‑93 (MQ=255)

                                                                                cAGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtg             <  1:1636279/93‑1 (MQ=255)

                                                                                 aGTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgc            <  1:1848653/93‑1 (MQ=255)

                                                                                  gTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgct           >  2:5188117/1‑93 (MQ=255)

                                                                                  gTCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgct           >  2:150324/1‑93 (MQ=255)

                                                                                   tCCTGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgctg          >  2:949790/1‑93 (MQ=255)

                                                                                    ccTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgctgc         <  2:506819/93‑1 (MQ=255)

                                                                                      tGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTgcgc       >  1:7301988/1‑93 (MQ=255)

                                                                                      tGGCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTgcgc       >  2:1212613/1‑93 (MQ=255)

                                                                                       ggCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGcc      >  2:688123/1‑93 (MQ=255)

                                                                                       ggCGCTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGcc      <  2:2973907/93‑1 (MQ=255)

                                                                                        gcgcTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCt     <  1:6069119/93‑1 (MQ=255)

                                                                                          gTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCa   >  2:3919708/1‑93 (MQ=255)

CCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAA  >  NC_000913/3503441‑3503624

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