BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	126,160	A→G	T156A (ACC→GCC)	aceF →	pyruvate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	126,160	0	A	G	99.7%	992.1 / ‑5.0	318	T156A (ACC→GCC)	aceF	pyruvate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (166/151); total (166/152)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 325 total aligned reads displayed.

AAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCG  >  NC_000913/126070‑126249

aaaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGcc                                                                                         >  1:7322081/1‑93 (MQ=255)

  aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGt                                                                                       <  2:4067085/93‑1 (MQ=255)

  aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGt                                                                                       <  2:4659403/93‑1 (MQ=255)

      ggCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGaaa                                                                                   <  2:182939/93‑1 (MQ=255)

         gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAagag                                                                                >  2:7831635/1‑93 (MQ=255)

         gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAagag                                                                                <  2:3564153/93‑1 (MQ=255)

         gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAagag                                                                                >  2:6715932/1‑93 (MQ=255)

          aTAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAgaga                                                                               <  2:5669330/93‑1 (MQ=255)

          aTAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAgaga                                                                               <  2:4785132/93‑1 (MQ=255)

           tAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAgagac                                                                              >  1:6165069/1‑92 (MQ=255)

             aaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCa                                                                            >  1:4361636/1‑93 (MQ=255)

              aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCaa                                                                           >  2:3652125/1‑93 (MQ=255)

               gTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCaaa                                                                          <  2:858009/93‑1 (MQ=255)

               gTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCaaa                                                                          <  1:4828354/93‑1 (MQ=255)

               gTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCaaa                                                                          <  2:5554507/93‑1 (MQ=255)

                ttGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAg                                                                         >  1:6068569/1‑93 (MQ=255)

                ttGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAg                                                                         >  1:5002101/1‑93 (MQ=255)

                ttGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAg                                                                         >  1:3092136/1‑93 (MQ=255)

                 tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGt                                                                        <  1:1488101/93‑1 (MQ=255)

                 tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGt                                                                        <  2:7172616/93‑1 (MQ=255)

                 tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGt                                                                        <  1:6321335/93‑1 (MQ=255)

                   aaGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGa                                                                      <  2:7644232/93‑1 (MQ=255)

                   aaGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGa                                                                      <  1:3519278/93‑1 (MQ=255)

                    aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGaa                                                                     <  1:3810458/93‑1 (MQ=255)

                    aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGaa                                                                     >  2:31980/1‑93 (MQ=255)

                    aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGaa                                                                     <  2:4635852/93‑1 (MQ=255)

                      cTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACg                                                                   <  1:503229/93‑1 (MQ=255)

                       tGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGt                                                                  <  2:2696079/93‑1 (MQ=255)

                       tGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGt                                                                  <  2:6442921/93‑1 (MQ=255)

                        gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTg                                                                 <  1:3991000/93‑1 (MQ=255)

                        gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTg                                                                 >  2:5728884/1‑93 (MQ=255)

                           cAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGt                                                              <  1:6433574/93‑1 (MQ=255)

                            aGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTg                                                             >  2:499969/1‑93 (MQ=255)

                            aGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTg                                                             >  2:2833702/1‑93 (MQ=255)

                            aGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATAAAAGTGAACGTGGGTg                                                             <  1:7619134/93‑1 (MQ=255)

                             gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGa                                                            >  2:1355792/1‑93 (MQ=255)

                             gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTAT‑GAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGAc                                                           <  1:3033894/93‑1 (MQ=255)

                              tCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGAc                                                           >  1:4055724/1‑93 (MQ=255)

                                gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACaa                                                         >  1:2469577/1‑93 (MQ=255)

                                 cTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACaaa                                                        <  2:3844280/93‑1 (MQ=255)

                                  tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                                                       >  2:439988/1‑93 (MQ=255)

                                   gATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgt                                                      >  1:2658816/1‑93 (MQ=255)

                                     tCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtgt                                                    <  2:257763/93‑1 (MQ=255)

                                      cacaGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTc                                                   <  2:610669/93‑1 (MQ=255)

                                      cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTc                                                   <  1:4136225/93‑1 (MQ=255)

                                       aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCt                                                  >  1:7748783/1‑93 (MQ=255)

                                       aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCt                                                  >  1:860750/1‑93 (MQ=255)

                                       aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCt                                                  >  2:5033141/1‑93 (MQ=255)

                                        ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa                                                 >  2:6491735/1‑93 (MQ=255)

                                        ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa                                                 <  1:5739918/93‑1 (MQ=255)

                                         cGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAc                                                >  1:2936959/1‑93 (MQ=255)

                                          gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAcc                                               <  1:5133584/93‑1 (MQ=255)

                                          gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAcc                                               >  2:3331156/1‑93 (MQ=255)

                                           tAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCg                                              >  1:3674752/1‑93 (MQ=255)

                                           tAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCg                                              >  2:2275659/1‑93 (MQ=255)

                                           tAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGACGTGAAATAGATCAAAGGAAACGCGGGGGACAAAGGTTCTACCg                                              >  1:7055921/1‑93 (MQ=18)

                                             gAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGc                                            >  2:899330/1‑93 (MQ=255)

                                             gAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGc                                            >  1:7827927/1‑93 (MQ=255)

                                              aaGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCt                                           <  1:1336062/93‑1 (MQ=255)

                                               aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTc                                          <  2:1641797/93‑1 (MQ=255)

                                                   gACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                      <  1:1035076/93‑1 (MQ=255)

                                                   gACAAGGCTTCTATGAAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                      <  1:2050642/93‑1 (MQ=255)

                                                    aCAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGa                                     <  1:2256609/93‑1 (MQ=255)

                                                    aCAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCCCTgg                                     >  1:4702829/1‑92 (MQ=255)

                                                     cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                                    >  2:6607665/1‑93 (MQ=255)

                                                       aGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAttt                                  >  2:2330764/1‑92 (MQ=255)

                                                       aGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTa                                  <  1:356931/93‑1 (MQ=255)

                                                        ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                                 <  2:2169759/93‑1 (MQ=255)

                                                        ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                                 <  2:5957585/93‑1 (MQ=255)

                                                          gttCTAGGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg                               <  1:6777460/92‑1 (MQ=255)

                                                          cTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg                               >  2:5361392/1‑93 (MQ=255)

                                                           ttCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGt                              >  1:3288583/1‑93 (MQ=255)

                                                           ttCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCCCCGATTATGGt                              >  2:5587162/1‑93 (MQ=37)

                                                            tCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTc                             <  2:2722437/93‑1 (MQ=255)

                                                             cTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATTGTCt                            >  2:3425611/1‑93 (MQ=255)

                                                             cTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCt                            <  1:581320/93‑1 (MQ=255)

                                                              tATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCtt                           <  2:6829204/93‑1 (MQ=255)

                                                              tATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCtt                           >  1:7536136/1‑93 (MQ=255)

                                                              tATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCtt                           <  2:1140599/93‑1 (MQ=255)

                                                                tGGAAGTTCCGGCTCCGTTTGCTTGCGCCGTGAAAGAGCTCAAAGTGAACGTTGGTGACAAAGTGTCTACCTGCTCCCTCAGTTTCGTCttgg                         >  2:2636329/1‑91 (MQ=255)

                                                                tGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCg                         <  1:3404357/93‑1 (MQ=255)

                                                                tGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGGCTTCg                         >  1:6986087/1‑93 (MQ=255)

                                                                    aGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGt                     >  1:5460142/1‑93 (MQ=255)

                                                                    aGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGt                     <  2:1913243/93‑1 (MQ=255)

                                                                        ccGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATTGTCTTCGGAGTcgcg                 >  2:2039829/1‑93 (MQ=38)

                                                                         cGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCgg                <  2:3061401/93‑1 (MQ=255)

                                                                          ggCTCCGTTTGCTGTCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg               <  2:7516210/93‑1 (MQ=255)

                                                                          ggCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg               <  2:6184371/93‑1 (MQ=255)

                                                                           gCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGt              >  1:3624738/1‑93 (MQ=255)

                                                                            cTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATTGTCTTCGGAGTCGCggggg             >  2:1502498/1‑91 (MQ=255)

                                                                               cGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAg          >  2:728012/1‑93 (MQ=255)

                                                                                gTTTGCTGGGGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGGAGTCGCGGGTGAAGc         >  2:5241987/1‑93 (MQ=255)

                                                                                gTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTTACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGc         <  1:5253360/93‑1 (MQ=255)

                                                                                 tttGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa        <  1:3381822/93‑1 (MQ=255)

                                                                                 tttGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa        <  1:7201068/93‑1 (MQ=255)

                                                                                  ttGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAg       <  1:7464598/93‑1 (MQ=255)

                                                                                  ttGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGGGAAGCAg       >  2:7403022/1‑93 (MQ=255)

                                                                                    gCTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgc     <  1:4081092/93‑1 (MQ=255)

                                                                                     cTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGGAGACGGGGGGGTAACAGgcg    >  1:6626764/1‑93 (MQ=37)

                                                                                     cTGGCGCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcg    >  1:1298843/1‑93 (MQ=255)

AAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCG  >  NC_000913/126070‑126249

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