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	98.2%	706.8 / 4.3	218	coding (1136/1386 nt)	hofB	HofB
Reads supporting (aligned to +/- strand): ref base . (3/1); new base G (109/105); total (112/106)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 255 total aligned reads displayed.

GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑GG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051

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

 aCGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑g                                                                                         >  1:4198874/1‑93 (MQ=255)

  cGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg                                                                                        <  1:6178417/93‑1 (MQ=255)

  cGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg                                                                                        <  1:1439335/93‑1 (MQ=255)

   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑g                                                                                      >  2:4732675/1‑93 (MQ=255)

   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑g                                                                                      >  2:364958/1‑93 (MQ=255)

   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑g                                                                                      <  2:6357237/93‑1 (MQ=255)

   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑g                                                                                      >  1:6114547/1‑93 (MQ=255)

   gAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑g                                                                                      <  2:2846227/93‑1 (MQ=255)

   gAATGACCGGCGTTATGGACAGAACTTCAAATAAGGCCGTACAACCATAAAACCCGTGGTAGCAATGTACCCACCCGGGTGCCTGCCAGT‑gg‑g                                                                                      <  1:6726289/93‑1 (MQ=255)

    aaTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑gg                                                                                     >  2:6616257/1‑93 (MQ=255)

     aTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATTAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑ggg                                                                                    >  1:5377392/1‑92 (MQ=255)

     aTGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gc‑ggg                                                                                    <  1:3524680/93‑2 (MQ=255)

      tGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑gggc                                                                                   <  2:1384060/93‑3 (MQ=255)

        aCCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑gg‑gggcag                                                                                 <  1:36650/93‑5 (MQ=255)

          cGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAgcg                                                                                >  2:3120082/1‑93 (MQ=255)

            gCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAgcggc                                                                              <  1:5005207/93‑1 (MQ=255)

              ggtATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGa                                                                            >  1:5177832/3‑93 (MQ=255)

              gTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGa                                                                            >  2:4902484/1‑93 (MQ=255)

              gTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGa                                                                            <  2:4961804/93‑1 (MQ=255)

               ttATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGAt                                                                           >  1:2267655/1‑93 (MQ=255)

                tATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATg                                                                          <  1:3381486/93‑1 (MQ=255)

                tATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATg                                                                          >  1:1026124/1‑93 (MQ=255)

                tATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑GGTGGCAGCGGCGATg                                                                          >  2:2599392/1‑93 (MQ=255)

                 aTGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATgg                                                                         >  1:5428218/1‑93 (MQ=255)

                   gggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGcc                                                                       >  2:1135720/1‑93 (MQ=255)

                    ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCa                                                                      <  1:722555/93‑1 (MQ=255)

                    ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCa                                                                      >  1:3547554/1‑93 (MQ=255)

                    ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCa                                                                      <  1:3862932/93‑1 (MQ=255)

                     gCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCAt                                                                     >  1:2466746/1‑93 (MQ=255)

                       aGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATAc                                                                   >  2:633293/1‑93 (MQ=255)

                        gAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACa                                                                  <  2:6391520/93‑1 (MQ=255)

                          aCTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACAtt                                                                >  2:1795168/1‑93 (MQ=255)

                           cTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTg                                                               <  2:2502545/93‑1 (MQ=255)

                           cTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTg                                                               >  1:3051633/1‑93 (MQ=255)

                            ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGt                                                              >  1:2762325/1‑93 (MQ=255)

                             tCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTc                                                             >  2:1353127/1‑93 (MQ=255)

                              cAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCt                                                            <  2:3945207/93‑1 (MQ=255)

                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTg                                                           <  1:129444/93‑1 (MQ=255)

                               aaaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTg                                                           >  2:879673/1‑93 (MQ=255)

                                aaTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTgg                                                          >  2:5636633/1‑93 (MQ=255)

                                 aTAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGa                                                         >  2:5484717/1‑93 (MQ=255)

                                  tAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGaa                                                        <  2:4119746/93‑1 (MQ=255)

                                  tAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGaa                                                        <  1:6107917/93‑1 (MQ=255)

                                  tAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGaa                                                        >  1:3745565/1‑93 (MQ=255)

                                   aaGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAAt                                                       <  2:5081102/93‑1 (MQ=255)

                                   aaGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAAt                                                       >  1:214168/1‑93 (MQ=255)

                                    aGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAAtg                                                      >  1:6438893/1‑93 (MQ=255)

                                     ggCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAAttg                                                     >  2:198616/1‑91 (MQ=255)

                                     ggCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAAtgt                                                     >  1:4938517/1‑93 (MQ=255)

                                      gCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAAtgtg                                                    <  1:398140/93‑1 (MQ=255)

                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  1:3278451/93‑1 (MQ=255)

                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 >  1:253500/1‑93 (MQ=255)

                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  2:1080895/93‑1 (MQ=255)

                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 <  2:6731489/93‑1 (MQ=255)

                                         gTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                                 >  2:4826975/1‑93 (MQ=255)

                                          tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATg                                                >  1:5677625/1‑93 (MQ=255)

                                           aCGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATgg                                               <  1:1351612/93‑1 (MQ=255)

                                             gACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATggcc                                             >  2:1821426/1‑91 (MQ=255)

                                             gACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGc                                             >  1:4317777/1‑93 (MQ=255)

                                                cATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcc                                          <  1:1214254/93‑1 (MQ=255)

                                                 atgaaaaCCGTGGTACCAATGTACTCAACCGGTTGCCTGCCAGTGGG‑GGCAGCGGCGATGTCCATACATTGTCTGGAATGTGGATGGGCTccc                                         <  1:6287032/90‑1 (MQ=255)

                                                  tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                        <  2:4607584/93‑1 (MQ=255)

                                                   aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCt                                       <  1:4653943/93‑1 (MQ=255)

                                                   aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCt                                       <  1:4550713/93‑1 (MQ=255)

                                                   aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCt                                       >  2:141436/1‑93 (MQ=255)

                                                   aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGATCCCCt                                       <  2:6083045/93‑1 (MQ=255)

                                                    aaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCtt                                      >  2:6129973/1‑93 (MQ=255)

                                                    aaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCtt                                      <  1:1820127/93‑1 (MQ=255)

                                                    aaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCtt                                      <  2:5733375/93‑1 (MQ=255)

                                                    aaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCtt                                      <  1:1874032/93‑1 (MQ=255)

                                                     aaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtg                                     >  1:1971253/1‑93 (MQ=255)

                                                      agccGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  2:1033052/91‑1 (MQ=255)

                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  2:2422353/93‑1 (MQ=255)

                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    <  2:3427697/93‑1 (MQ=255)

                                                      aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                                    >  1:4893894/1‑93 (MQ=255)

                                                       aCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgct                                   >  2:296592/1‑93 (MQ=255)

                                                          gTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGcc                                >  2:2113197/1‑93 (MQ=255)

                                                          gTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGcc                                <  1:5315407/93‑1 (MQ=255)

                                                           gggTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCg                               <  2:5270107/92‑1 (MQ=255)

                                                             gTAGCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGc                             >  2:5878869/1‑93 (MQ=255)

                                                                gCAATGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGAc                          >  2:2248667/1‑93 (MQ=255)

                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  2:4333130/93‑1 (MQ=255)

                                                                  aaTGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                        <  1:269463/93‑1 (MQ=255)

                                                                   aTGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAt                       >  2:1252620/1‑93 (MQ=255)

                                                                    tGTACACAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtg                      >  2:2089683/1‑93 (MQ=255)

                                                                       acacAACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTgg                   >  2:6034415/1‑93 (MQ=255)

                                                                        cacaACCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTggg                  <  2:4175366/93‑1 (MQ=255)

                                                                            aCCGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCaaa              >  1:5501464/1‑93 (MQ=255)

                                                                             ccGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAg             <  2:4445728/93‑1 (MQ=255)

                                                                              cGGGTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGt            <  2:3278883/93‑1 (MQ=255)

                                                                               gggTGCCTGCCAGTGGG‑GGTAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           >  1:6543987/1‑93 (MQ=255)

                                                                               gggTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           >  2:2223858/1‑93 (MQ=255)

                                                                               gggTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt           <  1:5918985/93‑1 (MQ=255)

                                                                                gctgcctgccAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGttt          <  1:2712381/91‑1 (MQ=255)

                                                                                 gTGCCTGCCAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTg         >  1:5280472/1‑93 (MQ=255)

                                                                                 gTGCCTGCAAGTGGG‑GGCAGCTGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTg         <  1:603229/93‑1 (MQ=255)

                                                                                      tgccAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc    >  1:1342371/1‑93 (MQ=255)

                                                                                       gccAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAcc   <  1:3084721/93‑1 (MQ=255)

                                                                                       gccAGTGGG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAcc   >  2:3759580/1‑93 (MQ=255)

GACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGT‑GG‑GGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCA  >  NC_000913/115872‑116051

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: `‑`