BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,533,492	A→T	intergenic (+112/‑272)	cysK → / → ptsH	cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,533,492	0	A	T	94.9%	721.5 / 19.0	218	intergenic (+112/‑272)	cysK/ptsH	cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
Reads supporting (aligned to +/- strand): ref base A (8/3); new base T (105/101); total (113/105)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

AATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGG  >  NC_000913/2533396‑2533587

aaTGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGTAACTGGTATTTAAc                                                                                               >  1:1670408/1‑99 (MQ=255)

aaTGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTAc                                                                                               <  2:3412065/99‑1 (MQ=255)

aaTGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTAc                                                                                               <  1:505292/99‑1 (MQ=255)

aaTGCGTAAAAAAGCACCTTATTAGGCGCTTTTTTGTGGGCTGCTCCAAACTTTCGCCCCTCCTGGCATTCATTCTGCCTGTCGGAACTGGTATTTTAc                                                                                               <  1:5270464/99‑1 (MQ=255)

  tGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCa                                                                                             <  2:279696/99‑1 (MQ=255)

    cGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGa                                                                                           >  2:4247325/1‑99 (MQ=255)

     gTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGAc                                                                                          >  1:3871497/1‑99 (MQ=255)

      tAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACt                                                                                         <  2:1575928/99‑1 (MQ=255)

        aaaaaGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTaa                                                                                       <  1:3219533/99‑1 (MQ=255)

          aaaGCACCTTTTTAGGTGCTTTTTTGTCGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGGATTTTACCAGCCTAatt                                                                                     >  2:25189/1‑99 (MQ=255)

           aaGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAatta                                                                                    >  1:5294681/1‑99 (MQ=255)

           aaGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAatta                                                                                    >  1:602489/1‑99 (MQ=255)

             gCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAattatt                                                                                  <  1:2850767/99‑1 (MQ=255)

               aCCTTTTTAGGTGCTTTTTTTTTTTTTTTTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTCtttt                                                                                >  2:1799217/1‑99 (MQ=255)

               aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTAtttt                                                                                <  2:4774668/99‑1 (MQ=255)

               aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTAtttt                                                                                <  2:4561991/99‑1 (MQ=255)

               aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTAtttt                                                                                <  1:4426449/99‑1 (MQ=255)

                ccTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTg                                                                               <  2:2740040/99‑1 (MQ=255)

                ccTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTg                                                                               >  1:1336754/1‑99 (MQ=255)

                 cTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGa                                                                              <  2:4453600/99‑1 (MQ=255)

                  tttttAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGAt                                                                             <  2:843830/99‑1 (MQ=255)

                   aaaaaGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATg                                                                            >  2:3299320/5‑99 (MQ=255)

                    tttAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATgc                                                                           <  2:5211181/99‑1 (MQ=255)

                      tAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATgcgc                                                                         <  1:1629097/99‑1 (MQ=255)

                         gTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTTATGCGCGaa                                                                      >  1:1548601/1‑99 (MQ=255)

                         gTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGaa                                                                      >  1:1404628/1‑99 (MQ=255)

                          tGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGaaa                                                                     <  2:2598857/99‑1 (MQ=255)

                           gCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAAt                                                                    >  1:3595285/1‑99 (MQ=255)

                               ttttGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAAt                                                                <  2:3948082/99‑1 (MQ=255)

                               ttttGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAAt                                                                <  2:1253024/99‑1 (MQ=255)

                                 ttGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCg                                                              <  1:881287/99‑1 (MQ=255)

                                  tgtgGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGt                                                             <  1:2292691/99‑1 (MQ=255)

                                    tgGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTa                                                           >  1:841168/1‑99 (MQ=255)

                                    tgGCCTGCTTCAAACGTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGTTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTa                                                           <  2:4369251/99‑1 (MQ=255)

                                     ggCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTAc                                                          <  1:4664846/99‑1 (MQ=255)

                                     ggCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTAc                                                          <  2:2906861/99‑1 (MQ=255)

                                      gCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACa                                                         <  1:58020/99‑1 (MQ=255)

                                       ccTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAg                                                        >  2:2346148/1‑99 (MQ=255)

                                        cTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCCGACTAATTATTTTGATGCGCGAAATTAATCGTTACAgg                                                       >  1:3838021/1‑99 (MQ=255)

                                           cTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaa                                                    >  2:5289573/1‑99 (MQ=255)

                                            ttaaaaCTTTCGCCCCTCCTGGCATTGATTAAGCCTGTCGGAACTGGTATTTTACCATACTCATTCGTTTGATGCGCGAAATTAATCGTTACACGaaaa                                                   >  2:1062114/4‑99 (MQ=255)

                                            ttaaaaCTTTCGCCCCTCCCGGCATTGATTCACCCTGTCGCAACTGGAATTTAACCAGAACAATTCTTTTGTTGCGCGAAATTCATTGCTACTGGaaaa                                                   >  2:1140677/4‑99 (MQ=255)

                                            ttCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaaa                                                   >  2:1180633/1‑99 (MQ=255)

                                            ttCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaaa                                                   >  1:1426699/1‑99 (MQ=255)

                                            ttCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaaa                                                   >  2:4844147/1‑99 (MQ=255)

                                            ttCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTCTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGAGCGAAATTAATCGTTACAGGaaaa                                                   >  2:3785914/1‑99 (MQ=255)

                                               aaaCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGcc                                                <  1:2592760/99‑1 (MQ=255)

                                               aaaCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGcc                                                >  1:1911282/1‑99 (MQ=255)

                                                aaCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCa                                               <  1:3541529/99‑1 (MQ=255)

                                                aaCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCa                                               <  2:421726/99‑1 (MQ=255)

                                                aaCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCa                                               >  2:2073610/1‑99 (MQ=255)

                                                 aCTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCACGAAAGTAATCGTTACAGGAAAAGCCaa                                              >  1:5125313/1‑99 (MQ=255)

                                                  cTTTCGCCCCTCCTGGCATTGATTCAGCCTTTCGAACCTGGTATTTTACCAGACTAATTATTTTGATGCGCGAGATTAATCGTTACAGGAAAAGCCaaa                                             >  2:5383505/1‑99 (MQ=255)

                                                   tttCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAg                                            >  2:3224032/1‑99 (MQ=255)

                                                   tttCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAg                                            >  1:3493562/1‑99 (MQ=255)

                                                   tttCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAg                                            <  1:1220451/99‑1 (MQ=255)

                                                     tCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCt                                          <  1:4891455/99‑1 (MQ=255)

                                                         cccTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAGATTAATCGTTACAGGAAAAGCCAAAGCTGAAt                                      >  1:5172101/1‑99 (MQ=255)

                                                         cccTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt                                      <  1:1977850/99‑1 (MQ=255)

                                                         cccTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt                                      >  1:4235218/1‑99 (MQ=255)

                                                         cccGCCTGGCATTAATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATATTGATGCGCGAAATTAATCTTTACAGGAAAAGCCAAAGCTGAAt                                      >  2:3522939/1‑99 (MQ=255)

                                                          cctcctGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATc                                     >  2:2690133/1‑99 (MQ=255)

                                                          cctcctGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATc                                     <  2:946158/99‑1 (MQ=255)

                                                           ctcctGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCg                                    >  1:3483300/1‑99 (MQ=255)

                                                             cctGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAt                                  <  1:4313813/99‑1 (MQ=255)

                                                              ctGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAtt                                 <  1:5645752/99‑1 (MQ=255)

                                                               tGGCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAttt                                >  1:1047467/1‑99 (MQ=255)

                                                                ggCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAtttt                               >  1:4993350/1‑99 (MQ=255)

                                                                 gCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTa                              <  1:1158553/99‑1 (MQ=255)

                                                                 gCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTa                              <  2:3079378/99‑1 (MQ=255)

                                                                 gCATTGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTa                              >  2:5050815/1‑99 (MQ=255)

                                                                  cATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTAt                             >  2:3853682/1‑99 (MQ=255)

                                                                    ttGATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGa                           <  1:3755892/99‑1 (MQ=255)

                                                                      gATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAtt                         <  1:2231305/99‑1 (MQ=255)

                                                                      gATTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAtt                         >  2:4638903/1‑99 (MQ=255)

                                                                       aTTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAttt                        >  2:3213397/1‑99 (MQ=255)

                                                                       aTTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAttt                        >  2:4926412/1‑99 (MQ=255)

                                                                       aTTCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAttt                        >  1:203159/1‑99 (MQ=255)

                                                                        ttCAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTg                       <  1:4583248/99‑1 (MQ=255)

                                                                          cAGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGt                     <  1:4367242/99‑1 (MQ=255)

                                                                           aGCCTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGtt                    <  2:5085140/99‑1 (MQ=255)

                                                                             ccTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAATCCAAAGCTGTATCGATTTTATGATTTGGTTCa                  >  2:3608277/1‑99 (MQ=255)

                                                                              cTGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCaa                 >  1:388334/1‑99 (MQ=255)

                                                                               tGTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAt                <  1:1942199/99‑1 (MQ=255)

                                                                                gTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAtt               >  2:3701799/1‑99 (MQ=255)

                                                                                gTCGGAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAtt               >  1:1243052/1‑99 (MQ=255)

                                                                                  cGGAACTCGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAttct             <  1:4745756/99‑1 (MQ=255)

                                                                                    gAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAttcttc           <  2:675053/99‑1 (MQ=255)

                                                                                    gAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAttcttc           >  1:1770728/1‑99 (MQ=255)

                                                                                    gAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAttcttc           >  1:4132449/1‑99 (MQ=255)

                                                                                    gAACTGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAttcttc           <  2:34353/99‑1 (MQ=255)

                                                                                        tGGTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCttt       >  2:76855/1‑99 (MQ=255)

                                                                                         ggTAATTTACCAGACTAATTATCTAGAGGCACGAAATTAATCGTTACAGGAACAGCCAAAGCCGAATCGATTTTATGATTTGGTTCAATTCGTCCTTTa      <  2:3561107/99‑1 (MQ=255)

                                                                                          gTATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAg     <  1:1854376/99‑1 (MQ=255)

                                                                                           tATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGc    >  1:2981936/1‑99 (MQ=255)

                                                                                           tATTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGc    >  1:3386576/1‑99 (MQ=255)

                                                                                            aTTTTACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCg   >  1:5448864/1‑99 (MQ=255)

                                                                                            aTTTTACCAGACTAATTATTTGGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCg   <  2:3750237/99‑1 (MQ=255)

                                                                                             ttttACCAGACTGATTATTTTGATGCACGAAATTAATCGTTACAGGAAAAGCCAAAGCAGAATCGATTTTATGATTTGGTTCAATTCGTCCTTTAGCgg  >  1:4967084/1‑99 (MQ=255)

                                                                                             ttttACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCGAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCgg  >  2:5194413/1‑99 (MQ=255)

AATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCCTTTAGCGG  >  NC_000913/2533396‑2533587

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: `‑`