BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,772,080	G→A	A90T (GCC→ACC)	ydiB →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,772,080	0	G	A	99.5%	719.8 / ‑5.3	199	A90T (GCC→ACC)	ydiB	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (92/106); total (93/106)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 212 total aligned reads displayed.

AAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGT  >  NC_000913/1771990‑1772169

 aaaTGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCa                                                                                        <  2:2664463/93‑1 (MQ=255)

 aaaTGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCa                                                                                        >  2:2954892/1‑93 (MQ=255)

 aaaTGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCa                                                                                        >  2:1508304/1‑93 (MQ=255)

  aaTGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCAt                                                                                       <  2:3775884/93‑1 (MQ=255)

   aTGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCAta                                                                                      >  2:3520038/1‑92 (MQ=255)

       gcgGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAaca                                                                                  <  2:1831152/93‑1 (MQ=255)

       gcgGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAaca                                                                                  <  1:1384713/93‑1 (MQ=255)

       gcgGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAaca                                                                                  <  1:81286/93‑1 (MQ=255)

        cgGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAacac                                                                                 <  2:5541671/93‑1 (MQ=255)

          gAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCa                                                                               >  1:581242/1‑93 (MQ=255)

          gAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCa                                                                               >  2:5847965/1‑93 (MQ=255)

           aaCTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCAt                                                                              <  2:3386502/93‑1 (MQ=255)

             cTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCg                                                                            <  1:6695691/93‑1 (MQ=255)

              tGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGt                                                                           >  2:5136968/1‑93 (MQ=255)

               ggTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGtt                                                                          <  1:4442134/93‑1 (MQ=255)

                  gtATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAat                                                                       <  2:2895475/93‑1 (MQ=255)

                    aTCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatga                                                                     <  1:2381723/93‑1 (MQ=255)

                     tCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgat                                                                    >  1:5968003/1‑93 (MQ=255)

                     tCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgat                                                                    >  2:4645058/1‑93 (MQ=255)

                      cGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACAGCAGCGGCCAAACTGGTGGGGAACATCAACACCATCGTTAatgatg                                                                   >  1:1736919/1‑93 (MQ=255)

                       gATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatga                                                                  <  2:1644595/93‑1 (MQ=255)

                       gATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatga                                                                  <  1:1557762/93‑1 (MQ=255)

                       gATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatga                                                                  <  1:6654282/93‑1 (MQ=255)

                       gATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatga                                                                  <  2:874060/93‑1 (MQ=255)

                        aTGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatgat                                                                 <  2:2620028/93‑1 (MQ=255)

                        aTGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatgat                                                                 >  1:2738772/1‑93 (MQ=255)

                        aTGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatgat                                                                 <  1:2398381/93‑1 (MQ=255)

                        aTGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatgat                                                                 <  1:517389/93‑1 (MQ=255)

                        aTGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAatgatgat                                                                 >  2:3618141/1‑93 (MQ=255)

                          gCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATgg                                                               <  2:5659578/93‑1 (MQ=255)

                          gCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATgg                                                               >  2:835522/1‑93 (MQ=255)

                           ccGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGc                                                              <  2:495202/93‑1 (MQ=255)

                             gAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTa                                                            >  1:367041/1‑93 (MQ=255)

                               acaaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATc                                                          <  2:2266528/93‑1 (MQ=255)

                               acaaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATc                                                          <  2:6264250/93‑1 (MQ=255)

                               acaaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATc                                                          >  2:1554158/1‑93 (MQ=255)

                                 aaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTg                                                        >  2:468800/1‑93 (MQ=255)

                                 aaacaaCTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTg                                                        >  1:4162470/1‑93 (MQ=255)

                                  aacaacTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATGGTTAATGATGATGGCTATTTGc                                                       >  1:4060097/1‑93 (MQ=255)

                                  aacaacTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGc                                                       >  2:1922204/1‑93 (MQ=255)

                                   acaacTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCg                                                      >  1:3225579/1‑93 (MQ=255)

                                    caacTGGCGTGTTAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGt                                                     <  1:4243183/93‑1 (MQ=255)

                                    caacTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGt                                                     <  2:3306370/93‑1 (MQ=255)

                                    caacTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGt                                                     <  1:1050093/93‑1 (MQ=255)

                                     aacTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTg                                                    <  1:5590301/93‑1 (MQ=255)

                                      acTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTgg                                                   >  2:5963193/1‑93 (MQ=255)

                                      acTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTgg                                                   >  1:3071477/1‑93 (MQ=255)

                                      acTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTgg                                                   >  1:5456942/1‑93 (MQ=255)

                                       cTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGc                                                  >  1:351160/1‑93 (MQ=255)

                                            gtgtGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATaa                                             >  1:6170645/1‑93 (MQ=255)

                                            gtgtGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATaa                                             >  2:3421526/1‑93 (MQ=255)

                                            gtgtGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATaa                                             <  1:2236729/93‑1 (MQ=255)

                                             tgtgAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAac                                            >  1:1268945/1‑93 (MQ=255)

                                               tgAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAacac                                          <  1:5740356/93‑1 (MQ=255)

                                               tgAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAacac                                          <  1:4353998/93‑1 (MQ=255)

                                                gAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACAcc                                         <  2:6487818/93‑1 (MQ=255)

                                                 aaTATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCg                                        >  1:2921002/1‑93 (MQ=255)

                                                  atatGTTGATGAATTAGCACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGa                                       >  1:3377866/1‑93 (MQ=255)

                                                  atatGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGa                                       <  2:4277377/93‑1 (MQ=255)

                                                     tGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACgg                                    <  2:4350018/93‑1 (MQ=255)

                                                       ttGATGAATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCa                                  >  2:3454647/1‑93 (MQ=255)

                                                        tgatgaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCAc                                 <  2:4530056/93‑1 (MQ=255)

                                                         gatgaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACg                                >  1:2615100/1‑93 (MQ=255)

                                                          atgaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACgg                               <  2:2488396/93‑1 (MQ=255)

                                                          atgaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACgg                               <  1:4654770/93‑1 (MQ=255)

                                                          atgaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACgg                               <  1:3167464/93‑1 (MQ=255)

                                                           tgaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACggg                              <  1:6053158/93‑1 (MQ=255)

                                                            gaATTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCAACGGCACGGGc                             <  2:682155/93‑1 (MQ=255)

                                                             aaTTAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGcc                            <  1:918376/93‑1 (MQ=255)

                                                               ttAACACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCat                          >  2:663656/1‑93 (MQ=255)

                                                                 aaCACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCatat                        >  2:3358186/1‑93 (MQ=255)

                                                                 aaCACCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCatat                        >  2:4044184/1‑93 (MQ=255)

                                                                  acacCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATAtt                       <  1:5423973/93‑1 (MQ=255)

                                                                   cacCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTc                      >  2:2098475/1‑93 (MQ=255)

                                                                   cacCAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTc                      >  1:2357093/1‑93 (MQ=255)

                                                                      cAGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTcgcg                   >  2:700047/1‑93 (MQ=255)

                                                                       aGCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCgcgc                  >  2:5280761/1‑93 (MQ=255)

                                                                        gCTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGcc                 >  1:4413609/1‑93 (MQ=255)

                                                                        gCTGCAAAACTGGTGGGGACCATCAACACCATCGTTAATGATGAGGGCTATCTGCGTGTCTATAACACCGACGGCACGGGCCATATTCGCGcc                 <  1:3459005/93‑1 (MQ=255)

                                                                        acTGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGcc                 <  2:1415854/92‑1 (MQ=255)

                                                                          tGCCAAGCTGGTGGGGACCATCAACCCCATCGTTAATGATGATGACTATCTGCGTGGCAATAACACCGACGGCCCGGGCCATATTCGCGCCAt               <  1:1575175/93‑1 (MQ=255)

                                                                          tGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCAt               <  1:1094911/93‑1 (MQ=255)

                                                                          tGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCAt               >  2:2514270/1‑93 (MQ=255)

                                                                          tGCCAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCAt               >  2:1808598/1‑93 (MQ=255)

                                                                           cccAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCAtt              <  2:237830/92‑1 (MQ=255)

                                                                             cAAACTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTaa            >  2:5418430/1‑93 (MQ=255)

                                                                              aaaCTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTaaa           <  2:5342089/93‑1 (MQ=255)

                                                                              aaaCTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTaaa           <  2:1089654/93‑1 (MQ=255)

                                                                               aaCTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAag          >  2:6152086/1‑93 (MQ=255)

                                                                                aCTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAaga         <  1:1139657/93‑1 (MQ=255)

                                                                                 cTGGTGGGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAagag        <  1:4906128/93‑1 (MQ=255)

                                                                                  tggtggGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACCCCGACGGCACGGGCCATATTCGCGCCATTAAagaga       <  1:2676731/93‑1 (MQ=255)

                                                                                  tggtggGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAagaga       <  1:5632843/93‑1 (MQ=255)

                                                                                  tggtggGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAagaga       <  2:6073214/93‑1 (MQ=255)

                                                                                  tggtggGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAagaga       <  1:4775532/93‑1 (MQ=255)

                                                                                   ggtggGGACCATCAACACCATCGTTAATGATGATGGCTAACTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAagagag      <  1:6425181/93‑1 (MQ=255)

                                                                                    gtggGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGc     <  1:4567025/93‑1 (MQ=255)

                                                                                    gtggGGACCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGc     <  1:3336400/93‑1 (MQ=255)

                                                                                      ggggaCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCgg   >  1:1865408/1‑93 (MQ=255)

                                                                                       gggaCCATCAACGCCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGt  >  2:4080149/1‑93 (MQ=255)

AAAATGCGCGGAACTGGTGTATCGATGCCGAACAAACAACTGGCGTGTGAATATGTTGATGAATTAACACCAGCTGCCAAACTGGTGGGGGCCATCAACACCATCGTTAATGATGATGGCTATCTGCGTGGCTATAACACCGACGGCACGGGCCATATTCGCGCCATTAAAGAGAGCGGT  >  NC_000913/1771990‑1772169

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