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	98.6%	660.8 / 1.5	209	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (121/85); total (123/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 233 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGt                                                                                          <  1:6090170/93‑1 (MQ=255)

ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGt                                                                                          <  2:1097265/93‑1 (MQ=255)

   gcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccc                                                                                      >  1:2229395/1‑93 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                     >  2:603873/1‑92 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                     <  2:5856681/93‑2 (MQ=255)

     gcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccca                                                                                    >  1:5837088/1‑91 (MQ=255)

      ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaa                                                                                   <  1:3009666/93‑4 (MQ=255)

       cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                  >  1:4125987/1‑89 (MQ=255)

         ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCCAaccc                                                                                <  1:3623994/93‑1 (MQ=255)

         ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCCAaccc                                                                                >  1:4936111/1‑93 (MQ=255)

         ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCCAaccc                                                                                >  2:5698688/1‑93 (MQ=255)

            ctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCAc                                                                              <  1:3608523/93‑1 (MQ=255)

             tgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACg                                                                             <  2:1518409/93‑1 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGc                                                                            >  2:1117417/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGcc                                                                           >  1:4580952/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGcc                                                                           >  2:160707/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCaa                                                                         <  2:1550896/93‑1 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCaa                                                                         >  1:3127188/1‑93 (MQ=255)

                  tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAAc                                                                        <  1:3713869/93‑1 (MQ=255)

                    ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCa                                                                      >  1:1795623/1‑93 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGc                                                                    <  1:2467660/93‑1 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGc                                                                    <  1:3365690/93‑1 (MQ=255)

                         cAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGcaca                                                                 >  2:466010/1‑93 (MQ=255)

                         cAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGcaca                                                                 >  2:5446395/1‑93 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACaa                                                                >  2:5941988/1‑93 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACaa                                                                <  1:714517/93‑1 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAata                                                              <  2:795312/93‑1 (MQ=255)

                             tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAatat                                                             <  2:5310474/93‑1 (MQ=255)

                             tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAatat                                                             >  1:2322579/1‑93 (MQ=255)

                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATg                                                            >  1:5804883/1‑93 (MQ=255)

                               aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGa                                                           <  1:4830845/93‑1 (MQ=255)

                                aGTAAGCCGGGCCTGCCCACTTCAACAATTCAAACGCAATCACTGAGCCGGAAAATAGCGTCCCCAACCCCACGCCACCCAGCACAATATGaa                                                          <  1:3463978/93‑1 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGaa                                                          <  1:1218787/93‑1 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGaa                                                          <  2:2663275/93‑1 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         >  2:1561341/1‑93 (MQ=255)

                                  tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATc                                                        <  2:4778785/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  1:795882/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     >  1:2711024/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  1:1802694/93‑1 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    >  1:3424449/1‑93 (MQ=255)

                                       cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAg                                                   >  2:4042146/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  2:4599788/93‑1 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  2:3831685/93‑1 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                <  2:4885833/93‑1 (MQ=255)

                                           ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                                               <  1:6222220/93‑1 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:4006076/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:4846759/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:2231650/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:5085574/1‑93 (MQ=255)

                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  2:23914/1‑93 (MQ=255)

                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            >  1:902688/1‑93 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  1:6163770/1‑93 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  1:3795264/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:5046223/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:5504425/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  2:4148717/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:4405322/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  1:2160179/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  1:5472596/1‑93 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGTCAGTCCGGTCTGAAGcc                                     <  1:4498817/93‑1 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  2:464024/93‑1 (MQ=255)

                                                        acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  >  1:5909138/1‑93 (MQ=255)

                                                        acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCAAGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  >  1:4257200/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAATCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAACa                                >  2:2036045/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:4448601/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               <  2:3537750/93‑1 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:3231709/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:2964284/1‑93 (MQ=255)

                                                              caaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAg                            >  1:1147538/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  2:691524/93‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  2:2824750/93‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:5221963/93‑1 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:3494475/1‑93 (MQ=255)

                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        <  2:2087279/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  1:1670701/1‑93 (MQ=255)

                                                                    cAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa                      >  1:2974199/1‑93 (MQ=255)

                                                                    cAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa                      <  1:1013363/93‑1 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     <  1:3668318/93‑1 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:1282579/1‑93 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:2731456/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  1:2038391/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  1:1961628/1‑93 (MQ=255)

                                                                        cACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGTCAccgc                  <  2:5242158/93‑1 (MQ=255)

                                                                          cTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGTCAccgcgc                <  2:6149457/93‑1 (MQ=255)

                                                                          cTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                <  2:3658048/93‑1 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               >  1:2758338/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               >  1:1568287/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               <  2:5502212/93‑1 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             >  2:1060119/1‑93 (MQ=255)

                                                                               cGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGa           >  2:3120994/1‑93 (MQ=255)

                                                                               cGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGa           >  1:5981351/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  1:824841/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  1:4325864/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  1:5104931/1‑93 (MQ=255)

                                                                                   aaaaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       >  1:5266389/1‑93 (MQ=255)

                                                                                    aaaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAcc      >  1:4466950/1‑93 (MQ=255)

                                                                                      aTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  1:1606922/93‑1 (MQ=255)

                                                                                      aTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  2:1410040/93‑1 (MQ=255)

                                                                                      aTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  2:4022724/93‑1 (MQ=255)

                                                                                       tAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   >  1:1401731/1‑93 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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