BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,417,173	A→G	S8S (AGT→AGC)	b2299 ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,417,173	0	A	G	84.6%	582.3 / 94.8	253	S8S (AGT→AGC)	b2299	putative regulator
Reads supporting (aligned to +/- strand): ref base A (26/13); new base G (95/119); total (121/132)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 254 total aligned reads displayed.

ACGATGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGATTAAT  >  NC_000913/2417083‑2417265

aCGATGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCt                                                                                            >  2:759644/1‑93 (MQ=255)

aCGATGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCt                                                                                            <  2:3446273/93‑1 (MQ=255)

 cGATGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCACTCAGTGCtt                                                                                           <  2:5488592/93‑1 (MQ=255)

   aTGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGc                                                                                         >  1:4954105/1‑93 (MQ=255)

   aTGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGc                                                                                         >  1:2907518/1‑93 (MQ=255)

       cGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCaaa                                                                                     <  1:6965754/93‑1 (MQ=255)

       cGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCaaa                                                                                     >  2:2389480/1‑93 (MQ=255)

        gCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAac                                                                                    <  2:222067/93‑1 (MQ=255)

        gCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAac                                                                                    <  2:6164253/93‑1 (MQ=255)

          aGACACTGTGCCCGCGTTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAacga                                                                                  <  2:1052142/93‑1 (MQ=255)

           gACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAacgac                                                                                 >  1:7023936/1‑93 (MQ=255)

            acacTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAacgacg                                                                                >  2:2965386/1‑93 (MQ=255)

            acacTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAacgacg                                                                                >  1:411362/1‑93 (MQ=255)

             cacTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGc                                                                               >  1:6657382/1‑93 (MQ=255)

               cTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTg                                                                             >  1:256149/1‑93 (MQ=255)

               cTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTg                                                                             >  2:3026295/1‑93 (MQ=255)

               cTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTg                                                                             >  1:688646/1‑93 (MQ=255)

                tgtgCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGt                                                                            <  2:6707016/93‑1 (MQ=255)

                  tgCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTc                                                                          <  1:3753673/93‑1 (MQ=255)

                  tgCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTc                                                                          >  2:5903311/1‑93 (MQ=255)

                   gCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTcc                                                                         <  1:5437097/93‑1 (MQ=255)

                    cccGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCa                                                                        <  2:3366229/93‑1 (MQ=255)

                     ccGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCAt                                                                       >  1:3374918/1‑93 (MQ=255)

                     ccGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCAt                                                                       >  2:265471/1‑93 (MQ=255)

                      cGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCAtt                                                                      >  2:74423/1‑93 (MQ=255)

                      cGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCAtt                                                                      >  2:2486203/1‑93 (MQ=255)

                          tttGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCa                                                                  >  1:6769821/1‑93 (MQ=255)

                          tttGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCa                                                                  >  1:4713123/1‑93 (MQ=255)

                          tttGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCa                                                                  >  1:1238738/1‑93 (MQ=255)

                           ttGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAg                                                                 <  1:1744014/93‑1 (MQ=255)

                             gTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGa                                                               <  2:494922/93‑1 (MQ=255)

                             gTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGa                                                               <  1:7138364/93‑1 (MQ=255)

                              ttCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGGTTGCCAAACGACGCTGTTCCATTATCAGGaa                                                              >  1:4492802/1‑93 (MQ=255)

                                 ccGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAAcc                                                           <  2:1160656/93‑1 (MQ=255)

                                 ccGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAAcc                                                           <  1:2248653/93‑1 (MQ=255)

                                 ccGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAAcc                                                           >  2:4886170/1‑93 (MQ=255)

                                 ccGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAAcc                                                           <  2:5523278/93‑1 (MQ=255)

                                 ccGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAAcc                                                           <  1:4194321/93‑1 (MQ=255)

                                  cGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACct                                                          >  1:2274944/1‑93 (MQ=255)

                                  cGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACct                                                          <  2:5139152/93‑1 (MQ=255)

                                   ggctggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACctt                                                         >  1:4354048/1‑93 (MQ=255)

                                    gctggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACcttc                                                        <  2:4180352/93‑1 (MQ=255)

                                    gctggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACcttc                                                        <  2:6906309/93‑1 (MQ=255)

                                     ctggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACcttct                                                       <  1:201726/93‑1 (MQ=255)

                                     ctggctTGTGCAATAACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACcttct                                                       <  1:6893779/93‑1 (MQ=255)

                                      tggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACcttctt                                                      >  2:693659/1‑93 (MQ=255)

                                      tggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACcttctt                                                      >  2:1181654/1‑93 (MQ=255)

                                       ggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCttt                                                     <  2:3314426/93‑1 (MQ=255)

                                       ggctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCttt                                                     <  2:3758167/93‑1 (MQ=255)

                                        gctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCtttt                                                    >  2:478618/1‑93 (MQ=255)

                                        gctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCtttt                                                    <  2:3890126/93‑1 (MQ=255)

                                         ctTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCttttt                                                   <  1:5032057/93‑1 (MQ=255)

                                          ttGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTc                                                  >  2:6292013/1‑93 (MQ=255)

                                          ttGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTc                                                  >  1:1093289/1‑93 (MQ=255)

                                          ttGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTc                                                  >  2:4268131/1‑93 (MQ=255)

                                          ttGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTc                                                  <  1:4024078/93‑1 (MQ=255)

                                           tgtgCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCt                                                 <  1:1425222/93‑1 (MQ=255)

                                                  tGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGctc                                          >  1:2298188/1‑93 (MQ=255)

                                                     cTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTa                                       <  2:1281515/93‑1 (MQ=255)

                                                     cTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTa                                       <  2:4738100/93‑1 (MQ=255)

                                                       tCGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACg                                     <  2:5332262/93‑1 (MQ=255)

                                                        cGTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACgg                                    >  1:5276195/1‑93 (MQ=255)

                                                         gTTCTCTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGgc                                   >  1:1270916/1‑93 (MQ=255)

                                                         gTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGgc                                   >  1:3547749/1‑93 (MQ=255)

                                                           tctcTTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGgcgc                                 >  2:5607111/1‑93 (MQ=255)

                                                           tctcTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGgcgc                                 >  2:2617771/1‑93 (MQ=255)

                                                            ctctTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCa                                >  1:3154894/1‑93 (MQ=255)

                                                            ctctTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCa                                >  1:3411911/1‑93 (MQ=255)

                                                            ctctTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCa                                <  2:5618086/93‑1 (MQ=255)

                                                             tctTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCAt                               >  2:5224271/1‑93 (MQ=255)

                                                             tctTCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCAt                               <  2:3521261/93‑1 (MQ=255)

                                                             tctTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCAt                               >  2:6254319/1‑93 (MQ=255)

                                                               ttCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCAttt                             >  1:6327335/1‑93 (MQ=255)

                                                               ttCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCAttt                             <  1:3113578/93‑1 (MQ=255)

                                                               ttCATTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCAttt                             <  1:6033129/93‑1 (MQ=255)

                                                                  aTTGACAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGa                          <  1:912028/93‑1 (MQ=255)

                                                                  aTTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGa                          >  2:4007197/1‑93 (MQ=255)

                                                                       cAATATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTtgt                     >  1:5802335/1‑93 (MQ=255)

                                                                        aaTATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTtgtg                    <  1:751102/93‑1 (MQ=255)

                                                                        aaTATCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTtgtg                    <  2:7261206/93‑1 (MQ=255)

                                                                         atatCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTgg                   <  1:5338457/93‑1 (MQ=255)

                                                                          tatCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTggg                  >  1:3482195/1‑93 (MQ=255)

                                                                            tCCACCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTa                <  2:1653152/93‑1 (MQ=255)

                                                                              cACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTaaa              <  1:940670/93‑1 (MQ=255)

                                                                               aCCCATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAAt             >  2:3133972/1‑93 (MQ=255)

                                                                                 ccATTCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTa           >  1:6329158/1‑93 (MQ=255)

                                                                                  cATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTAc          >  1:365645/1‑93 (MQ=255)

                                                                                   aTTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTTGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTatg         >  2:6505028/1‑91 (MQ=255)

                                                                                   aTTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACg         >  1:2848028/1‑93 (MQ=255)

                                                                                    ttCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACgg        >  2:1396995/1‑93 (MQ=255)

                                                                                    ttCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACgg        >  1:7190621/1‑93 (MQ=255)

                                                                                     tCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGa       <  2:6117272/93‑1 (MQ=255)

                                                                                     tCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGa       <  1:6161159/93‑1 (MQ=255)

                                                                                     tCAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGa       <  2:2758858/93‑1 (MQ=255)

                                                                                     tCAGTACTTGACAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGa       <  1:4865905/93‑1 (MQ=255)

                                                                                      cAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACggga      >  1:2429897/1‑91 (MQ=255)

                                                                                      cAGTGCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGAt      <  1:5413612/93‑1 (MQ=255)

                                                                                      cAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAATTTACGGAt      <  2:4860458/93‑1 (MQ=255)

                                                                                       aGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGAtt     >  1:1013584/1‑93 (MQ=255)

                                                                                          aCTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGATTAAt  >  1:1449745/1‑93 (MQ=255)

ACGATGACGCAGACACTGTGCCCGCATTTGTTCCCGGCTGGCTTGTGCAATGACTTCGTTCTCTTCATTGACAATATCCACCCATTCAGTACTTGCCAAACGACGCTGTTCCATTATCAGGAAACCTTCTTTTTCTGGCGCTCTTACGGCGCATTTTGAGTTGTGGGTAAATTACGGATTAAT  >  NC_000913/2417083‑2417265

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