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.7%	1364.1 / ‑2.9	301	coding (541/609 nt)	b4556	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (163/137); minor base A (1/0); total (164/137)
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.63e-01

Only 100 of 338 total aligned reads displayed.

TTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTC  >  NC_000913/4021591‑4021771

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGGATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          <  1:4174857/93‑1 (MQ=255)

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          <  1:2924140/93‑1 (MQ=255)

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          >  2:3125199/1‑93 (MQ=255)

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          >  1:4505952/1‑93 (MQ=255)

 tcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcg                                                                                         >  2:6306899/1‑92 (MQ=255)

  cttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcga                                                                                        <  1:3332959/93‑3 (MQ=255)

  cttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcga                                                                                        >  1:4061657/1‑91 (MQ=255)

   ttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                                       <  1:97463/93‑4 (MQ=255)

   ttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                                       >  2:4847814/1‑90 (MQ=255)

    tcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgacg                                                                                      >  1:3066416/1‑89 (MQ=255)

     cATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgc                                                                                    <  2:6984685/93‑1 (MQ=255)

      aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgcg                                                                                   <  2:1178969/93‑1 (MQ=255)

      aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgcg                                                                                   >  1:3517448/1‑93 (MQ=255)

      aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgcg                                                                                   <  1:3697648/93‑1 (MQ=255)

       ttGCTTTACCGGCACACGCGAAGAGATTCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCgg                                                                                  <  1:4467702/93‑1 (MQ=255)

         gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                <  1:2252355/93‑1 (MQ=255)

         gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                <  2:7437517/93‑1 (MQ=255)

          cTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAcg                                                                               >  1:4885888/1‑92 (MQ=255)

          cTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACt                                                                               >  2:4385191/1‑93 (MQ=255)

             tACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGa                                                                            <  1:4728980/93‑1 (MQ=255)

             tACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGa                                                                            <  2:1799948/93‑1 (MQ=255)

                cGGCACACGCGAAGAGATGCAGGCGTGCGGGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCt                                                                         >  2:3311140/1‑93 (MQ=255)

                  gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGc                                                                       <  2:1976535/93‑1 (MQ=255)

                      acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGa                                                                   >  1:5503853/1‑93 (MQ=255)

                       cgcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGaa                                                                  <  2:6761894/93‑1 (MQ=255)

                        gcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAAc                                                                 <  2:5525744/93‑1 (MQ=255)

                          gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                                               >  2:2177307/1‑93 (MQ=255)

                           aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGGACttt                                                              >  1:5847865/1‑93 (MQ=255)

                           aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACttt                                                              >  1:267292/1‑93 (MQ=255)

                           aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACttt                                                              <  1:5465078/93‑1 (MQ=255)

                              agaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                                           <  2:5142633/93‑1 (MQ=255)

                                aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                                         >  2:3486193/1‑93 (MQ=255)

                                 tGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGt                                                        >  1:1622840/1‑93 (MQ=255)

                                 tGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGt                                                        >  1:51722/1‑93 (MQ=255)

                                  gCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                                       <  2:6039383/93‑1 (MQ=255)

                                   cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTg                                                      >  1:1450082/1‑93 (MQ=255)

                                    aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                     <  2:5919480/93‑1 (MQ=255)

                                    aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                     <  2:7827081/93‑1 (MQ=255)

                                    aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                     <  1:5687079/93‑1 (MQ=255)

                                     ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                                    >  2:274509/1‑93 (MQ=255)

                                     ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                                    >  1:5419797/1‑93 (MQ=255)

                                      gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                                   <  2:1078796/93‑1 (MQ=255)

                                      gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                                   <  2:2523292/93‑1 (MQ=255)

                                          gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcgg                                               >  1:2542329/1‑93 (MQ=255)

                                           cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                              <  1:4376330/93‑1 (MQ=255)

                                                 gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                        >  2:3992534/1‑93 (MQ=255)

                                                 gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                        >  2:3246015/1‑93 (MQ=255)

                                                  gaatgAATTTATTACGGCCATACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                       <  1:1752534/89‑1 (MQ=255)

                                                  cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                       >  1:3762393/1‑93 (MQ=255)

                                                  cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                       <  2:7328269/93‑1 (MQ=255)

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                     >  2:5042924/1‑93 (MQ=255)

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                     >  1:6754658/1‑93 (MQ=255)

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                     >  2:7114389/1‑93 (MQ=255)

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAta                                     >  1:784310/1‑92 (MQ=255)

                                                     ggAATTTATATCGGCATTACCGGTTGGGTTTGCGTTGAA‑CGACGCGGACTGGAGCTGCGGGAACTCTTGCCGCTGATTCCGGCGGAAAAATTa                                    >  2:2652682/1‑93 (MQ=255)

                                                     ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                                    <  1:1565411/93‑1 (MQ=255)

                                                      gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                   >  1:1158170/1‑93 (MQ=255)

                                                      gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                   <  1:7534147/93‑1 (MQ=255)

                                                       aaTTTATATCGGCATTACCGGTTGTGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCTGGAACTTTTGTCGTTGATTCCGGCGGAAAAATTACt                                  >  1:6926961/1‑93 (MQ=255)

                                                       aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                                  <  2:6104652/93‑1 (MQ=255)

                                                        aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                 <  1:701469/93‑1 (MQ=255)

                                                         tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                >  1:4579437/1‑93 (MQ=255)

                                                         tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                >  2:5718345/1‑93 (MQ=255)

                                                          ttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                               <  1:521498/93‑1 (MQ=255)

                                                            atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGGAAAATTACTGATCg                             >  1:1050405/1‑93 (MQ=255)

                                                            atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg                             >  2:959660/1‑93 (MQ=255)

                                                               tCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGgaa                          <  2:6387753/93‑1 (MQ=255)

                                                                cGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAc                         <  2:1488297/93‑1 (MQ=255)

                                                                   cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                      <  2:200348/93‑1 (MQ=255)

                                                                   cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                      >  1:6098082/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  2:359444/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  2:7216449/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  2:1404244/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  1:7809211/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  1:6449523/1‑93 (MQ=255)

                                                                      tACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc                   >  1:716362/1‑93 (MQ=255)

                                                                        ccGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcgc                 <  1:7472238/93‑1 (MQ=255)

                                                                        ccGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcgc                 >  1:2004434/1‑93 (MQ=255)

                                                                        ccGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcgc                 <  2:560195/93‑1 (MQ=255)

                                                                           gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGt              >  2:329044/1‑93 (MQ=255)

                                                                           gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGt              >  1:38875/1‑93 (MQ=255)

                                                                           gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGt              >  1:842318/1‑93 (MQ=255)

                                                                           gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGt              >  2:4076983/1‑93 (MQ=255)

                                                                            ttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                      >  2:1672082/1‑84 (MQ=255)

                                                                            ttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTa             <  2:6798727/93‑1 (MQ=255)

                                                                              gggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc           <  2:4990736/93‑1 (MQ=255)

                                                                               ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt          >  2:2904351/1‑93 (MQ=255)

                                                                               ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt          <  1:1218983/93‑1 (MQ=255)

                                                                               ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt          >  1:626153/1‑93 (MQ=255)

                                                                               ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAATGATGCGCCGTATCt          <  2:5326028/93‑1 (MQ=255)

                                                                                gTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTg         >  1:4779141/1‑93 (MQ=255)

                                                                                 tttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGc        <  1:3589034/93‑1 (MQ=255)

                                                                                  ttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCt       >  2:5823873/1‑93 (MQ=255)

                                                                                   tGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc      <  1:3245303/93‑1 (MQ=255)

                                                                                   tGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc      >  1:7371621/1‑93 (MQ=255)

                                                                                    gCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTcc     >  2:146838/1‑93 (MQ=255)

                                                                                    gCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTcc     <  2:4285074/93‑1 (MQ=255)

                                                                                      gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCt   <  2:4585788/93‑1 (MQ=255)

                                                                                      gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCt   >  1:1857264/1‑93 (MQ=255)

                                                                                       aTGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTc  <  2:3829614/93‑1 (MQ=255)

TTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTC  >  NC_000913/4021591‑4021771

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