BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,974,806	C→T	A306V (GCG→GTG)	wzxE →	putative cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,974,806	0	C	T	90.6%	704.1 / 44.5	235	A306V (GCG→GTG)	wzxE	putative cytochrome
Reads supporting (aligned to +/- strand): ref base C (13/9); new base T (94/118); total (108/127)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

Only 100 of 269 total aligned reads displayed.

TGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCTTTGCCT  >  NC_000913/3974716‑3974895

tgctgcCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTggcggtgg                                                                                         >  2:2435329/1‑93 (MQ=255)

    gcCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAg                                                                                     >  2:1367904/1‑93 (MQ=255)

     cccACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGCTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGt                                                                                    >  1:6034031/1‑93 (MQ=255)

      ccACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTTGCGTGtt                                                                                   >  1:1235093/1‑93 (MQ=255)

      ccACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAggt                                                                                   >  1:4640546/1‑91 (MQ=255)

      ccACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAggt                                                                                   >  1:2852535/1‑91 (MQ=255)

      ccACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGtt                                                                                   <  1:3559936/93‑1 (MQ=255)

       cACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGttt                                                                                  >  1:4547398/1‑93 (MQ=255)

       cACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGttt                                                                                  <  1:2688872/93‑1 (MQ=255)

        aCGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGtttt                                                                                 >  1:248283/1‑93 (MQ=255)

          gTTGTCGCTGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGG‑TGCGAGTTTTAcc                                                                              >  2:1002463/1‑93 (MQ=255)

          gTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGTAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTAc                                                                               >  2:385591/1‑93 (MQ=255)

           ttGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGGTTTAcc                                                                              >  2:4542760/1‑93 (MQ=255)

            tGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCg                                                                             >  1:5083627/1‑93 (MQ=255)

            tGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGGTTTACCg                                                                             >  1:5616042/1‑93 (MQ=255)

               cgcgGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCTCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCt                                                                          <  2:4170449/93‑1 (MQ=255)

               cgcgGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTTGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCt                                                                          <  1:4953069/93‑1 (MQ=255)

                gcgGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTg                                                                         <  1:3233612/93‑1 (MQ=255)

                gcgGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTg                                                                         <  1:5447020/93‑1 (MQ=255)

                 cgGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTgg                                                                        >  2:5134901/1‑93 (MQ=255)

                 cgGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTgg                                                                        >  2:1488955/1‑93 (MQ=255)

                  ggCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgc                                                                       <  1:2287186/93‑1 (MQ=255)

                  ggCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgc                                                                       <  2:6122738/93‑1 (MQ=255)

                  ggCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgc                                                                       >  1:4450821/1‑93 (MQ=255)

                  ggCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGgc                                                                       <  2:3483198/93‑1 (MQ=255)

                   gCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgct                                                                      >  2:2700072/1‑93 (MQ=255)

                   gCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgct                                                                      <  2:164581/93‑1 (MQ=255)

                   gCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgct                                                                      >  2:1737053/1‑93 (MQ=255)

                   gCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGGTTTACCGTCTGgct                                                                      >  2:3718552/1‑93 (MQ=255)

                    cTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgctg                                                                     >  2:5294946/1‑93 (MQ=255)

                    cTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGGTTTACCGTCTGgctg                                                                     >  2:5831097/1‑93 (MQ=255)

                     tAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgctgc                                                                    >  1:4300480/1‑93 (MQ=255)

                     tAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgctgc                                                                    <  1:88741/93‑1 (MQ=255)

                     tAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgctgc                                                                    >  1:3674888/1‑93 (MQ=255)

                     tAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgctgc                                                                    <  2:4835967/93‑1 (MQ=255)

                      aaCGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGgctgct                                                                   >  1:638103/1‑93 (MQ=255)

                       aCGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGctgctg                                                                  >  1:2140573/1‑93 (MQ=255)

                        cGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCtgctgc                                                                 <  2:323369/93‑1 (MQ=255)

                         ggAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCg                                                                >  2:735210/1‑93 (MQ=255)

                          gAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGt                                                               <  2:2476285/93‑1 (MQ=255)

                            aaaGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGa                                                             <  1:6294662/93‑1 (MQ=255)

                             aaGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGAt                                                            <  2:6561191/93‑1 (MQ=255)

                              aGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGAtt                                                           >  1:4199985/1‑93 (MQ=255)

                               gcgcGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGAttt                                                          <  1:945590/93‑1 (MQ=255)

                                 gcgATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTg                                                        >  1:5958279/1‑93 (MQ=255)

                                  cgATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGc                                                       <  1:1415677/93‑1 (MQ=255)

                                  cgATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGc                                                       <  2:4973160/93‑1 (MQ=255)

                                   gATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCt                                                      >  2:4252359/1‑93 (MQ=255)

                                   gATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCCGTGGCGAGGTTTACCGTCTGGCTGCTGCGTGATTTTGCt                                                      >  2:3157881/1‑93 (MQ=255)

                                    gtatCCCCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTa                                                     <  2:2981014/92‑1 (MQ=255)

                                      atCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATc                                                   <  2:3395638/93‑1 (MQ=255)

                                      atCACCCGGGAAGTGCTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATc                                                   <  1:1597314/93‑1 (MQ=255)

                                        cACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTg                                                 <  2:650113/93‑1 (MQ=255)

                                          cccGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTGACCGTCTGGCTGCTGCGTGATTTTGCTATCTGgc                                               <  1:3382364/93‑1 (MQ=255)

                                           ccGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGgct                                              >  2:3538415/1‑93 (MQ=255)

                                               gAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGctgctg                                          >  2:4046561/1‑93 (MQ=255)

                                                aaGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCtgt                                         <  2:4635020/93‑1 (MQ=255)

                                                 aGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGGGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCtgtt                                        <  1:5452439/93‑1 (MQ=255)

                                                  gTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCtgttg                                       >  2:2579717/1‑93 (MQ=255)

                                                  gTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCtgttg                                       <  2:5405649/93‑1 (MQ=255)

                                                     gTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCg                                    <  1:1112431/93‑1 (MQ=255)

                                                     gTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCg                                    <  1:3461903/93‑1 (MQ=255)

                                                      ttAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGa                                   <  1:3281087/93‑1 (MQ=255)

                                                      ttAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGa                                   >  1:414987/1‑93 (MQ=255)

                                                       tAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGaa                                  <  2:2758661/93‑1 (MQ=255)

                                                       tAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGaa                                  >  2:4688711/1‑93 (MQ=255)

                                                        aaaTCTCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGGGATTTTGCTATCTGGCTGCTGTTGTCGAAt                                 <  2:627483/93‑1 (MQ=255)

                                                          aTCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATaa                               <  2:2173898/93‑1 (MQ=255)

                                                           tCGCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATaaa                              <  2:1948297/93‑1 (MQ=255)

                                                           tCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATaaa                              >  1:215497/1‑93 (MQ=255)

                                                             gCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAAtt                            >  1:164292/1‑93 (MQ=255)

                                                             gCTGAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAAtt                            >  2:3604120/1‑93 (MQ=255)

                                                                gAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTAc                         >  2:798194/1‑93 (MQ=255)

                                                                gAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTAc                         <  2:5552291/93‑1 (MQ=255)

                                                                gAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTAc                         <  1:3266290/93‑1 (MQ=255)

                                                                gAAATTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTAc                         <  2:4891729/93‑1 (MQ=255)

                                                                  aaTTCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCg                       >  1:3931756/1‑93 (MQ=255)

                                                                   aTTCGTCTTACCGGCAGTGTCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGc                      <  1:5223278/93‑1 (MQ=255)

                                                                     tCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTa                    <  2:2881857/93‑1 (MQ=255)

                                                                     tCGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTa                    <  2:5173182/93‑1 (MQ=255)

                                                                      cGTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTAt                   <  2:2037705/93‑1 (MQ=255)

                                                                       gTCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATg                  >  2:4231693/1‑93 (MQ=255)

                                                                        tCTTACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATgc                 >  2:569626/1‑93 (MQ=255)

                                                                          ttACCGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATgcgc               >  2:1829245/1‑93 (MQ=255)

                                                                              cGGCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGAtc           >  2:1161702/1‑93 (MQ=255)

                                                                               ggCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGAtct          >  1:3692367/1‑93 (MQ=255)

                                                                               ggCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGAtct          <  1:1520248/93‑1 (MQ=255)

                                                                                gCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGAtctc         >  1:5116012/1‑93 (MQ=255)

                                                                                gCAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGAtctc         <  2:5754285/93‑1 (MQ=255)

                                                                                 cAGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATctct        <  2:6101255/93‑1 (MQ=255)

                                                                                 cAGTGGAGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATctct        <  1:679811/93‑1 (MQ=255)

                                                                                  aGTGGCGGTGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCtt       <  1:1777409/93‑1 (MQ=255)

                                                                                   gtggcggtggcgAGTTTTATCGTCTGGCTGCTGCGTGGTTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCttt      >  1:6490509/1‑93 (MQ=255)

                                                                                   gtggcggtggcgAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCttt      <  1:2588647/93‑1 (MQ=255)

                                                                                   gtggcggtggcgAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCttt      >  2:3772289/1‑93 (MQ=255)

                                                                                   gTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCttt      >  2:2298345/1‑93 (MQ=255)

                                                                                     ggcggtggcGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCTTTGc    <  1:1852205/93‑1 (MQ=255)

                                                                                       cggtggcgAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCTTTGCCt  <  1:750687/93‑1 (MQ=255)

                                                                                       cggtggcgAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCTTTGCCt  <  2:3972807/93‑1 (MQ=255)

                                                                                       cggtggcgAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCTTTGCCt  >  1:1041536/1‑93 (MQ=255)

TGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATCACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGGCGGCGGCGAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAATTTACCGCTATGCGCGATCTCTTTGCCT  >  NC_000913/3974716‑3974895

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