BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	115,788	A→G	C438R (TGC→CGC)	hofB ←	HofB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	115,788	0	A	G	99.2%	792.3 / ‑3.5	259	C438R (TGC→CGC)	hofB	HofB
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (147/110); total (148/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 317 total aligned reads displayed.

CCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAA  >  NC_000913/115698‑115877

 caGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGcc                                                                                        <  1:7459067/93‑1 (MQ=255)

  cggcccAGAGTTGCTTACTCGCCATGCGGCACCCCCAGTACGCGGATTACCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCg                                                                                       <  2:6806132/89‑1 (MQ=255)

    cgcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGtt                                                                                     <  2:1880358/93‑1 (MQ=255)

     gcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGttt                                                                                    >  2:1473189/1‑93 (MQ=255)

      ccAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGtttt                                                                                   >  2:4573769/1‑93 (MQ=255)

       cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTc                                                                                  >  2:5980142/1‑93 (MQ=255)

            tttcTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaaa                                                                             <  1:2068182/90‑1 (MQ=255)

            ttGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaaa                                                                             <  2:6140853/93‑1 (MQ=255)

             tGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAg                                                                            <  2:3895841/93‑1 (MQ=255)

             tGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAg                                                                            <  1:4013140/93‑1 (MQ=255)

               cTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCg                                                                          <  2:6043608/93‑1 (MQ=255)

               cTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCg                                                                          >  1:6794624/1‑93 (MQ=255)

                ttACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGt                                                                         >  2:5674681/1‑93 (MQ=255)

                ttACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGt                                                                         >  1:7880006/1‑93 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  2:2783069/93‑1 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  1:6688691/93‑1 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  2:4220264/93‑1 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  1:2850717/93‑1 (MQ=255)

                 tACTCCCCATGCTGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCCGGCGGCTGTTTTCAAAAAGCGTa                                                                        <  2:121292/93‑1 (MQ=255)

                  aCTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTAc                                                                       >  2:2012687/1‑93 (MQ=255)

                       ccATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATa                                                                  >  1:2684122/1‑93 (MQ=255)

                          tGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTTCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAccc                                                               >  1:3404000/1‑93 (MQ=255)

                          tGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAccc                                                               >  1:3687436/1‑93 (MQ=255)

                            cGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGc                                                             <  2:5541236/93‑1 (MQ=255)

                            cGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGc                                                             >  2:3864204/1‑93 (MQ=255)

                              gCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCt                                                           <  1:2473859/93‑1 (MQ=255)

                              gCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCt                                                           <  1:7273319/93‑1 (MQ=255)

                                aTCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGt                                                         >  2:5420774/1‑93 (MQ=255)

                                aTCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGt                                                         >  2:2957361/1‑93 (MQ=255)

                                aTCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGt                                                         >  2:635435/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  1:3031551/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  1:3273825/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  2:6229610/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  1:7679798/1‑93 (MQ=255)

                                   cccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgt                                                      >  2:5458403/1‑93 (MQ=255)

                                    ccAGTACGCGGATTAACTCTTCACAGGTGGTTAAGCCTTGCTCTACGGTCAGGCGGCCGTTTTCAAAAAGCGTACGCATACACGCCTGTcgtg                                                     >  1:187433/1‑93 (MQ=255)

                                     cAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgc                                                    >  2:5180991/1‑93 (MQ=255)

                                       gtcggcgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgt                                                  <  2:6672440/89‑1 (MQ=255)

                                        tACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgtg                                                 >  1:4559010/1‑93 (MQ=255)

                                        tACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgtg                                                 >  2:7280276/1‑93 (MQ=255)

                                        tACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgtg                                                 >  1:3778832/1‑93 (MQ=255)

                                           gcgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTgcgtgcgt                                              >  2:2530885/1‑93 (MQ=255)

                                            cgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGtt                                             <  1:6127048/93‑1 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  2:7167255/1‑93 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  1:4165451/1‑93 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  1:617841/1‑93 (MQ=255)

                                              gATTACCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc                                           <  2:583251/93‑1 (MQ=255)

                                              gATTAACTCTTCAAAGGTGGTTAAGCTTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc                                           <  2:4626808/93‑1 (MQ=255)

                                              gATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc                                           >  1:5868329/1‑93 (MQ=255)

                                               cttAACTCTACAAGGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTcc                                          <  1:4716297/92‑1 (MQ=255)

                                               aTTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTcc                                          >  2:4777198/1‑93 (MQ=255)

                                                ttAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCa                                         >  2:3665295/1‑93 (MQ=255)

                                                 tAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAg                                        >  1:1032266/1‑93 (MQ=255)

                                                    ctctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGa                                     >  2:1319298/1‑93 (MQ=255)

                                                    ctctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGa                                     >  2:178473/1‑93 (MQ=255)

                                                      ctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAtt                                   >  2:14034/1‑93 (MQ=255)

                                                       ttCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                  <  2:3503983/93‑1 (MQ=255)

                                                        tCAAAGGTGGGTAAGCCTTGCTCTACGGCCAGGCGGCGGTTTTCAAAAACCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCa                                 <  2:4832512/93‑1 (MQ=255)

                                                         cAAAGGTGGTTAATCCTTGCTCTACGGCCAGGCGGCCTTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCaa                                <  1:4796549/93‑1 (MQ=255)

                                                         aaaaGGTGGTTAAGCCTTGCTCTACGGCCGGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCaa                                <  1:5605305/92‑1 (MQ=255)

                                                         aaaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGCGCGTGCGTTTCCAGCGATTCaa                                <  1:4941396/92‑1 (MQ=255)

                                                          aaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                               >  1:1388723/1‑93 (MQ=255)

                                                            aGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGt                             <  2:1153375/93‑1 (MQ=255)

                                                            aGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGt                             >  1:6064663/1‑93 (MQ=255)

                                                             ggtggtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTc                            >  1:356737/1‑93 (MQ=255)

                                                              gtggtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                           <  2:705919/93‑1 (MQ=255)

                                                              gtggtTAAGCCTGGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                           <  1:1411625/93‑1 (MQ=255)

                                                               tggtGAAGCGTTGCGTTCCGGCCAGGCGGCCTTTTTCAAAAAGCGTACGCATTCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCgg                          <  1:921046/93‑1 (MQ=255)

                                                               gggtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCgg                          <  1:5145595/92‑1 (MQ=255)

                                                                 gtTCAGCCTTGCTCACCGGCCAGGTGGGTTTTTTAAAAAAGCGTCCGCACACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  2:4451599/93‑1 (MQ=255)

                                                                 gtTAGGCCTTGCTCTACGGCCAGGGCGCTGGTTTCAAAAAGCGTTCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  1:6008858/93‑1 (MQ=255)

                                                                 gtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        >  2:5911302/1‑93 (MQ=255)

                                                                 gtTAAGCCTTCCTCTGCGGCGAGGCGCCCGTTTTCAAAATGCGTCCCCATACCCGCCTGTCGCGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  2:1281621/93‑1 (MQ=255)

                                                                  ttAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAt                       >  2:2144928/1‑93 (MQ=255)

                                                                   tAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAtt                      >  2:2283128/1‑93 (MQ=255)

                                                                   tAAGCCTTGCTCTACGGCAGGGCGGCTTTTTTCAAAAAGCTTACTCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAtt                      <  2:4052078/93‑1 (MQ=255)

                                                                    aaGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                     <  2:2739096/93‑1 (MQ=255)

                                                                     aGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAg                    >  2:6544268/1‑93 (MQ=255)

                                                                     aGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAg                    <  2:438044/93‑1 (MQ=255)

                                                                      gCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGc                   >  2:1562456/1‑93 (MQ=255)

                                                                      gCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGAATTAgg                   >  1:2332652/1‑92 (MQ=255)

                                                                         ttGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGa                >  1:6446942/1‑93 (MQ=255)

                                                                           gCTCTACGGCCAGGCGGCCTTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              <  1:6339160/93‑1 (MQ=255)

                                                                           gCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              <  2:3131609/93‑1 (MQ=255)

                                                                           gCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              >  1:2746499/1‑93 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            >  2:5876138/1‑93 (MQ=255)

                                                                             cctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            <  1:5432998/92‑1 (MQ=255)

                                                                              ctACGGCCGGGCGGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATaa           <  2:5754507/93‑1 (MQ=255)

                                                                               tccGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAg          <  2:7480938/91‑1 (MQ=255)

                                                                               tACGGCCGGGCGGCTGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAg          <  1:820510/93‑1 (MQ=255)

                                                                                 cGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        <  2:2894489/93‑1 (MQ=255)

                                                                                 cGGCCAGGCGGCCGTTTTAAAAATGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        <  2:4148143/93‑1 (MQ=255)

                                                                                   gCCAGGCGGCCTTTTTCAAAAAGCGTACCCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa      <  1:1842506/93‑1 (MQ=255)

                                                                                   gCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa      >  1:699942/1‑93 (MQ=255)

                                                                                    ccAGGCGGCCGTTTTCAAAAAGCTTACCCTTCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGAc     <  1:6554904/93‑1 (MQ=255)

                                                                                      aGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGa   >  1:596306/1‑93 (MQ=255)

                                                                                      aGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGa   >  2:7587899/1‑93 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGaa  <  1:5893521/93‑1 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGaa  >  2:6479767/1‑93 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGaa  >  1:1795885/1‑93 (MQ=255)

CCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAA  >  NC_000913/115698‑115877

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