BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,614	(A)_7→6	coding (569/726 nt)	yjjQ →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,608	0	A	.	99.8%	2059.0 / ‑3.0	451	coding (563/726 nt)	yjjQ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (240/210); total (241/210)
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.97e-01

Only 100 of 530 total aligned reads displayed.

AGAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTT  >  NC_000913/4601516‑4601696

agagCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCa                                                                                          <  1:1850380/93‑1 (MQ=255)

agagCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGAGCCCa                                                                                          >  2:6420511/1‑93 (MQ=255)

     gcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaa                                                                                     >  2:9208998/1‑93 (MQ=255)

       aaGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaca                                                                                   <  1:1413650/93‑3 (MQ=255)

       aaGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaca                                                                                   <  1:4450364/93‑3 (MQ=255)

       aaGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaca                                                                                   <  2:4852114/93‑3 (MQ=255)

        aGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacag                                                                                  >  1:4871177/1‑90 (MQ=255)

         ggTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacagt                                                                                 >  1:5809728/1‑89 (MQ=255)

         ggTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaacagt                                                                                 >  2:4472174/1‑89 (MQ=255)

          gTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTg                                                                               >  2:3014933/1‑93 (MQ=255)

          gTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTg                                                                               <  2:2536957/93‑1 (MQ=255)

          gTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTg                                                                               <  1:2828409/93‑1 (MQ=255)

          gTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTg                                                                               <  2:1871225/93‑1 (MQ=255)

          gTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTg                                                                               >  1:4259553/1‑93 (MQ=255)

           tACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGc                                                                              <  1:5657189/93‑1 (MQ=255)

            aCTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCg                                                                             >  1:919869/1‑93 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  1:6506178/93‑1 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  2:7599445/93‑1 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  1:5988324/93‑1 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  1:3778605/93‑1 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  1:9339829/93‑1 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  1:336612/93‑1 (MQ=255)

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                            <  2:5473705/93‑1 (MQ=255)

              tGCGGCTATTAGGTAAAGGATGGGGCATCAGCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           <  2:7058987/93‑1 (MQ=255)

              tGCGGCTATTAGATAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                           >  2:7456104/1‑93 (MQ=255)

                cGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGc                                                                         >  2:4253656/1‑93 (MQ=255)

                 ggCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCg                                                                        >  2:5165638/1‑93 (MQ=255)

                 ggCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCg                                                                        <  2:8210029/93‑1 (MQ=255)

                  gCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGt                                                                       <  1:6812325/93‑1 (MQ=255)

                     aTTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGa                                                                    <  1:6397946/93‑1 (MQ=255)

                     aTTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGa                                                                    >  2:1616756/1‑93 (MQ=255)

                     aTTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGa                                                                    >  1:8558159/1‑93 (MQ=255)

                      ttAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGAc                                                                   <  2:4768606/93‑1 (MQ=255)

                       tAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACt                                                                  <  2:2478865/93‑1 (MQ=255)

                          gTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGc                                                               >  2:9293695/1‑93 (MQ=255)

                            aaaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCGC‑AAAAACCAGGGCGATGCGTACACTGGCaa                                                             >  2:1007336/1‑93 (MQ=255)

                            aaaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGGaa                                                             >  1:3886369/1‑93 (MQ=255)

                             aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAt                                                            <  2:8089077/93‑1 (MQ=255)

                             aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAAt                                                            <  2:1246023/93‑1 (MQ=255)

                               ggATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCTTCGACTGGCAATTc                                                          <  2:3579673/93‑1 (MQ=255)

                                gATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTca                                                         >  2:5647875/1‑93 (MQ=255)

                                 aTGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTcac                                                        >  2:6033722/1‑93 (MQ=255)

                                  tGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTcaca                                                       >  2:9696171/1‑93 (MQ=255)

                                   ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                      >  2:8800423/1‑93 (MQ=255)

                                   ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                      >  1:5501612/1‑93 (MQ=255)

                                    gggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGc                                                     >  2:4047402/1‑93 (MQ=255)

                                    gggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGc                                                     >  2:2321792/1‑93 (MQ=255)

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                >  1:5889796/1‑93 (MQ=255)

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                >  2:9097137/1‑93 (MQ=255)

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                >  2:5342339/1‑93 (MQ=255)

                                         tCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGc                                                <  2:198820/93‑1 (MQ=255)

                                          cAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCt                                               <  1:520677/93‑1 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAACACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             <  1:8831875/93‑1 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             <  1:9698508/93‑1 (MQ=255)

                                            aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                                             >  1:9487173/1‑93 (MQ=255)

                                             ccAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGaa                                            >  2:5406008/1‑93 (MQ=255)

                                              cAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGaaa                                           <  1:6584934/93‑1 (MQ=255)

                                               aGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAAt                                          >  1:7926813/1‑93 (MQ=255)

                                               aGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAAt                                          >  2:7752819/1‑93 (MQ=255)

                                                gATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATg                                         >  2:6629116/1‑93 (MQ=255)

                                                gATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATg                                         >  1:3259427/1‑93 (MQ=255)

                                                 aTAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGt                                        <  2:3304497/93‑1 (MQ=255)

                                                    gCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATg                                     <  2:3757109/93‑1 (MQ=255)

                                                    gCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATg                                     >  2:6258044/1‑93 (MQ=255)

                                                      ttCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCa                                   >  1:3620803/1‑93 (MQ=255)

                                                      atCATTGCTTAAGAAAAGTCGTCACACTGTCAGCGCCCAAAAGAA‑AGAACGCTGCAACGAAAGGCCATTCACAACAAAGCAGAAATAAATGCa                                   <  1:7723344/92‑1 (MQ=255)

                                                       tCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCAt                                  <  2:7332640/93‑1 (MQ=255)

                                                       tCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCAt                                  <  2:9742633/93‑1 (MQ=255)

                                                       tCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCAt                                  >  1:837329/1‑93 (MQ=255)

                                                        cATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATg                                 >  2:4660831/1‑93 (MQ=255)

                                                         aTTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATgg                                <  2:433787/93‑1 (MQ=255)

                                                          ttGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGa                               >  1:1003631/1‑93 (MQ=255)

                                                          ttGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGa                               >  1:2604644/1‑93 (MQ=255)

                                                          ttGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGa                               >  1:817757/1‑93 (MQ=255)

                                                           tGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGGAACGCTGAAATGTATGCATGGAt                              >  2:1439774/1‑93 (MQ=255)

                                                            gCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATa                             >  1:2346794/1‑93 (MQ=255)

                                                             cTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaa                            <  1:2995498/93‑1 (MQ=255)

                                                             cTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaa                            >  1:6130251/1‑93 (MQ=255)

                                                              ttAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa                           >  2:5934376/1‑93 (MQ=255)

                                                              ttAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa                           >  1:2123478/1‑93 (MQ=255)

                                                              ttAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa                           <  2:1918396/93‑1 (MQ=255)

                                                               tAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAAt                          >  1:5077225/1‑93 (MQ=255)

                                                               tAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAAt                          <  1:7309271/93‑1 (MQ=255)

                                                                aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         <  2:5072215/93‑1 (MQ=255)

                                                                aaGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATa                         <  1:1755869/93‑1 (MQ=255)

                                                                  gAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgc                       <  2:7903011/93‑1 (MQ=255)

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  1:766298/1‑93 (MQ=255)

                                                                    aaaGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc                     >  2:843872/1‑93 (MQ=255)

                                                                      aGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcgc                   <  1:7821566/93‑1 (MQ=255)

                                                                       gTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCa                  <  2:6153344/93‑1 (MQ=255)

                                                                          ataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAggg               <  2:6825305/93‑1 (MQ=255)

                                                                           taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgggg              >  1:987800/1‑92 (MQ=255)

                                                                            aaaaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGCCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTg             <  2:2894062/93‑1 (MQ=255)

                                                                                 tATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAaga        <  2:2755337/93‑1 (MQ=255)

                                                                                 tATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAaga        >  2:9469260/1‑93 (MQ=255)

                                                                                 tATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAaga        >  1:1587138/1‑93 (MQ=255)

                                                                                 tATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAaga        >  2:1684393/1‑93 (MQ=255)

                                                                                  aTCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAagag       >  1:9375938/1‑93 (MQ=255)

                                                                                   tCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAgaga      >  2:6746787/1‑93 (MQ=255)

                                                                                     aGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaac    <  1:7080564/93‑1 (MQ=255)

AGAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTT  >  NC_000913/4601516‑4601696

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