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.7%	1560.6 / ‑4.2	342	coding (563/726 nt)	yjjQ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (178/163); total (178/164)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 376 total aligned reads displayed.

AGAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTA  >  NC_000913/4601516‑4601697

agagCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAAGATAACTATCAGCGCCCa                                                                                           >  1:7042009/1‑93 (MQ=255)

 gagCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaa                                                                                          >  1:519775/1‑93 (MQ=255)

   gcgcAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaa                                                                                        >  1:2137061/1‑93 (MQ=255)

       aaGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaca                                                                                    >  2:2211901/1‑91 (MQ=255)

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

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

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

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

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

             cTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGa                                                                             >  1:6102479/1‑93 (MQ=255)

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

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

              tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                            >  1:7121781/1‑93 (MQ=255)

              tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                            >  2:804863/1‑93 (MQ=255)

              tGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGAt                                                                            >  1:3856887/1‑93 (MQ=255)

               gCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATg                                                                           <  1:6592388/93‑1 (MQ=255)

               gCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATg                                                                           <  2:903498/93‑1 (MQ=255)

               gCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATg                                                                           <  1:7237348/93‑1 (MQ=255)

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

                   cTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTc                                                                       >  1:5655926/1‑93 (MQ=255)

                   cTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTc                                                                       >  1:187352/1‑93 (MQ=255)

                   cTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTc                                                                       >  1:4707314/1‑93 (MQ=255)

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

                     aTTAGGTAAAGGAGGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGa                                                                     <  2:5856697/93‑1 (MQ=255)

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

                      ttAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGAc                                                                    >  2:6394026/1‑93 (MQ=255)

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

                         ggTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTgg                                                                 <  2:3119115/93‑1 (MQ=255)

                            aaaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCaa                                                              >  2:3382936/1‑93 (MQ=255)

                            aaaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCaa                                                              >  1:2372588/1‑93 (MQ=255)

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

                                gATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTca                                                          <  2:6961634/93‑1 (MQ=255)

                                  tGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTcaca                                                        <  1:3142346/93‑1 (MQ=255)

                                   ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                       <  1:1047359/93‑1 (MQ=255)

                                   ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                       <  2:3253581/93‑1 (MQ=255)

                                   ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                       <  2:7530126/93‑1 (MQ=255)

                                   ggggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAg                                                       <  1:1553019/93‑1 (MQ=255)

                                    gggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGc                                                      <  2:1355997/93‑1 (MQ=255)

                                      gCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCaa                                                    <  2:3159808/93‑1 (MQ=255)

                                        aTCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACg                                                  >  1:6909967/1‑93 (MQ=255)

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

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

                                          cAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCt                                                >  2:3741380/1‑93 (MQ=255)

                                          cAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCAC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCt                                                <  1:7902762/93‑1 (MQ=255)

                                           aaCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTg                                               >  2:7811582/1‑93 (MQ=255)

                                              cAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGaaa                                            >  2:6064664/1‑93 (MQ=255)

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

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

                                                 aTAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCCCAGCAACGCTGAAATGt                                         >  1:2088718/1‑93 (MQ=255)

                                                 aTAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGt                                         >  2:5162305/1‑93 (MQ=255)

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

                                                     cTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGc                                     <  1:5955522/93‑1 (MQ=255)

                                                     cTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGc                                     <  2:6306821/93‑1 (MQ=255)

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

                                                      ttCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCa                                    <  2:2376946/93‑1 (MQ=255)

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

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

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

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

                                                        taTTGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATg                                  <  1:4229444/92‑1 (MQ=255)

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

                                                          ttGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGa                                <  2:2604554/93‑1 (MQ=255)

                                                           tGCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGAt                               >  1:5022565/1‑93 (MQ=255)

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

                                                            gCTTAAGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATa                              <  1:894552/93‑1 (MQ=255)

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

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

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

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

                                                               tAAGAAAAGTAAAAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAAt                           <  2:7158672/93‑1 (MQ=255)

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

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

                                                                 aGAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAg                         >  1:4669619/1‑93 (MQ=255)

                                                                  gAAAAGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgc                        >  1:1714146/1‑93 (MQ=255)

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

                                                                  gAAAAGGAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgc                        <  1:3059944/93‑1 (MQ=255)

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

                                                                      aGTAATAAAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcgc                    >  2:7463050/1‑93 (MQ=255)

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

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

                                                                        taatatCACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCTATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAg                  <  2:7042008/93‑1 (MQ=255)

                                                                        taataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAg                  <  1:6201618/93‑1 (MQ=255)

                                                                         aataaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAgg                 >  1:5518347/1‑93 (MQ=255)

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

                                                                           taaAACTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGt               <  1:5120321/93‑1 (MQ=255)

                                                                              aaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCa            >  2:4261247/1‑93 (MQ=255)

                                                                              aaCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCa            >  1:5320648/1‑93 (MQ=255)

                                                                               aCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCaa           >  1:3769321/1‑93 (MQ=255)

                                                                               aCTATCAGCGCCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCaa           >  2:2787502/1‑93 (MQ=255)

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

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

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

                                                                                      gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCCAGAGaact    >  2:3409522/1‑93 (MQ=255)

                                                                                      gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact    <  2:1024287/93‑1 (MQ=255)

                                                                                      gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact    >  2:3075817/1‑93 (MQ=255)

                                                                                      gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact    <  2:4047862/93‑1 (MQ=255)

                                                                                      gcgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaact    <  1:1383096/93‑1 (MQ=255)

                                                                                       cgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaactt   >  1:1253128/1‑93 (MQ=255)

                                                                                       cgcCC‑AAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaactt   >  2:3860119/1‑93 (MQ=255)

                                                                                         ccccAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGaactta  <  2:2538623/93‑1 (MQ=255)

AGAGCGCAAGGTACTGCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTA  >  NC_000913/4601516‑4601697

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