BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,785,180	Δ1 bp	coding (413/495 nt)	b2655 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,785,180	0	T	.	99.5%	1019.2 / ‑2.8	217	coding (413/495 nt)	b2655	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (132/84); total (133/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

ACCTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCAA  >  NC_000913/2785092‑2785267

cccTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                                    <  2:862969/92‑1 (MQ=255)

aCCTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                                    >  1:2082002/1‑93 (MQ=255)

 ccTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                                   >  2:1377019/1‑93 (MQ=255)

 ccTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                                   >  1:3328107/1‑93 (MQ=255)

 ccTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                                   >  2:723962/1‑93 (MQ=255)

 ccTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                                   >  2:5861073/1‑93 (MQ=255)

  cTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACgg                                                                                  <  1:151674/93‑1 (MQ=255)

  cTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACgg                                                                                  <  1:6048286/93‑1 (MQ=255)

   tACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGa                                                                                 >  2:1844905/1‑93 (MQ=255)

   tACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGa                                                                                 >  1:3219061/1‑93 (MQ=255)

   tACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGa                                                                                 >  2:525455/1‑93 (MQ=255)

   tACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑ACGACGGa                                                                                 >  1:5821443/1‑93 (MQ=255)

      cagcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCa                                                                              <  2:1584816/93‑1 (MQ=255)

      cagcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCa                                                                              >  1:2819736/1‑93 (MQ=255)

       agcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                                             <  2:4137968/93‑1 (MQ=255)

        gcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGa                                                                            >  2:3669365/1‑93 (MQ=255)

          agATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGaaa                                                                          <  1:765480/93‑1 (MQ=255)

            aTCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACt                                                                        >  1:4362455/1‑93 (MQ=255)

            aTCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACt                                                                        >  2:2845548/1‑93 (MQ=255)

             acAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTc                                                                       <  1:1084727/92‑1 (MQ=255)

              cAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTcc                                                                      >  2:4648486/1‑93 (MQ=255)

                 agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGGATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAAATCCCga                                                                   >  2:3563183/1‑93 (MQ=255)

                  gaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGC‑TCCGACGGACCAGAAACTCCCgaa                                                                  >  1:3506374/1‑93 (MQ=255)

                    aaaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaaga                                                                >  2:5393800/1‑93 (MQ=255)

                     aaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                                               <  2:1684373/93‑1 (MQ=255)

                       tatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaaga                                                             >  1:3093497/1‑93 (MQ=255)

                         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaa                                                           <  1:4928199/93‑1 (MQ=255)

                         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaa                                                           >  1:4086086/1‑93 (MQ=255)

                          gTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaaa                                                          <  1:1193383/93‑1 (MQ=255)

                          gTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaaa                                                          >  2:4508178/1‑93 (MQ=255)

                           tCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAg                                                         >  2:2512906/1‑93 (MQ=255)

                            cTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                                        >  2:5729076/1‑93 (MQ=255)

                             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGATGGACCAGAAACTCCCGAAGAAGAAAAggg                                                       <  1:1626729/93‑1 (MQ=255)

                             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                                       >  1:668426/1‑93 (MQ=255)

                              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACGACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGt                                                      <  1:2314945/93‑1 (MQ=255)

                              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGt                                                      <  1:1134352/93‑1 (MQ=255)

                               aaTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGtt                                                     >  1:5018673/1‑93 (MQ=255)

                               aaTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAACAAGAAAAGGGtt                                                     >  1:221354/1‑93 (MQ=255)

                                aTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTa                                                    >  1:5657455/1‑93 (MQ=255)

                                aTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTa                                                    <  1:1938038/93‑1 (MQ=255)

                                 tttGCACGTTGCCATTATGGGGATATACAAATCCATCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGACAAGGGttct                                                   >  2:1731512/1‑91 (MQ=255)

                                    gCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTgg                                                >  2:5566131/1‑93 (MQ=255)

                                     cACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGt                                               >  2:711877/1‑93 (MQ=255)

                                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGtt                                              >  1:4519014/1‑93 (MQ=255)

                                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGtt                                              <  2:3868195/93‑1 (MQ=255)

                                       cGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGttt                                             >  1:1394276/1‑93 (MQ=255)

                                         ttGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCa                                           <  2:2073082/93‑1 (MQ=255)

                                          tGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCAt                                          <  1:1187203/93‑1 (MQ=255)

                                          tGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCAt                                          >  2:4375649/1‑93 (MQ=255)

                                          tGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCAt                                          <  1:3368939/93‑1 (MQ=255)

                                           gCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATg                                         >  1:4535225/1‑93 (MQ=255)

                                           gCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATg                                         >  1:3749502/1‑93 (MQ=255)

                                            ccATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGc                                        >  1:5336912/1‑93 (MQ=255)

                                               ttATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTa                                     >  2:5470652/1‑93 (MQ=255)

                                                tATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTAt                                    >  2:3205401/1‑93 (MQ=255)

                                                 aTGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc                                   >  2:2986137/1‑93 (MQ=255)

                                                 aTGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc                                   >  1:4228481/1‑93 (MQ=255)

                                                  tGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCa                                  >  2:6148617/1‑93 (MQ=255)

                                                  tGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCa                                  >  2:4825416/1‑93 (MQ=255)

                                                    ggTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAc                                >  2:4564859/1‑93 (MQ=255)

                                                     gTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAcc                               <  2:6020852/93‑1 (MQ=255)

                                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCAGAAGAAGAAAAGGGATACTGGTTTCATGCTTATCAAcca                              >  1:923863/1‑92 (MQ=255)

                                                        tataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa                            >  2:5551411/1‑93 (MQ=255)

                                                        tataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa                            >  2:962999/1‑93 (MQ=255)

                                                        tataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa                            >  1:6082446/1‑93 (MQ=255)

                                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGaa                           >  1:2203767/1‑93 (MQ=255)

                                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGaa                           <  1:3901903/93‑1 (MQ=255)

                                                          taCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAg                          >  2:1006256/1‑93 (MQ=255)

                                                             aaaTCCCTCCACTGACAACTTTCCGCG‑CCGAGGGACCAGAAACTCCCGAAGAAGAAAAGAGTTACTGGTTTCATGCTTATCAACCCGAAGATc                       <  1:1466813/93‑1 (MQ=255)

                                                              aaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCt                      >  1:3003431/1‑93 (MQ=255)

                                                               aTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCtt                     <  1:2411718/93‑1 (MQ=255)

                                                                 cccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTg                   <  2:213443/93‑1 (MQ=255)

                                                                 cccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTg                   >  1:4189758/1‑93 (MQ=255)

                                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGCCCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGt                  >  1:1791285/1‑93 (MQ=255)

                                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGt                  >  1:3401395/1‑93 (MQ=255)

                                                                    tCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAc                >  2:1127291/1‑93 (MQ=255)

                                                                     ccAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcc               >  1:1060137/1‑93 (MQ=255)

                                                                     ccAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcc               >  2:2851178/1‑93 (MQ=255)

                                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcca              <  2:1871502/93‑1 (MQ=255)

                                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcca              <  2:6059091/93‑1 (MQ=255)

                                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcca              <  1:2728719/93‑1 (MQ=255)

                                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccat             >  1:5094865/1‑93 (MQ=255)

                                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTCTAccat             >  2:5162955/1‑93 (MQ=255)

                                                                         tGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatcc           >  1:5903170/1‑93 (MQ=255)

                                                                          gACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatcca          <  1:2604145/93‑1 (MQ=255)

                                                                          gACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatcca          >  2:5896928/1‑93 (MQ=255)

                                                                          gACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatcca          >  1:3673783/1‑93 (MQ=255)

                                                                            cAACTTTCCGCG‑CCGACGTACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTTTACCATCCATa        <  1:2267116/93‑1 (MQ=255)

                                                                            cAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATa        <  1:1546587/93‑1 (MQ=255)

                                                                             aaCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATaa       <  2:3417529/93‑1 (MQ=255)

                                                                             aaCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATaa       >  1:4582540/1‑93 (MQ=255)

                                                                              aCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAAt      >  1:3297979/1‑93 (MQ=255)

                                                                              aCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAAt      >  1:2130907/1‑93 (MQ=255)

                                                                              aCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAAt      >  1:5112952/1‑93 (MQ=255)

                                                                               cTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATc     <  1:576208/93‑1 (MQ=255)

                                                                               cTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATc     <  1:3815566/93‑1 (MQ=255)

                                                                                tttCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATcc    >  1:1315012/1‑93 (MQ=255)

                                                                                tttCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGGTTCATGCTTATTAACCCGAAGATCTTTGTACCATCCATAATcc    >  2:4308575/1‑93 (MQ=255)

                                                                                 ttCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCa   >  1:870615/1‑93 (MQ=255)

                                                                                  tCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCaa  <  2:2406811/93‑1 (MQ=255)

ACCTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCAA  >  NC_000913/2785092‑2785267

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