BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,431,925	A→G	I79V (ATC→GTC)	b4212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,431,925	0	A	G	98.6%	658.3 / ‑2.8	217	I79V (ATC→GTC)	b4212	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (120/93); total (122/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 224 total aligned reads displayed.

AGCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAC  >  NC_000913/4431835‑4432014

 gCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcg                                                                                        >  1:2720304/1‑93 (MQ=255)

 gCGACCTGCGGCGCAAAATCGGTGGGGGGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcg                                                                                        >  1:2400746/1‑93 (MQ=255)

   gACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGCcgcg                                                                                      >  2:346451/1‑93 (MQ=255)

       tGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATc                                                                                  <  1:1895572/93‑1 (MQ=255)

       tGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATc                                                                                  >  2:2141863/1‑93 (MQ=255)

        gcggcgCAAAATCGGTGGGGTGAGTCAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATcc                                                                                 >  2:3830739/1‑93 (MQ=255)

          ggcgCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCgg                                                                               <  2:3386601/93‑1 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggt                                                                              >  1:3132228/1‑93 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggt                                                                              <  1:1140475/93‑1 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggt                                                                              <  1:1775584/93‑1 (MQ=255)

            cgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtg                                                                             <  2:1390515/93‑1 (MQ=255)

             gcAAAATCGGTGTGGTGAGTAAAAAGATGCTTGCTCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtgg                                                                            <  2:1935463/93‑1 (MQ=255)

              cAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtggt                                                                           >  1:3393261/1‑93 (MQ=255)

                aaaTCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgc                                                                         <  2:230106/93‑1 (MQ=255)

                aaaTCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgc                                                                         >  2:1275348/1‑93 (MQ=255)

                  aTCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccg                                                                       >  2:389161/1‑93 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      <  1:3931963/93‑1 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      >  1:2057098/1‑93 (MQ=255)

                    cGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgcc                                                                     >  1:755292/1‑93 (MQ=255)

                       tGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGGGCCGCCGCa                                                                  >  1:4173673/1‑93 (MQ=255)

                        ggggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAt                                                                 <  1:3771830/93‑1 (MQ=255)

                        ggggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAt                                                                 <  2:3244293/93‑1 (MQ=255)

                         gggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtg                                                                >  2:2332847/1‑93 (MQ=255)

                         gggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtg                                                                <  1:1038294/93‑1 (MQ=255)

                         gggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtg                                                                >  2:4379606/1‑93 (MQ=255)

                         gggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtg                                                                >  1:295365/1‑93 (MQ=255)

                           gtgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtgtg                                                              <  2:1118303/93‑1 (MQ=255)

                           gtgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtgtg                                                              >  1:1169358/1‑93 (MQ=255)

                           gtgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtgtg                                                              <  1:302140/93‑1 (MQ=255)

                             gagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGa                                                            <  1:2530571/93‑1 (MQ=255)

                              agtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGaa                                                           <  1:861027/93‑1 (MQ=255)

                              agcgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGaa                                                           <  2:3737824/90‑1 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        >  2:4024796/1‑93 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        >  2:2699481/1‑93 (MQ=255)

                                  aaaaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTCTCGCGTATCCGGTGGTGCCGCCGCATGTGGAAtata                                                       <  2:3549215/93‑1 (MQ=255)

                                    aaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGc                                                     >  1:1306996/1‑93 (MQ=255)

                                     aGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAgac                                                    >  1:327871/1‑91 (MQ=255)

                                     aGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGcc                                                    >  2:1785944/1‑93 (MQ=255)

                                     aGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGcc                                                    <  1:3344060/93‑1 (MQ=255)

                                      gATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGTAATATAGCCt                                                   >  1:4358530/1‑93 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  <  1:4607513/93‑1 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  >  1:21569/1‑93 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  >  1:4285474/1‑93 (MQ=255)

                                        tGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCa                                                 <  2:709111/93‑1 (MQ=255)

                                           ttGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCGGCCGCATGTGGAATATAGCCTCAcgc                                              >  2:2161440/1‑93 (MQ=255)

                                           ttGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgc                                              >  1:1209024/1‑93 (MQ=255)

                                               gcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCt                                          >  1:162385/1‑93 (MQ=255)

                                                 aGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTgg                                        >  1:2374024/1‑93 (MQ=255)

                                                 aGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTgg                                        <  2:458602/93‑1 (MQ=255)

                                                   tCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGc                                      >  2:2947140/1‑93 (MQ=255)

                                                     gTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGa                                    >  1:2104896/1‑93 (MQ=255)

                                                     gTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGa                                    >  1:597257/1‑93 (MQ=255)

                                                     gTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGa                                    >  1:3468852/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:520982/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:2021158/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:1952770/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  1:3513423/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:1640345/1‑93 (MQ=255)

                                                       tACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCCGGGCGAGc                                  <  2:4572811/93‑1 (MQ=255)

                                                       tACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGc                                  >  1:1326796/1‑93 (MQ=255)

                                                        aCAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCa                                 <  1:4441024/93‑1 (MQ=255)

                                                           gtcGTAAGATCAGGATGGTTTTCTTAATCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGt                              <  2:4673988/91‑1 (MQ=255)

                                                           ggCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAggg                              >  1:562510/1‑92 (MQ=255)

                                                             cGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGa                            >  2:1850259/1‑93 (MQ=255)

                                                              gTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGACTCACGCCGCTGGGCGAGCAGGTGAg                           >  1:3544010/1‑93 (MQ=255)

                                                              gTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGGGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAg                           >  2:1232460/1‑93 (MQ=255)

                                                                  gAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaa                       >  1:2822200/1‑93 (MQ=255)

                                                                  gAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaa                       >  2:1103054/1‑93 (MQ=255)

                                                                   aaCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaa                      >  1:2771004/1‑93 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaaa                     >  2:1908166/1‑93 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaaa                     <  2:569234/93‑1 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaaa                     >  2:2690520/1‑93 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGCCGAGCAGGGGAGCGaaaa                     >  1:1263704/1‑93 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGCAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaaa                     >  1:95854/1‑93 (MQ=255)

                                                                     cAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAg                    <  2:3094004/93‑1 (MQ=255)

                                                                      aGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAgg                   <  2:1306589/93‑1 (MQ=255)

                                                                      aGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAgg                   <  2:3625942/93‑1 (MQ=255)

                                                                      aGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAgg                   <  1:3980371/93‑1 (MQ=255)

                                                                      aGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGCGAGCGAAAAgg                   >  2:144846/1‑93 (MQ=255)

                                                                        gATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 <  1:4627030/93‑1 (MQ=255)

                                                                        gATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 <  1:3840584/93‑1 (MQ=255)

                                                                         aTGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtg                >  2:4113507/1‑93 (MQ=255)

                                                                         aTGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtg                <  2:733107/93‑1 (MQ=255)

                                                                         aTGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtg                >  2:2418177/1‑93 (MQ=255)

                                                                         aTGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtg                <  1:324758/93‑1 (MQ=255)

                                                                         aTGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtg                >  1:2496315/1‑93 (MQ=255)

                                                                          tGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgc               <  2:2846233/93‑1 (MQ=255)

                                                                          tGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgc               <  2:59311/93‑1 (MQ=255)

                                                                          tGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgc               >  1:298839/1‑93 (MQ=255)

                                                                          tGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgc               <  2:4418428/93‑1 (MQ=255)

                                                                          tGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgc               >  1:2917091/1‑93 (MQ=255)

                                                                             ttttCTTAACGGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgctgc            <  2:4072367/93‑1 (MQ=255)

                                                                              tttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCa           <  1:1068986/93‑1 (MQ=255)

                                                                               ttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCAc          <  1:30962/93‑1 (MQ=255)

                                                                               ttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCAc          >  2:283261/1‑93 (MQ=255)

                                                                                 cTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTg        >  2:2588605/1‑93 (MQ=255)

                                                                                  ttAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTgg       >  1:2026958/1‑93 (MQ=255)

                                                                                     aCCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCgg    >  1:1515531/1‑93 (MQ=255)

                                                                                      ccGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCggc   >  1:2662063/1‑92 (MQ=255)

                                                                                       cGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAc  <  1:21991/93‑1 (MQ=255)

AGCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAC  >  NC_000913/4431835‑4432014

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