BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,494,690	A→G	A17A (GCA→GCG)	nirC →	Nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,494,690	0	A	G	99.3%	866.0 / ‑3.7	279	A17A (GCA→GCG)	nirC	Nitrite transporter
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (239/38); total (241/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 290 total aligned reads displayed.

TTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAC  >  NC_000913/3494600‑3494779

ttGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGcgcg                                                                                         >  2:4113874/1‑93 (MQ=255)

 tGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTgcgcgc                                                                                        >  1:1026899/1‑93 (MQ=255)

   ttttATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCt                                                                                      >  1:6151289/1‑93 (MQ=255)

      tATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTAAAGACACTATTAATAAGTGTGCGCCTAACGCTGCGCGCCTTGCGCGCCTGtg                                                                                   >  1:4177319/1‑92 (MQ=255)

       aTTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCg                                                                                  >  2:5618757/1‑93 (MQ=255)

          ttttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCa                                                                               >  2:3607382/1‑93 (MQ=255)

           tttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACATACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCaa                                                                              >  2:3422604/1‑93 (MQ=255)

           tttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCaa                                                                              >  2:2741055/1‑93 (MQ=255)

           tttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCaa                                                                              >  1:6968422/1‑93 (MQ=255)

             tGTTTTTATTTTTTAGAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAAt                                                                            <  2:3921563/93‑1 (MQ=255)

             tGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAAt                                                                            >  1:2820809/1‑93 (MQ=255)

             tGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAAt                                                                            >  2:2619762/1‑93 (MQ=255)

              gTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATa                                                                           >  2:2251666/1‑93 (MQ=255)

              gTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATa                                                                           >  2:1176500/1‑93 (MQ=255)

               tttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATaa                                                                          >  1:5259772/1‑93 (MQ=255)

               tttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATaa                                                                          >  1:7383682/1‑93 (MQ=255)

                ttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAc                                                                         >  2:6906341/1‑93 (MQ=255)

                ttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAc                                                                         >  1:561915/1‑93 (MQ=255)

                ttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAc                                                                         >  2:6345256/1‑93 (MQ=255)

                 tttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAcc                                                                        <  1:2436498/93‑1 (MQ=255)

                 tttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGCGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAcc                                                                        <  1:5931379/93‑1 (MQ=255)

                     ttttttAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCt                                                                    >  2:731045/1‑93 (MQ=255)

                     ttttttAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCt                                                                    >  2:5557419/1‑93 (MQ=255)

                      tttttAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATCACCCGCTc                                                                   >  1:5814118/1‑93 (MQ=255)

                       ttttAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCg                                                                  >  2:444934/1‑93 (MQ=255)

                         ttAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGc                                                                >  2:6091250/1‑93 (MQ=255)

                          tAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCt                                                               >  2:1275616/1‑93 (MQ=255)

                          tAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCt                                                               >  1:1808255/1‑93 (MQ=255)

                          tAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCt                                                               >  1:6199868/1‑93 (MQ=255)

                           aaaGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCtt                                                              >  2:7928285/1‑93 (MQ=255)

                            aaGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCttt                                                             >  1:5307743/1‑93 (MQ=255)

                            aaGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCttt                                                             >  2:3494486/1‑93 (MQ=255)

                            aaGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCttt                                                             >  1:3706727/1‑93 (MQ=255)

                             aGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCtttt                                                            >  1:3984904/1‑93 (MQ=255)

                             aGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCtttt                                                            <  1:3288888/93‑1 (MQ=255)

                             aGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCtttt                                                            >  2:738857/1‑93 (MQ=255)

                              ggATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                           >  2:5174468/1‑93 (MQ=255)

                               gATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTgg                                                          >  2:6507013/1‑93 (MQ=255)

                                aTAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTggg                                                         >  1:7260352/1‑93 (MQ=255)

                                 tAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATCACCCGCTCGGCTTTTGGGt                                                        >  1:341487/1‑93 (MQ=255)

                                 tAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGt                                                        >  2:3057555/1‑93 (MQ=255)

                                 tAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGt                                                        >  2:1015285/1‑93 (MQ=255)

                                    tCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAg                                                     >  1:7908994/1‑93 (MQ=255)

                                    tCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAg                                                     >  1:1204851/1‑93 (MQ=255)

                                    tCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAg                                                     >  2:7604834/1‑93 (MQ=255)

                                      aaaTGTTTACAGACACTATTAATAAGTGGGCGGCTAATGTGGCGCTCACGTCGCGCCTGTCGGTAAATAACCCGTCCGGCTTTTGGGTTAgtt                                                   >  1:6038371/1‑91 (MQ=255)

                                       aaTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTc                                                  >  1:4402/1‑93 (MQ=255)

                                         tGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCg                                                >  1:1493696/1‑93 (MQ=255)

                                         tGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCg                                                >  2:2916092/1‑93 (MQ=255)

                                          gTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGc                                               >  1:6498098/1‑93 (MQ=255)

                                          gTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGc                                               >  1:5085627/1‑93 (MQ=255)

                                          gTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGc                                               >  1:3622168/1‑93 (MQ=255)

                                              aCAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATg                                           <  2:7757577/93‑1 (MQ=255)

                                               cAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATgg                                          >  1:1894668/1‑93 (MQ=255)

                                                 gACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcg                                        >  1:6896406/1‑93 (MQ=255)

                                                   cacTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcggg                                      >  2:5410065/1‑93 (MQ=255)

                                                    acTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgggg                                     >  2:1421536/1‑92 (MQ=255)

                                                     cTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgggcg                                    >  2:4891523/1‑93 (MQ=255)

                                                     cTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgggcg                                    >  1:4077100/1‑93 (MQ=255)

                                                     cTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgggcg                                    >  1:5534315/1‑93 (MQ=255)

                                                      tATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGgcgc                                   >  1:2287923/1‑93 (MQ=255)

                                                      tATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGgcgc                                   >  2:4252956/1‑93 (MQ=255)

                                                       aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGcgcg                                  >  2:5899323/1‑93 (MQ=255)

                                                       aTTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGcgcg                                  >  1:6330801/1‑93 (MQ=255)

                                                         taataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTa                                <  1:177011/93‑1 (MQ=255)

                                                          aataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAt                               <  1:763448/93‑1 (MQ=255)

                                                          aataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAt                               >  1:3044484/1‑93 (MQ=255)

                                                          aataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAt                               >  1:5375002/1‑93 (MQ=255)

                                                           ataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtg                              >  1:2667063/1‑93 (MQ=255)

                                                           ataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtg                              >  2:1608256/1‑93 (MQ=255)

                                                            taaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgg                             >  1:7090526/1‑92 (MQ=255)

                                                             aaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgtg                            >  1:4732575/1‑93 (MQ=255)

                                                              aGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTgg                           <  1:2145419/93‑1 (MQ=255)

                                                               gtgtGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGGATGTggg                          >  1:4857248/1‑93 (MQ=255)

                                                               gagtgCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTACGCggg                          >  1:391089/3‑93 (MQ=255)

                                                                 gtgCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTc                        <  1:7687381/93‑1 (MQ=255)

                                                                 gtgCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTc                        <  2:7067886/93‑1 (MQ=255)

                                                                 gtgCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTc                        >  2:641055/1‑93 (MQ=255)

                                                                    cGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTg                     >  1:6408446/1‑93 (MQ=255)

                                                                     ggCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTgg                    >  1:6575762/1‑93 (MQ=255)

                                                                     ggCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTgg                    >  2:3750679/1‑93 (MQ=255)

                                                                       cTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTgggg                  >  2:2000447/1‑92 (MQ=255)

                                                                       cTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTgggg                  >  2:1759886/1‑92 (MQ=255)

                                                                       cTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGa                  >  1:821555/1‑93 (MQ=255)

                                                                         aaCGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGatc                >  1:4666957/1‑93 (MQ=255)

                                                                         aaCGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGatc                >  2:6896707/1‑93 (MQ=255)

                                                                             cTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATcc            >  1:7512010/1‑93 (MQ=255)

                                                                              tGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCATATGTGGGTCTTGGGATCATCCt           >  2:5641342/1‑93 (MQ=255)

                                                                               gcgcgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTg          >  1:4247311/1‑93 (MQ=255)

                                                                               gcgcgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTg          >  1:7923469/1‑93 (MQ=255)

                                                                               gcgcgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGGATGTGGGTCTTGGGATCATCCTg          >  2:871462/1‑93 (MQ=255)

                                                                                  cgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGAtt       >  1:2034101/1‑93 (MQ=255)

                                                                                   gcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGAttt      <  1:3968988/93‑1 (MQ=255)

                                                                                   gcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGAttt      <  1:3888244/93‑1 (MQ=255)

                                                                                    cATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGAtttt     >  2:1892822/1‑93 (MQ=255)

                                                                                     aTTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGCCTTTGGATCATCCTGATTTTc    >  2:4700761/1‑93 (MQ=255)

                                                                                      ttGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGTTTATCCTGAttttgc   >  2:6648848/1‑91 (MQ=255)

                                                                                      ttGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCa   <  1:3547138/93‑1 (MQ=255)

                                                                                      ttGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCa   >  2:7162/1‑93 (MQ=255)

                                                                                       tGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAc  >  2:582239/1‑93 (MQ=255)

TTGTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAC  >  NC_000913/3494600‑3494779

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