BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,665,031	A→G	V279A (GTG→GCG)	hcaT ←	HcaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,665,031	0	A	G	98.3%	366.7 / ‑3.1	119	V279A (GTG→GCG)	hcaT	HcaT
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (56/60); total (57/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 129 total aligned reads displayed.

AGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGC  >  NC_000913/2664941‑2665122

aGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCt                                                                                           >  1:2416303/1‑93 (MQ=255)

aGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCt                                                                                           >  1:4158641/1‑93 (MQ=255)

aGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCt                                                                                           >  1:2271945/1‑93 (MQ=255)

 ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTACCAACCACGGCAACGCCGGCGTTGCTCCCATAATGCCCCAGCGCGCTa                                                                                          >  1:1294230/1‑93 (MQ=255)

  gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAc                                                                                         >  1:4695773/1‑93 (MQ=255)

        gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCa                                                                                   >  2:2065512/1‑93 (MQ=255)

        gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCa                                                                                   >  1:19469/1‑93 (MQ=255)

        gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCa                                                                                   <  1:9881/93‑1 (MQ=255)

        gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCa                                                                                   >  2:4389860/1‑93 (MQ=255)

         aCCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCaa                                                                                  >  2:3905355/1‑93 (MQ=255)

          ccGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCaaa                                                                                 >  2:3511029/1‑93 (MQ=255)

           cGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAAt                                                                                >  2:4330234/1‑93 (MQ=255)

            gTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATc                                                                               >  2:122049/1‑93 (MQ=255)

            gTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATc                                                                               <  1:1254924/93‑1 (MQ=255)

            gTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATc                                                                               >  2:870894/1‑93 (MQ=255)

            gTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATc                                                                               <  1:3443779/93‑1 (MQ=255)

             tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCg                                                                              >  1:557316/1‑93 (MQ=255)

             tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCg                                                                              >  1:890989/1‑93 (MQ=255)

               aaGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGcc                                                                            >  2:3634068/1‑93 (MQ=255)

               aaGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCACGCTACGCCGCAAATCGcc                                                                            >  1:1410708/1‑93 (MQ=255)

               aaGGCGCCGCACTGCAGAATTTGCAACACTATCAACCACGGCAACGCCGAAGTTGCTCCCATAATGCCCCAGCGAGCTACGCCTCAAATCGcc                                                                            >  2:262500/1‑93 (MQ=255)

                aGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCg                                                                           >  2:36473/1‑93 (MQ=255)

                     ccGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATc                                                                      >  1:3597536/1‑93 (MQ=255)

                      cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCa                                                                     <  2:1911946/93‑1 (MQ=255)

                        cAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAAc                                                                   >  2:3555771/1‑93 (MQ=255)

                        cAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAAc                                                                   >  2:2205535/1‑93 (MQ=255)

                          aTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAg                                                                 >  2:2084080/1‑93 (MQ=255)

                            gCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCa                                                               <  1:297936/93‑1 (MQ=255)

                             cAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCat                                                              <  2:4581753/93‑1 (MQ=255)

                              aGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCata                                                             >  2:3556729/1‑93 (MQ=255)

                               gAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCatat                                                            >  1:34132/1‑93 (MQ=255)

                               gAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCatat                                                            >  2:640730/1‑93 (MQ=255)

                               gAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCatat                                                            <  1:4191255/93‑1 (MQ=255)

                                 aTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcg                                                          >  1:4650382/1‑93 (MQ=255)

                                  tttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTATGCCGCAAATCGCCGAGATCAACAGCATATcgc                                                         >  1:3078914/1‑93 (MQ=255)

                                  tttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcgc                                                         <  1:4082376/93‑1 (MQ=255)

                                    tGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGt                                                       <  2:1059316/93‑1 (MQ=255)

                                    tGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGt                                                       <  1:3924796/93‑1 (MQ=255)

                                    tGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGt                                                       <  2:2778294/93‑1 (MQ=255)

                                     gCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTg                                                      <  1:3591307/93‑1 (MQ=255)

                                      caccacTATCAACCACGGCAACTCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGc                                                     <  2:894432/93‑1 (MQ=255)

                                      caccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGc                                                     >  2:1522400/1‑93 (MQ=255)

                                       accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                    <  1:2077662/93‑1 (MQ=255)

                                           cTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTAc                                                <  2:3569095/93‑1 (MQ=255)

                                           cTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTAc                                                <  2:2441011/93‑1 (MQ=255)

                                            tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                               >  1:1750447/1‑93 (MQ=255)

                                            tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                               >  1:3233742/1‑93 (MQ=255)

                                            tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                               >  1:2159409/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCTAGATCAACAGCATATCGCGTGCACTACaa                                              >  1:4196694/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                              >  1:1896415/1‑93 (MQ=255)

                                               cAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcg                                            <  1:3078089/93‑1 (MQ=255)

                                               cAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcg                                            <  2:1287294/93‑1 (MQ=255)

                                                aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCTAGATCAACAGCATATCGCGTGCACTACAAcgg                                           <  1:5233/93‑1 (MQ=255)

                                                aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACCAcgg                                           >  1:593623/1‑93 (MQ=255)

                                                aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcgg                                           <  2:215393/93‑1 (MQ=255)

                                                 aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                          <  1:233838/93‑1 (MQ=255)

                                                 aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                          <  1:233837/93‑1 (MQ=255)

                                                  ccACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcg                                         <  2:3189227/93‑1 (MQ=255)

                                                    aCGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCTCTACAACGGCGGa                                       >  1:3151387/1‑93 (MQ=255)

                                                      ggCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaa                                     <  2:2697125/93‑1 (MQ=255)

                                                        cAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAg                                   <  2:2230856/93‑1 (MQ=255)

                                                         aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGt                                  <  1:3370578/93‑1 (MQ=255)

                                                          aCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAggt                                 >  2:2819261/1‑91 (MQ=255)

                                                          aCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGtt                                 <  2:4448655/93‑1 (MQ=255)

                                                           cGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGgtt                                >  2:3941504/1‑93 (MQ=255)

                                                            gCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTtta                               <  2:1327805/93‑1 (MQ=255)

                                                            gCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTtta                               <  2:4109216/93‑1 (MQ=255)

                                                             ccGTAGTTGTTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATAAACAGCATATCGCGTGCAATACAACAGCGGAAAAGtttct                              >  2:1936469/1‑91 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                              <  1:2084769/93‑1 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                              >  1:1752434/1‑93 (MQ=255)

                                                              cGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                             >  2:1216356/1‑93 (MQ=255)

                                                              cGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                             <  1:1925981/93‑1 (MQ=255)

                                                               gTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAGCAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                            <  2:2336401/93‑1 (MQ=255)

                                                                tAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTTTTAc                           >  1:3722926/1‑93 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                         <  2:4582090/93‑1 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                         <  1:1406468/93‑1 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                         <  1:3784106/93‑1 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                         <  1:854876/93‑1 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                         <  1:563002/93‑1 (MQ=255)

                                                                   ttGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCa                        >  1:4546787/1‑93 (MQ=255)

                                                                     gCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgc                      >  2:1610854/1‑93 (MQ=255)

                                                                          cATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaaa                 <  2:356995/93‑1 (MQ=255)

                                                                          cATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaaa                 >  1:4346680/1‑93 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg                >  2:3919895/1‑93 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg                <  1:1686317/93‑1 (MQ=255)

                                                                            tAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGa               <  1:3637106/93‑1 (MQ=255)

                                                                            tAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGa               <  1:3437813/93‑1 (MQ=255)

                                                                            tAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGa               <  2:34567/93‑1 (MQ=255)

                                                                              aTGCCCCAGCGCGCTACGCCGCATATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATa             <  2:2589043/93‑1 (MQ=255)

                                                                                gCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt           <  1:1959340/93‑1 (MQ=255)

                                                                                gCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt           <  2:4523586/93‑1 (MQ=255)

                                                                                 ccccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATg          <  1:4018201/93‑1 (MQ=255)

                                                                                  cccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGa         <  2:3663388/93‑1 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt      <  2:2994802/93‑1 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTCAGCGCAAAGATAATGACtt      >  2:3172938/1‑93 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc     >  1:2915041/1‑93 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc     <  2:2689887/93‑1 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTCAGCGCAAAGATAATGACTTc     >  1:3378085/1‑93 (MQ=255)

                                                                                       cgcgCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTcc    >  1:3107249/1‑93 (MQ=255)

                                                                                         cACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGc  >  2:2368311/1‑93 (MQ=255)

AGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGC  >  NC_000913/2664941‑2665122

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