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.1%	337.1 / ‑4.5	107	A17A (GCA→GCG)	nirC	Nitrite transporter
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (82/24); total (82/25)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 108 total aligned reads displayed.

GTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAC  >  NC_000913/3494602‑3494779

gTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGcc                                                                                       <  1:594709/93‑1 (MQ=255)

    tATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTc                                                                                   >  1:2864204/1‑93 (MQ=255)

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

       tttttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGc                                                                                >  2:960734/1‑93 (MQ=255)

       tttttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGc                                                                                <  2:4141762/93‑1 (MQ=255)

       tttttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGc                                                                                >  1:3352602/1‑93 (MQ=255)

       tttttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGc                                                                                >  2:1624462/1‑93 (MQ=255)

       tttttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGc                                                                                >  1:4727395/1‑93 (MQ=255)

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

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

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

          ttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCaaa                                                                             >  1:340895/1‑93 (MQ=255)

          ttGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCaaa                                                                             >  2:4006651/1‑93 (MQ=255)

           tGTTTTTATTTTTTAAAGGTTAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAAt                                                                            >  1:2849130/1‑93 (MQ=255)

            gTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATa                                                                           <  1:4412419/93‑1 (MQ=255)

            gTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATa                                                                           <  1:4412418/93‑1 (MQ=255)

             tttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATaa                                                                          <  1:4644621/93‑1 (MQ=255)

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

             tttttATTTTTTAAAGGATAATCAAATGTTTACAGACACTAATAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATaa                                                                          >  2:1228375/1‑93 (MQ=255)

               tttATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAcc                                                                        >  1:3548294/1‑93 (MQ=255)

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

                 tATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAAcccc                                                                      >  1:1510510/1‑92 (MQ=255)

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

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

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

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

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

                            ggATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                           <  1:4099093/93‑1 (MQ=255)

                            ggATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                           <  2:3835642/93‑1 (MQ=255)

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

                            ggATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                           <  2:3558245/93‑1 (MQ=255)

                            ggATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCACATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTg                                                           >  2:3919617/1‑93 (MQ=255)

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

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

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

                             gATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCCCTCGGCTTTTgg                                                          >  1:2932472/1‑93 (MQ=255)

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

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

                                aaTCAAATGTTTAGAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTc                                                       >  1:3137873/1‑93 (MQ=255)

                                aaTCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTc                                                       >  1:4671273/1‑93 (MQ=255)

                                aaTCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTc                                                       >  2:1651133/1‑93 (MQ=255)

                                aaTCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTc                                                       <  1:3983305/93‑1 (MQ=255)

                                 aTCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCa                                                      >  2:2907303/1‑93 (MQ=255)

                                  tCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTccg                                                     >  2:1320329/1‑91 (MQ=255)

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

                                   cAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAgt                                                    >  1:2598300/1‑92 (MQ=255)

                                   cAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGc                                                    >  1:2598301/1‑93 (MQ=255)

                                    aaaTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCaccc                                                   >  1:1432853/1‑90 (MQ=255)

                                    aaaTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCt                                                   >  2:3671484/1‑93 (MQ=255)

                                    aaaTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCt                                                   >  2:523154/1‑93 (MQ=255)

                                    aaaTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCt                                                   >  2:3671485/1‑93 (MQ=255)

                                    aaaTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTAGGCAAATAACCCGCTCGGCTTTTGGGTCAGCt                                                   >  2:4038669/1‑93 (MQ=255)

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

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

                                      aTGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTcc                                                 >  1:2455539/1‑93 (MQ=255)

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

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

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

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

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

                                              aGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggc                                         >  2:4054399/1‑93 (MQ=255)

                                              aGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggc                                         >  2:2900526/1‑93 (MQ=255)

                                               gACACTATTAATAAGTGTGCGGCTAACGCTGTGCGCATTGTTCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcg                                        >  2:3844369/1‑93 (MQ=255)

                                               gACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcg                                        <  2:1910956/93‑1 (MQ=255)

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

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

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

                                                acacTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgg                                       <  2:2374793/93‑1 (MQ=255)

                                                acacTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATggcgg                                       >  2:250348/1‑93 (MQ=255)

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

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

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

                                                       taataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTa                                >  1:37497/1‑93 (MQ=255)

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

                                                        aataaGTGTGCGGCTAACGCTGTGCGCATTGCGCGCCTGTTGTTAAATAACCTGCTCGGCGTTTGTGTCAGCTCCGCCATGGCGGTcgcggtg                               >  2:4199134/1‑90 (MQ=255)

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

                                                         ataaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGCGTCAGCTCCGCCATGGCGGGCGCGTAtg                              >  2:117264/1‑93 (MQ=255)

                                                          taaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                             >  2:3123121/1‑93 (MQ=255)

                                                          taaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                             <  2:3205388/93‑1 (MQ=255)

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

                                                          aaaGTTTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTGTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                             <  1:3464245/92‑1 (MQ=255)

                                                          aaaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgt                             <  1:3464246/92‑1 (MQ=255)

                                                           aaGTGTGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTAtgtg                            <  1:3233892/93‑1 (MQ=255)

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

                                                             gtgtGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTggg                          >  2:2927255/1‑93 (MQ=255)

                                                             gtgtGCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTggg                          <  2:1546108/93‑1 (MQ=255)

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

                                                                tgCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCt                       >  1:3600861/1‑93 (MQ=255)

                                                                tgCGGCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCt                       <  1:3706789/93‑1 (MQ=255)

                                                                   ggCTAACGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGGCTTgg                    >  2:399149/1‑93 (MQ=255)

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

                                                                        aCGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGatca               >  2:4121187/1‑93 (MQ=255)

                                                                         cGCTGCGCGCATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGatcat              <  2:3880604/93‑1 (MQ=255)

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

                                                                              cgcgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGTCGCGTATGTGGGTCTTGGGATCATCCTgt         >  1:3893571/1‑92 (MQ=255)

                                                                               gcgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGAt        >  2:2675260/1‑93 (MQ=255)

                                                                               gcgcATTGCGCGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGAt        <  1:574465/93‑1 (MQ=255)

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

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

                                                                                     tGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAc  <  1:1554081/93‑1 (MQ=255)

GTTTTATTTTTTGTTTTTATTTTTTAAAGGATAATCAAATGTTTACAGACACTATTAATAAGTGTGCGGCTAACGCTGCGCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCAC  >  NC_000913/3494602‑3494779

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