BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	115,968	+G	coding (1132/1386 nt)	hofB ←	HofB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	115,964	1	.	G	93.5%	374.7 / 18.8	124	coding (1136/1386 nt)	hofB	HofB
Reads supporting (aligned to +/- strand): ref base . (6/2); new base G (52/64); total (58/66)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 142 total aligned reads displayed.

GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑G‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051

 aCGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g                                                                                           <  1:1598656/93‑1 (MQ=255)

 aCGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g                                                                                           <  1:1306256/93‑1 (MQ=255)

  cGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g                                                                                         <  2:2478749/93‑1 (MQ=255)

  cGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g                                                                                         <  2:2478750/93‑1 (MQ=255)

  cGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g                                                                                         >  1:1399776/1‑93 (MQ=255)

   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g                                                                                       <  1:3365656/93‑1 (MQ=255)

    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g‑g                                                                                     <  1:3807995/93‑1 (MQ=255)

    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g‑g                                                                                     <  2:4299052/93‑1 (MQ=255)

     aTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g‑gg                                                                                    <  1:2163521/93‑2 (MQ=255)

       gACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g‑ggca                                                                                  <  2:3830096/93‑4 (MQ=255)

        aCCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g‑ggcag                                                                                 <  2:3851832/93‑5 (MQ=255)

        aCCGCCGTTATCGCCACAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g‑g‑g‑ggcag                                                                                 <  2:3851831/93‑5 (MQ=255)

          cGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑G‑G‑G‑GCAgcgg                                                                               <  1:4106227/93‑1 (MQ=255)

           ggCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAgcgg                                                                               >  2:1680534/1‑93 (MQ=255)

            gCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAACCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAgcggc                                                                              <  1:2499839/93‑1 (MQ=255)

            gCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAgcggc                                                                              <  2:1948744/93‑1 (MQ=255)

            gCGTTATGGGCAGAACTTCAAATAAGGCCGCACGACCATAAAAACCCTGGTAGCAATGTACACAACCGGGTCCCTTCCAGTGG‑G‑G‑GCAgcggg                                                                              >  2:3809243/1‑92 (MQ=255)

              gTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGa                                                                            >  2:777836/1‑93 (MQ=255)

              gTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGa                                                                            >  2:1334138/1‑93 (MQ=255)

               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACGCAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGAt                                                                           >  2:2157659/1‑93 (MQ=255)

               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGGGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGag                                                                           >  1:3876703/1‑92 (MQ=255)

               gtATGGGCGGATCTTCAAATAAGGCCGTACGACCATAAACGCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGAt                                                                           <  2:1908300/92‑1 (MQ=255)

                   gggCAGAACTTCAAATAAGGCCGTCCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGcc                                                                       <  2:3395965/93‑1 (MQ=255)

                   gggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGcc                                                                       >  2:396783/1‑93 (MQ=255)

                    ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCa                                                                      <  2:3149105/93‑1 (MQ=255)

                      cAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATa                                                                    <  2:2424875/93‑1 (MQ=255)

                       aGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATAc                                                                   >  2:1702399/1‑93 (MQ=255)

                           cTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTg                                                               >  1:4657860/1‑93 (MQ=255)

                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGt                                                              >  1:939013/1‑93 (MQ=255)

                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGt                                                              >  1:1069874/1‑93 (MQ=255)

                             tCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTc                                                             <  1:3107847/93‑1 (MQ=255)

                             tCAAATAAGGCCGTACGACCAGAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTc                                                             <  1:3801455/93‑1 (MQ=255)

                              cAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCt                                                            <  2:2790172/93‑1 (MQ=255)

                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTg                                                           <  2:10126/93‑1 (MQ=255)

                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTg                                                           >  1:1972790/1‑93 (MQ=255)

                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTgg                                                          >  1:4716593/1‑93 (MQ=255)

                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑G‑GTG‑GCAGCGGCGTTGGCCATACATTGTCTgg                                                          >  1:18197/1‑93 (MQ=255)

                                 aTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGa                                                         <  1:596899/93‑1 (MQ=255)

                                 aTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGa                                                         >  2:1559425/1‑93 (MQ=255)

                                  tAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGaa                                                        <  1:1668767/93‑1 (MQ=255)

                                   aaGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAAt                                                       <  1:4179646/93‑1 (MQ=255)

                                     ggCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAAtgt                                                     >  1:2838120/1‑93 (MQ=255)

                                     ggCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAAtgt                                                     <  1:1205226/93‑1 (MQ=255)

                                      gCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAAtgtg                                                    <  1:2925503/93‑1 (MQ=255)

                                       cagcaCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTgg                                                   <  2:325106/89‑1 (MQ=255)

                                        cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTggt                                                  >  2:2853832/1‑92 (MQ=255)

                                        cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGa                                                  <  1:1147999/93‑1 (MQ=255)

                                         gTACGACCATAAAAGCCGGGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  2:1284029/93‑1 (MQ=255)

                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  2:547942/93‑1 (MQ=255)

                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                >  2:2078316/1‑93 (MQ=255)

                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                <  2:1212469/93‑1 (MQ=255)

                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                <  1:4050621/93‑1 (MQ=255)

                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                <  2:2622022/93‑1 (MQ=255)

                                           aCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATgg                                               <  1:1174156/93‑1 (MQ=255)

                                           aCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATgg                                               <  1:3549644/93‑1 (MQ=255)

                                           aCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATgg                                               <  1:3300625/93‑1 (MQ=255)

                                             gACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGc                                             <  2:2397414/93‑1 (MQ=255)

                                             gACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGc                                             >  2:4118256/1‑93 (MQ=255)

                                                  tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                        >  1:895537/1‑93 (MQ=255)

                                                  tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                        >  2:3744626/1‑93 (MQ=255)

                                                    aaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCtt                                      >  2:303529/1‑93 (MQ=255)

                                                     aaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtg                                     >  1:3830000/1‑93 (MQ=255)

                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  1:3118435/93‑1 (MQ=255)

                                                       aCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgct                                   >  2:4588013/1‑93 (MQ=255)

                                                         cGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgctgc                                 >  2:4120541/1‑93 (MQ=255)

                                                          gTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGcc                                >  1:121126/1‑93 (MQ=255)

                                                             gTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGc                             >  2:4297906/1‑93 (MQ=255)

                                                             gTAGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGc                             >  2:2387419/1‑93 (MQ=255)

                                                               aGCAATGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGa                           <  2:1622338/93‑1 (MQ=255)

                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  2:618707/93‑1 (MQ=255)

                                                                   aTGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGTGACAAt                       <  2:1168259/93‑1 (MQ=255)

                                                                   aTGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAt                       >  2:3692576/1‑93 (MQ=255)

                                                                   aTGTACACAACCGGGTGCCTGCCAGT‑GTG‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAt                       >  1:3731505/1‑93 (MQ=255)

                                                                   aTGTACACAACCGGGTGCCTGCCAGT‑G‑G‑GTGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGAAAAt                       >  2:1035444/1‑93 (MQ=255)

                                                                    tGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      >  2:2719939/1‑93 (MQ=255)

                                                                    tGTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      <  1:4309936/93‑1 (MQ=255)

                                                                     gTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgt                     <  1:4244767/93‑1 (MQ=255)

                                                                     gTACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGAGACAAtgt                     <  2:1152145/93‑1 (MQ=255)

                                                                      tACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgtg                    <  1:3011191/93‑1 (MQ=255)

                                                                      tACACAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgtg                    >  1:1806846/1‑93 (MQ=255)

                                                                       acacAACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTgg                   <  1:1689264/93‑1 (MQ=255)

                                                                        cacaACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTggg                  >  2:3533551/1‑93 (MQ=255)

                                                                        cacaACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTggg                  <  1:771676/93‑1 (MQ=255)

                                                                        cacaACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTggg                  >  1:467749/1‑93 (MQ=255)

                                                                          caACCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCa                >  1:1710875/1‑93 (MQ=255)

                                                                            aCCGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCaaa              <  2:3084588/93‑1 (MQ=255)

                                                                             ccGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAg             >  1:1650155/1‑93 (MQ=255)

                                                                             ccGGGTGCCTGCCAGT‑G‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGt            >  2:4497872/1‑93 (MQ=255)

                                                                              cGGGTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGt            <  2:2100391/93‑1 (MQ=255)

                                                                               gggTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           >  2:1673223/1‑93 (MQ=255)

                                                                               gggTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           <  1:2036769/93‑1 (MQ=255)

                                                                                ggTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGttt          <  2:4544987/93‑1 (MQ=255)

                                                                                 gTGCCTGCCAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTg         <  1:854029/93‑1 (MQ=255)

                                                                                  tgcctgccAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGc        <  2:4451199/93‑1 (MQ=255)

                                                                                   gcctgccAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCg       >  1:1389540/1‑93 (MQ=255)

                                                                                   gcctgccAGTGG‑G‑G‑GCAGCGGCGATGGCCACACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCg       <  2:4637760/93‑1 (MQ=255)

                                                                                      tgccAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    <  2:4610501/93‑1 (MQ=255)

                                                                                      tgccAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    <  1:716119/93‑1 (MQ=255)

                                                                                       gccAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAcc   <  1:1856211/93‑1 (MQ=255)

                                                                                        ccAGTGG‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCa  <  1:1825257/93‑1 (MQ=255)

GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑G‑G‑G‑GCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051

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