BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	96.6%	618.1 / 13.9	204	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/4); new base C (123/74); total (126/78)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 217 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

 gcgcgcCACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                         >  1:95659/1‑93 (MQ=255)

    cgcCACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑ccc                                                                                     <  1:2979590/93‑2 (MQ=255)

    cgcCACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑ccc                                                                                     >  2:2217605/1‑92 (MQ=255)

     gcCACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑ccca                                                                                    <  2:2726216/93‑3 (MQ=255)

      ccACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑cccaa                                                                                   >  2:4515322/1‑90 (MQ=255)

       cACTGCTGGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c‑cccaac                                                                                  >  2:602966/1‑89 (MQ=255)

         ctgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑CCCAaccc                                                                                <  2:1116935/93‑1 (MQ=255)

         ctgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑CCCAaccc                                                                                >  1:3281494/1‑93 (MQ=255)

          tgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAAcccc                                                                                >  1:3123142/1‑93 (MQ=255)

          tgctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAAcccc                                                                                <  1:2610419/93‑1 (MQ=255)

           gctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCa                                                                               >  2:3945825/1‑93 (MQ=255)

            ctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCAc                                                                              >  1:4398980/1‑93 (MQ=255)

            ctgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCAc                                                                              <  2:2409627/93‑1 (MQ=255)

             tgGATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACg                                                                             >  1:3664733/1‑93 (MQ=255)

              ggATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGc                                                                            <  1:2209185/93‑1 (MQ=255)

              ggATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGc                                                                            <  2:1693955/93‑1 (MQ=255)

              ggATTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGc                                                                            >  2:3367005/1‑93 (MQ=255)

                aTTTCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCa                                                                          <  1:4172719/93‑1 (MQ=255)

                 tttccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCaa                                                                         <  1:930833/90‑1 (MQ=255)

                    ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCa                                                                      >  1:3269869/1‑93 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGca                                                                   >  1:2435188/1‑93 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACaa                                                                <  1:3454943/93‑1 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACaa                                                                >  2:1203520/1‑93 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAat                                                               <  1:1488548/93‑1 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAat                                                               <  2:95652/93‑1 (MQ=255)

                             tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAatat                                                             >  2:2207526/1‑93 (MQ=255)

                             tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAatat                                                             <  1:1940572/93‑1 (MQ=255)

                               aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGa                                                           >  2:4360732/1‑93 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGaa                                                          >  1:3972437/1‑93 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGaa                                                          <  2:2868079/93‑1 (MQ=255)

                                aGTAAGCCGCCCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGaa                                                          <  1:848701/93‑1 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAAt                                                         >  1:2774154/1‑93 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAAt                                                         <  1:2349668/93‑1 (MQ=255)

                                  tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATc                                                        <  1:3883260/93‑1 (MQ=255)

                                  tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATc                                                        <  2:861228/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  2:3177569/1‑93 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      <  1:2766670/93‑1 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAGCGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCCCAATATGAATCGCCa                                                    <  2:620974/93‑1 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:2517806/93‑1 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:2235183/93‑1 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:670668/1‑93 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              <  1:3398684/93‑1 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCACCCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           <  2:4668560/93‑1 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  1:4205572/1‑93 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  2:3570197/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  2:1074445/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  2:2028740/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  1:798267/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  1:573068/93‑1 (MQ=255)

                                                   cttcaacactccaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       <  2:3201618/93‑1 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  2:3150835/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  1:3346866/1‑93 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGT‑CTCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  2:1645066/1‑93 (MQ=255)

                                                      caacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                                    <  1:1295254/93‑1 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  2:4139969/93‑1 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  2:2677053/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  2:507802/1‑93 (MQ=255)

                                                        acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  >  1:4213107/1‑93 (MQ=255)

                                                         cacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                                 >  2:3048994/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:3758646/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:1356009/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:1448759/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  1:475895/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  1:1007337/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:1962759/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:4215970/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:1813019/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              <  2:280082/93‑1 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              <  1:3119975/93‑1 (MQ=255)

                                                              caaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCGCAATATGAAGCGCCAGTCCGGTCTGAAGCCCAGCAATAg                            <  1:1613837/93‑1 (MQ=255)

                                                              caaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg                            <  1:3563991/93‑1 (MQ=255)

                                                              caaACGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg                            >  1:2394219/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  2:3823727/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:2001100/93‑1 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:3635162/1‑93 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  1:4596595/93‑1 (MQ=255)

                                                                 aCGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACg                         >  2:752650/1‑93 (MQ=255)

                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        <  2:4291834/93‑1 (MQ=255)

                                                                    cAATCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa                      >  2:1237668/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:698495/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:1795416/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:3186396/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGACCAGCAATAGACGCCAc                     >  2:1432598/1‑93 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    <  2:3115327/93‑1 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:3368203/1‑93 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  2:2362437/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   <  2:2916787/93‑1 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGT‑CACCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  1:462568/1‑93 (MQ=255)

                                                                        cACTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  <  1:1580893/93‑1 (MQ=255)

                                                                         aCTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  1:4283912/1‑93 (MQ=255)

                                                                         aCTGAGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  2:1178860/1‑93 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             >  1:565983/1‑93 (MQ=255)

                                                                              gCGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  1:198543/1‑93 (MQ=255)

                                                                               cGGGAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGa           <  1:3712851/93‑1 (MQ=255)

                                                                                gggAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  1:4345231/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  1:1079310/1‑93 (MQ=255)

                                                                                   aaaaTAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       >  1:2640296/1‑93 (MQ=255)

                                                                                    aaaTAGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAcc      >  2:2565420/1‑93 (MQ=255)

                                                                                        aGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  >  1:4650737/1‑93 (MQ=255)

                                                                                        aGCGTCC‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  >  1:2903883/1‑93 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑C‑CCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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