BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,354,431	T→C	D235G (GAT→GGT)	sapA ←	SapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,354,431	0	T	C	98.9%	908.3 / ‑3.1	263	D235G (GAT→GGT)	sapA	SapA
Reads supporting (aligned to +/- strand): ref base T (0/3); new base C (132/128); total (132/131)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 280 total aligned reads displayed.

CCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTCC  >  NC_000913/1354341‑1354520

ccGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGATCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtcatca                                                                                         >  2:6325890/1‑93 (MQ=255)

   gTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGa                                                                                      <  2:7762222/93‑1 (MQ=255)

    tCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGAc                                                                                     >  1:3109334/1‑93 (MQ=255)

     cAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACg                                                                                    >  2:5457130/1‑93 (MQ=255)

      cggAGTTTCGAACGGCGCCCGGCGCCGCCGGAGCCAAAATCCGCCACTACCTGCGGCATTACCGGGTTACCGCGCCAGAAGTCACCATCACGt                                                                                   <  2:7411665/92‑1 (MQ=255)

      aGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGt                                                                                   <  1:8649053/93‑1 (MQ=255)

      aGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGt                                                                                   >  1:6557392/1‑93 (MQ=255)

       ggAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGtt                                                                                  <  1:8009994/93‑1 (MQ=255)

       ggAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGtt                                                                                  >  2:5133523/1‑93 (MQ=255)

        gAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTg                                                                                 <  1:6334642/93‑1 (MQ=255)

        gAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTg                                                                                 <  1:5621211/93‑1 (MQ=255)

        gAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTg                                                                                 <  1:482813/93‑1 (MQ=255)

         aGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGt                                                                                <  1:7304031/93‑1 (MQ=255)

          gatacGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTa                                                                               <  1:609718/89‑1 (MQ=255)

          gTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTa                                                                               <  2:5763414/93‑1 (MQ=255)

             tCGACAGGCGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGc                                                                            <  1:504472/93‑1 (MQ=255)

               gACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGa                                                                          >  1:3516429/1‑93 (MQ=255)

                      gTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGGCGTTGTAGGCGAATAAAtt                                                                   <  2:2468346/93‑1 (MQ=255)

                      gTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAAtt                                                                   <  2:1953329/93‑1 (MQ=255)

                      gTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAAtt                                                                   >  2:4632975/1‑93 (MQ=255)

                        ccGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGc                                                                 >  2:2875562/1‑93 (MQ=255)

                          ggTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGccc                                                               <  2:3595896/93‑1 (MQ=255)

                            tGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCgg                                                             <  2:1589996/93‑1 (MQ=255)

                            tGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCgg                                                             >  1:2067416/1‑93 (MQ=255)

                              tcgccgGAGCGTAAATCCTCCACTACCTGCGGCATTATCGGTTTACCGCGCAAAAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGgcg                                                           <  2:5500882/92‑1 (MQ=255)

                                gccgGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGcgcg                                                         <  2:5802146/93‑1 (MQ=255)

                                   ggAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTa                                                      >  2:8528299/1‑93 (MQ=255)

                                    gAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTAt                                                     <  2:3904069/93‑1 (MQ=255)

                                     aGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTAtt                                                    >  1:1153651/1‑93 (MQ=255)

                                      gCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTc                                                   <  1:5234750/93‑1 (MQ=255)

                                      gCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTc                                                   <  2:767271/93‑1 (MQ=255)

                                      gCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGCATTc                                                   <  2:5142594/93‑1 (MQ=255)

                                       ccTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTcc                                                  >  1:3933476/1‑93 (MQ=255)

                                        cTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCg                                                 <  1:5605374/93‑1 (MQ=255)

                                        cTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCg                                                 >  1:5355498/1‑93 (MQ=255)

                                        cTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCg                                                 <  2:6761826/93‑1 (MQ=255)

                                        cTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGCATTCCg                                                 >  2:4350445/1‑93 (MQ=255)

                                        cTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGCATTCCg                                                 >  2:40518/1‑93 (MQ=255)

                                         tAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGa                                                >  2:2121895/1‑93 (MQ=255)

                                         tAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGa                                                <  2:1062848/93‑1 (MQ=255)

                                         tAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGa                                                >  2:1819662/1‑93 (MQ=255)

                                           aaTCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACa                                              >  2:4571637/1‑93 (MQ=255)

                                            aTCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACaa                                             >  1:672448/1‑93 (MQ=255)

                                            aTCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACaa                                             >  2:2668676/1‑93 (MQ=255)

                                             tCCACCACTACCTGAGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAAc                                            <  2:8519499/93‑1 (MQ=255)

                                                accacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGa                                         <  1:1122105/93‑1 (MQ=255)

                                                accacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGa                                         <  1:3967248/93‑1 (MQ=255)

                                                accacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGa                                         >  2:8393191/1‑93 (MQ=255)

                                                 ccacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGAt                                        >  2:1696796/1‑93 (MQ=255)

                                                 ccacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGAt                                        <  1:2570504/93‑1 (MQ=255)

                                                  cacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATa                                       >  2:7201357/1‑93 (MQ=255)

                                                  cacTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATa                                       >  2:7554124/1‑93 (MQ=255)

                                                   cctcccTGCGGCATTAACGGTTTACCGCGCCAGAAGTCGCCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATAc                                      <  2:6415340/89‑1 (MQ=255)

                                                      aCCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGTCCTGCGCTGCATTCCTACAACTGATACggt                                   >  2:7367615/1‑92 (MQ=255)

                                                       ccTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGAc                                  <  1:1962862/93‑1 (MQ=255)

                                                        cTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGAcc                                 >  1:453916/1‑93 (MQ=255)

                                                        cTGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGAcc                                 <  2:4822936/93‑1 (MQ=255)

                                                         tGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCg                                >  2:4923569/1‑93 (MQ=255)

                                                         tGCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCg                                >  2:4883292/1‑93 (MQ=255)

                                                          gCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCgg                               <  1:8115692/93‑1 (MQ=255)

                                                          gCGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCgg                               <  2:4189762/93‑1 (MQ=255)

                                                           cGGCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGCATTCCGACAACTGATCCGGACCggg                              >  2:7969636/1‑92 (MQ=255)

                                                            ggCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTg                             <  1:3309733/93‑1 (MQ=255)

                                                             gCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGc                            <  1:2275729/93‑1 (MQ=255)

                                                             gCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGc                            >  2:7574802/1‑93 (MQ=255)

                                                             gCATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGc                            <  1:4009366/93‑1 (MQ=255)

                                                              cATTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATCCGGACCGGTGcc                           >  1:2800985/1‑93 (MQ=255)

                                                               aTTAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCAACAACTGATACGGACCGGTGCCg                          >  1:7739171/1‑93 (MQ=255)

                                                                ttAACGGTTTACCGCGTCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGCATTCAGACAACCGATCCGGCCCGGTGCCgc                         >  2:5931696/1‑92 (MQ=255)

                                                                ttAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCgc                         >  2:6522441/1‑92 (MQ=255)

                                                                ttAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGa                         <  2:6841599/93‑1 (MQ=255)

                                                                ttAACGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGCATTCCGACAACTGATACGGACCGATGCCgc                         >  1:7260397/1‑92 (MQ=255)

                                                                   aCGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCg                      >  1:6146016/1‑93 (MQ=255)

                                                                    cGGTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCgg                     >  2:831113/1‑93 (MQ=255)

                                                                     ggTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGt                    <  1:8584042/93‑1 (MQ=255)

                                                                     ggTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGt                    <  2:6232156/93‑1 (MQ=255)

                                                                      gTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGtt                   >  2:3679290/1‑93 (MQ=255)

                                                                      gTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGtt                   <  1:1305583/93‑1 (MQ=255)

                                                                      gTTTACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGtt                   >  1:1692104/1‑93 (MQ=255)

                                                                         tACCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGAc                <  2:1459564/93‑1 (MQ=255)

                                                                          aCCGCGCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACg               <  1:6247275/93‑1 (MQ=255)

                                                                           acgcgCCAGTAGTCATCATGTTCTTGTAGGCTAATGAAGTGCACGGCTTGGTATTCCGACAACTGACACTGACCGGTGCCGACCGGTTGACgg              <  1:6007654/92‑1 (MQ=255)

                                                                            cgcgCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGt             <  1:142743/93‑1 (MQ=255)

                                                                            cgcgCCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGt             <  1:8756117/93‑1 (MQ=255)

                                                                               gcCAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGa          >  2:6100450/1‑93 (MQ=255)

                                                                                ccAGTAGTCGCCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAg         <  2:123705/93‑1 (MQ=255)

                                                                                 cAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGt        <  1:7449498/93‑1 (MQ=255)

                                                                                 cAGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGt        <  2:1884409/93‑1 (MQ=255)

                                                                                 cAGAAGTCACCATGACGGTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGt        <  2:7123981/93‑1 (MQ=255)

                                                                                  aGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGtt       <  1:3247369/93‑1 (MQ=255)

                                                                                  aGAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGtt       <  2:843194/93‑1 (MQ=255)

                                                                                   gAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTg      <  1:8355466/93‑1 (MQ=255)

                                                                                   gAAGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTg      >  2:3300669/1‑93 (MQ=255)

                                                                                    aaGTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGc     <  1:6843489/93‑1 (MQ=255)

                                                                                      gTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTc   >  1:4257506/1‑93 (MQ=255)

                                                                                      gTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTc   <  1:6654404/93‑1 (MQ=255)

                                                                                      gTCACCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTc   >  2:4246373/1‑93 (MQ=255)

                                                                                       tcaccaTGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTcc  <  1:999488/93‑1 (MQ=255)

                                                                                       tcaccaTGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTcc  <  2:7262719/93‑1 (MQ=255)

                                                                                       tcaccaTGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTcc  <  2:5348107/93‑1 (MQ=255)

CCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGGCGCGGTATTCCGACAACTGATACGGACCGGTGCCGACCGGTTGACGGTCGAGTTGCTCC  >  NC_000913/1354341‑1354520

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