BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,029,120	+C	coding (73/1191 nt)	b0967 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,029,116	1	.	C	98.0%	594.6 / 3.9	195	coding (77/1191 nt)	b0967	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (101/90); total (105/90)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029024‑1029203

gCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACgg                                                                                           >  2:1878415/1‑93 (MQ=255)

 cGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑c                                                                                         >  2:2310042/1‑93 (MQ=255)

  gAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc                                                                                        >  2:4219914/1‑93 (MQ=255)

  gAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc                                                                                        <  1:1736104/93‑1 (MQ=255)

  gAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc                                                                                        >  1:1981525/1‑93 (MQ=255)

   aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc‑c                                                                                      >  1:6439941/1‑93 (MQ=255)

   aaTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc‑c                                                                                      >  1:3146692/1‑93 (MQ=255)

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc‑cc                                                                                     >  1:6598367/1‑93 (MQ=255)

    aTAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc‑cc                                                                                     >  1:7650604/1‑93 (MQ=255)

     tAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc‑ccc                                                                                    >  1:2877774/1‑92 (MQ=255)

     tAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑cc‑ccc                                                                                    >  1:88039/1‑92 (MQ=255)

       aGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACCG‑cc‑cccgg                                                                                  >  2:2983829/1‑90 (MQ=255)

          gccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CC‑CCCGgaaa                                                                               >  2:2520115/1‑93 (MQ=255)

           ccgcGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGaaaa                                                                               <  2:7701530/93‑1 (MQ=255)

             gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGa                                                                             >  1:4706174/1‑93 (MQ=255)

             gcgTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGa                                                                             >  1:2352253/1‑93 (MQ=255)

                tGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGAccc                                                                          >  2:1585508/1‑93 (MQ=255)

                tGCTAACCATTTTCCCTGATGATCAACAATATCCATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGAcca                                                                          >  1:7084397/1‑92 (MQ=255)

                 gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCa                                                                         >  1:7456225/1‑93 (MQ=255)

                 gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCa                                                                         <  2:3451106/93‑1 (MQ=255)

                 gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCa                                                                         <  2:4060263/93‑1 (MQ=255)

                 gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCa                                                                         <  2:860606/93‑1 (MQ=255)

                  cTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCAc                                                                        <  2:1634519/93‑1 (MQ=255)

                   tAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACg                                                                       <  1:7312259/93‑1 (MQ=255)

                   tAACCATTTGCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACg                                                                       <  1:1849425/93‑1 (MQ=255)

                     aCCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGa                                                                     >  1:3868812/1‑93 (MQ=255)

                      ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGAt                                                                    >  1:4546275/1‑93 (MQ=255)

                      ccATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGAt                                                                    >  2:3564116/1‑93 (MQ=255)

                       taTTTTCACTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATg                                                                   <  1:792124/92‑1 (MQ=255)

                       cATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATg                                                                   <  1:1309006/93‑1 (MQ=255)

                       cATTTTCCCTAATGATCAACAATATCGATGGTTTCACCGAGGCTGTCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCTACGGATg                                                                   <  2:7761457/93‑1 (MQ=255)

                          tttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGcg                                                                <  1:1029576/93‑1 (MQ=255)

                           ttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGcga                                                               <  1:3514867/93‑1 (MQ=255)

                            tCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGcgac                                                              >  2:1159345/1‑93 (MQ=255)

                               cTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGaa                                                           >  1:5237342/1‑93 (MQ=255)

                                tgatgaTCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAg                                                          <  2:3686594/93‑1 (MQ=255)

                                tgatgaTCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CCACCGGAAAAGACCCACGGATGGCGACGAAg                                                          >  1:5391727/1‑93 (MQ=255)

                                  atgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTa                                                        >  1:705450/1‑93 (MQ=255)

                                  atgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTa                                                        <  2:3022459/93‑1 (MQ=255)

                                    gatCATCAATATCGATGGTTACACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCCCGGATGGCGACGAAGTAAt                                                      <  2:5371482/93‑1 (MQ=255)

                                    gatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAAt                                                      >  1:1893585/1‑93 (MQ=255)

                                    gatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAAt                                                      >  2:2369473/1‑93 (MQ=255)

                                     atCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                     <  1:3945259/93‑1 (MQ=255)

                                     atCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                     >  1:6863160/1‑93 (MQ=255)

                                      tCAACAATATCGATTGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                    <  2:6665831/93‑1 (MQ=255)

                                      tCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGa                                                    >  1:2632532/1‑93 (MQ=255)

                                        aacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGaat                                                  >  1:2981465/1‑91 (MQ=255)

                                        aacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGAtt                                                  <  2:443592/93‑1 (MQ=255)

                                            atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                              >  2:5981032/1‑93 (MQ=255)

                                            atatCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTc                                              >  1:2427887/1‑93 (MQ=255)

                                              atCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgc                                            >  1:4035258/1‑93 (MQ=255)

                                               tCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTcgcc                                           >  1:6267181/1‑92 (MQ=255)

                                                 gATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGcc                                         >  2:211292/1‑93 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        >  2:6305176/1‑93 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        <  1:6959763/93‑1 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        >  2:4907880/1‑93 (MQ=255)

                                                  aTGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGccc                                        >  1:3493147/1‑93 (MQ=255)

                                                     gTTTCACGGAGGCTGGTTTTCCCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     <  2:5575038/93‑1 (MQ=255)

                                                     gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     >  1:7151470/1‑93 (MQ=255)

                                                     gTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                     <  2:482306/93‑1 (MQ=255)

                                                      tttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTg                                    >  1:3312218/1‑93 (MQ=255)

                                                      gttCACCGAGGCTGGCTTGACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGGAGTAATGATTTTTCGCGCCCTTTg                                    <  1:4919450/92‑1 (MQ=255)

                                                       ttCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTgg                                   <  1:5732406/93‑1 (MQ=255)

                                                          aCCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTa                                >  1:5179555/1‑93 (MQ=255)

                                                           ccGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTaa                               <  2:6066424/93‑1 (MQ=255)

                                                            cGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAac                              >  1:435825/1‑93 (MQ=255)

                                                              aGGCTGGCTTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAacac                            >  1:1419476/1‑93 (MQ=255)

                                                               ggctggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                           <  2:3207308/93‑1 (MQ=255)

                                                               ggctggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                           <  2:6488256/93‑1 (MQ=255)

                                                               ggctggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACt                           <  1:5807010/93‑1 (MQ=255)

                                                                 ctggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaa                         <  1:6277150/93‑1 (MQ=255)

                                                                 ctggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaa                         <  2:2574354/93‑1 (MQ=255)

                                                                  tggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTaaa                        >  2:7399935/1‑93 (MQ=255)

                                                                   ggctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                       >  2:6275485/1‑93 (MQ=255)

                                                                    gctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACg                      >  2:7381731/1‑93 (MQ=255)

                                                                    gctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACg                      >  2:1782033/1‑93 (MQ=255)

                                                                    gctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACg                      >  2:4920554/1‑93 (MQ=255)

                                                                    gctTTACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACg                      <  1:6556433/93‑1 (MQ=255)

                                                                      tttACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTa                    >  1:6206868/1‑93 (MQ=255)

                                                                       ttACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTac                   >  1:1534152/1‑93 (MQ=255)

                                                                        tACCTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTaca                  <  2:7129813/93‑1 (MQ=255)

                                                                           cTTCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTc               >  2:7926015/1‑93 (MQ=255)

                                                                            ttCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAAACCTAAACGTacaccaa              >  2:3858919/1‑90 (MQ=255)

                                                                            ttCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCCCGCCCTTTGGCTAACACTAAACGTACACTCa              >  2:831533/1‑93 (MQ=255)

                                                                             tCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCAt             <  1:1224995/93‑1 (MQ=255)

                                                                             tCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCAt             <  2:1478713/93‑1 (MQ=255)

                                                                             tCCATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCAt             >  1:2069552/1‑93 (MQ=255)

                                                                              ccATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATa            <  1:4655130/93‑1 (MQ=255)

                                                                              ccATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATa            >  2:577986/1‑93 (MQ=255)

                                                                              ccATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATa            >  2:6971042/1‑93 (MQ=255)

                                                                               cATGCGGGCAAGGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           <  1:750670/93‑1 (MQ=255)

                                                                               cATGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATaa           <  2:388920/93‑1 (MQ=255)

                                                                                aTGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAt          >  2:629958/1‑93 (MQ=255)

                                                                                 tGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAtt         >  1:7857793/1‑93 (MQ=255)

                                                                                 tGCGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAAtt         >  2:6121993/1‑93 (MQ=255)

                                                                                   cGGGCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTa       <  1:6369651/93‑1 (MQ=255)

                                                                                     ggCAACGGCCC‑CCGTAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt     <  1:5537842/93‑1 (MQ=255)

                                                                                     ggCAACGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTAGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACt     >  1:1912338/1‑93 (MQ=255)

                                                                                        aaCGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  >  2:6626847/1‑93 (MQ=255)

                                                                                        aaCGGCCC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACtttt  >  1:2293034/1‑93 (MQ=255)

GCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGG‑CC‑CCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTT  >  NC_000913/1029024‑1029203

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