BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,021,684	Δ1 bp	coding (542/609 nt)\|coding (1/243 nt)	b4556 → tatD →	hypothetical protein Mg‑dependent DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,021,683	0	C	.	99.6%	2233.7 / ‑0.6	467	coding (541/609 nt)	b4556	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base . (189/276); total (191/276)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 504 total aligned reads displayed.

CGGTTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCG  >  NC_000913/4021588‑4021774

cGGTTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                             <  1:9799719/96‑1 (MQ=255)

 ggTTTTTCATTTCTTTACCGGCAAACGCGAAGAGAGGCGGGCGTGCGTGGCGCATGGAATTTATATCGGCA‑TACCGGTTGGGTTTGCGATGAAcga                                                                                           <  2:3216857/96‑3 (MQ=255)

 ggTTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcg                                                                                            >  1:7768226/1‑95 (MQ=255)

 ggTTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGGTTGCGATGAAcg                                                                                            >  1:2451161/1‑95 (MQ=255)

   ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                                          >  1:3482411/1‑93 (MQ=255)

    tcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgacg                                                                                         <  2:4989236/96‑5 (MQ=255)

        cATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGa                                                                                    <  2:8362308/96‑1 (MQ=255)

         aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                   <  2:6042381/96‑1 (MQ=255)

           tGCTTTGCCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                                                 <  2:7660541/96‑1 (MQ=255)

            gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTgg                                                                                <  1:4869339/96‑1 (MQ=255)

             cTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGa                                                                               <  1:2582979/96‑1 (MQ=255)

                 aCCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTg                                                                           <  2:7550813/96‑1 (MQ=255)

                  ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGc                                                                          <  2:5491509/96‑1 (MQ=255)

                   cGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCg                                                                         >  1:100041/1‑96 (MQ=255)

                   cGGCACACGCGAAGAGATGCAGGCGTGCGTGGCCCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCg                                                                         <  1:686943/96‑1 (MQ=255)

                    ggCACACGCGAAGAGATGCAGGCGTGCGTGGGGCATGGAAGTTATATCGGCATTACCGGTTGGGTTTCCGATGAA‑CGACGCGGACTGGAGCTGCgg                                                                        <  2:13911792/96‑1 (MQ=255)

                    ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCgg                                                                        <  1:1439086/96‑1 (MQ=255)

                      cacaCGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGa                                                                      <  2:7134205/96‑1 (MQ=255)

                       acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGAGTTTGCGATGGA‑CGACGCGGACTGGAGCTGCGGGaa                                                                     >  1:12703596/1‑96 (MQ=255)

                        cacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAAc                                                                    <  1:12035429/96‑1 (MQ=255)

                         acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACt                                                                   >  2:12658430/1‑96 (MQ=255)

                          cgcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                                                  <  2:960558/96‑1 (MQ=255)

                            cgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtttt                                                                <  1:11017249/96‑1 (MQ=255)

                             gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTg                                                               >  2:8725793/1‑96 (MQ=255)

                              aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                                              <  2:6409014/96‑1 (MQ=255)

                              aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                                              <  1:12265820/96‑1 (MQ=255)

                               agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGcc                                                             <  2:2942176/96‑1 (MQ=255)

                                gagaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGACCTTTTGCCg                                                            <  1:4071488/96‑1 (MQ=255)

                                 agaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGt                                                           >  2:1036236/1‑96 (MQ=255)

                                  gaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                                          >  1:2200863/1‑96 (MQ=255)

                                     gCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                                       >  1:12484105/1‑96 (MQ=255)

                                      cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                                      <  1:5947321/96‑1 (MQ=255)

                                        ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                                    <  1:8334167/96‑1 (MQ=255)

                                         gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                   <  1:8804689/96‑1 (MQ=255)

                                         gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                   <  1:7378658/96‑1 (MQ=255)

                                         gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                   <  2:4205984/96‑1 (MQ=255)

                                          cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcgg                                                  <  2:12060563/96‑1 (MQ=255)

                                          cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcgg                                                  <  2:3584637/96‑1 (MQ=255)

                                          cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTAGCCGTTGATTCcgg                                                  <  1:6679829/96‑1 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                                 <  2:12205991/96‑1 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                                 <  2:10691894/96‑1 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                                 >  1:13114059/1‑96 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                                 <  2:9845957/96‑1 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                                 <  2:14216576/96‑1 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTGGATTCcggc                                                 >  1:9046311/1‑96 (MQ=255)

                                           gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGATGGGTTTGCGATG‑AGCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                                 >  1:1941742/1‑96 (MQ=255)

                                              cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                                              >  2:13556681/1‑96 (MQ=255)

                                              cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTGTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                                              >  1:10162597/1‑96 (MQ=255)

                                               gtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaa                                             >  2:7997312/1‑96 (MQ=255)

                                                tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaa                                            >  1:14109991/1‑96 (MQ=255)

                                                tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaa                                            >  2:1444999/1‑96 (MQ=255)

                                                 ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                           <  2:9208878/96‑1 (MQ=255)

                                                  gcccATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                          <  1:13907549/93‑1 (MQ=255)

                                                     cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                                       >  2:3831691/1‑96 (MQ=255)

                                                     cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                                       <  2:1594940/96‑1 (MQ=255)

                                                       tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                                     <  2:11500288/96‑1 (MQ=255)

                                                        ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACta                                    <  1:12941092/96‑2 (MQ=255)

                                                        ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                    <  2:11051124/96‑1 (MQ=255)

                                                         gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                   >  2:9713237/1‑96 (MQ=255)

                                                         gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                   <  2:229878/96‑1 (MQ=255)

                                                          aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                                  <  1:4389915/96‑1 (MQ=255)

                                                          aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                                  <  1:13209370/96‑1 (MQ=255)

                                                            tttATATCGGCATTACCGGTTGGGTTTGCGATGGA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg                                >  1:5022108/1‑96 (MQ=255)

                                                            tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg                                >  2:5774581/1‑96 (MQ=255)

                                                             ttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAGTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGacaga                               >  1:4311012/1‑92 (MQ=255)

                                                               atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaga                             <  1:11730468/96‑3 (MQ=255)

                                                               atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaaa                             <  1:8567193/96‑1 (MQ=255)

                                                                 atCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACt                           <  1:13313166/96‑1 (MQ=255)

                                                                 atCGGCATTACCGGGTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACt                           <  1:2900086/96‑1 (MQ=255)

                                                                   cGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                         <  2:4772838/96‑1 (MQ=255)

                                                                    ggCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTgct                        <  1:10649068/96‑3 (MQ=255)

                                                                     gCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                       >  2:11749546/1‑96 (MQ=255)

                                                                     gCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                       <  2:4911306/96‑1 (MQ=255)

                                                                      cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc                      <  2:13070739/96‑1 (MQ=255)

                                                                      cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc                      >  2:845921/1‑96 (MQ=255)

                                                                      cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc                      >  2:9359061/1‑96 (MQ=255)

                                                                      cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCAGGCGGAAAAATTACTGATCGAAACTGATgc                      <  1:5442939/96‑1 (MQ=255)

                                                                       tttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcg                     <  2:7416129/95‑1 (MQ=255)

                                                                        ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcgc                    >  1:4667591/1‑96 (MQ=255)

                                                                         tACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGcc                   <  1:6242725/96‑1 (MQ=255)

                                                                            cGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTa                >  1:451282/1‑96 (MQ=255)

                                                                             gtttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTAt               >  1:10500472/3‑96 (MQ=255)

                                                                             ggTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTAt               >  1:13584629/1‑96 (MQ=255)

                                                                              gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc              <  1:1564104/96‑1 (MQ=255)

                                                                              gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc              >  2:7077938/1‑96 (MQ=255)

                                                                               ttGTGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt             <  1:9608184/96‑1 (MQ=255)

                                                                               ttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt             <  2:4360709/96‑1 (MQ=255)

                                                                               ttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt             <  1:11900636/96‑1 (MQ=255)

                                                                                tGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTg            <  2:10911961/96‑1 (MQ=255)

                                                                                  ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCt          >  1:7562005/1‑96 (MQ=255)

                                                                                   gTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGCTGCGCCGTATCTGCTc         >  1:3860032/1‑96 (MQ=255)

                                                                                   gTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc         <  2:12291518/96‑1 (MQ=255)

                                                                                   gTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc         <  1:794474/96‑1 (MQ=255)

                                                                                    tttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCtac        >  1:9350048/1‑94 (MQ=255)

                                                                                    tttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTcc        >  2:669939/1‑96 (MQ=255)

                                                                                       gCGATGAG‑CGACGCGGACTGTAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGTTCGAACCTGATGCGCCGTATCTGCTCCCTc     <  1:571398/96‑1 (MQ=255)

                                                                                       gCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCtt     <  1:11324733/96‑2 (MQ=255)

                                                                                       gCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTc     <  2:1884636/96‑1 (MQ=255)

                                                                                         gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTcgc   <  2:11261688/96‑1 (MQ=255)

                                                                                         gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTcgc   <  2:10031255/96‑1 (MQ=255)

CGGTTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCG  >  NC_000913/4021588‑4021774

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`