BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,078,245	C→G	G470G (GGG→GGC)	tktA ←	transketolase 1 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,078,245	0	C	G	99.6%	746.2 / NA	229	G470G (GGG→GGC)	tktA	transketolase 1 isozyme
Reads supporting (aligned to +/- strand): ref base C (0/0); major base G (105/123); minor base T (0/1); total (105/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

GCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078156‑3078334

gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGgcct                                                                                        >  1:6364722/1‑92 (MQ=255)

 cgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGt                                                                                       >  2:4729554/1‑93 (MQ=255)

 cgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGt                                                                                       <  1:4469159/93‑1 (MQ=255)

  ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                                                      >  1:4660117/1‑93 (MQ=255)

   gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCt                                                                                     <  1:4786703/93‑1 (MQ=255)

   gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCt                                                                                     <  1:4726877/93‑1 (MQ=255)

   gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCt                                                                                     >  2:5778796/1‑93 (MQ=255)

   gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCt                                                                                     <  2:5348219/93‑1 (MQ=255)

    aTTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCtt                                                                                    <  2:6258577/93‑1 (MQ=255)

        aaCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGcc                                                                                <  2:14484/93‑1 (MQ=255)

        aaCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGcc                                                                                >  1:2416156/1‑93 (MQ=255)

          ccTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGGCGGGCCGTCTTCGCCCa                                                                              >  2:4206375/1‑93 (MQ=255)

           cTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAg                                                                             <  2:516415/93‑1 (MQ=255)

                cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                                        <  1:2751867/93‑1 (MQ=255)

                cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                                        >  1:3558888/1‑93 (MQ=255)

                 acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                                                       <  2:2147144/93‑1 (MQ=255)

                 acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                                                       >  2:2057470/1‑93 (MQ=255)

                   acGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATgg                                                                     >  2:4852096/1‑93 (MQ=255)

                     ggACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAg                                                                   <  1:1458535/93‑1 (MQ=255)

                     ggACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAg                                                                   >  2:2241306/1‑93 (MQ=255)

                      gACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGt                                                                  >  2:3953005/1‑93 (MQ=255)

                      gACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGt                                                                  <  1:1300819/93‑1 (MQ=255)

                       aCGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGCGTc                                                                 >  1:3586482/1‑93 (MQ=255)

                         gCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                                                               <  2:313751/93‑1 (MQ=255)

                          ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                                              <  1:6092436/93‑1 (MQ=255)

                          ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                                              <  2:3380002/93‑1 (MQ=255)

                           cATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                                             <  1:6119687/93‑1 (MQ=255)

                           cATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                                             >  2:5566741/1‑93 (MQ=255)

                            aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCCGACCGATGGAGTCGTggg                                                            >  1:5731911/1‑93 (MQ=255)

                            aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            >  1:3704429/1‑93 (MQ=255)

                            aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            >  1:1597146/1‑93 (MQ=255)

                            aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            <  2:6306033/93‑1 (MQ=255)

                            aTGTAGACATGTTCGGGGCTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGCGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            >  1:5921684/1‑93 (MQ=255)

                              gTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtg                                                          <  1:2050162/93‑1 (MQ=255)

                               tAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtgt                                                         <  2:80961/93‑1 (MQ=255)

                               tAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtgt                                                         >  1:3069023/1‑93 (MQ=255)

                                 gACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTaa                                                       >  1:765900/1‑93 (MQ=255)

                                  aCATGTTCGGGGTGACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTaaa                                                      <  1:3870853/93‑1 (MQ=255)

                                   cATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                                                     >  1:5833538/1‑93 (MQ=255)

                                      gTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGCGTCGTGGGGGTGAACCAt                                                  >  2:3312262/1‑93 (MQ=255)

                                      gTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCAt                                                  <  2:6043714/93‑1 (MQ=255)

                                       ttCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATc                                                 <  2:2463833/93‑1 (MQ=255)

                                         cGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAc                                               <  2:1217836/93‑1 (MQ=255)

                                          ggggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAcc                                              <  2:2856178/93‑1 (MQ=255)

                                           gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             <  1:6019763/93‑1 (MQ=255)

                                           gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             <  2:4812262/93‑1 (MQ=255)

                                             gTTACGCTCAGAGAAGCGACCTGCTCAACCGGCGGGTGAGTCGGGCCGTCTTGGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGa                                           <  1:2560518/93‑1 (MQ=255)

                                             gTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGa                                           >  2:4212778/1‑93 (MQ=255)

                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  2:4177739/93‑1 (MQ=255)

                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  1:3146273/93‑1 (MQ=255)

                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  2:4898160/93‑1 (MQ=255)

                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          >  1:3951429/1‑93 (MQ=255)

                                               tACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACg                                         <  1:5145577/93‑1 (MQ=255)

                                                aCGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGc                                        <  1:3403775/93‑1 (MQ=255)

                                                 cgcgCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCt                                       <  1:639914/93‑1 (MQ=255)

                                                  gcgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAGACGATTAACTGACGCTg                                      >  2:4448694/1‑93 (MQ=255)

                                                    gcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGtt                                    <  2:1059107/93‑1 (MQ=255)

                                                    gcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGtt                                    >  1:51123/1‑93 (MQ=255)

                                                     cAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGttt                                   <  1:5527407/93‑1 (MQ=255)

                                                       gagaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca                                 <  1:4924189/93‑1 (MQ=255)

                                                         gaAGCGCCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc                               <  1:4396680/93‑1 (MQ=255)

                                                         gaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc                               >  2:651967/1‑93 (MQ=255)

                                                           aGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAg                             <  2:334942/93‑1 (MQ=255)

                                                            gCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgc                            <  2:1266993/93‑1 (MQ=255)

                                                            gCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgc                            >  1:5373640/1‑93 (MQ=255)

                                                              gACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCGTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcgc                          <  2:720660/93‑1 (MQ=255)

                                                               aCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCa                         >  1:2857781/1‑93 (MQ=255)

                                                               aCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCa                         >  2:4004044/1‑93 (MQ=255)

                                                                 cTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGc                       >  1:835486/1‑93 (MQ=255)

                                                                 cTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCCTGGGTTTAAACCATCACCTGACGCTGTTTCATAAGCGCAGc                       <  1:6271107/93‑1 (MQ=255)

                                                                  tGCTCAACCGGCTGGTGAGTCGGTCCGTCTTCGCCCAGACCGATGGAGTCGTGTGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGcc                      <  2:3694480/93‑1 (MQ=255)

                                                                  tGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGATGCTGTTTCATCAGCGCAGcc                      >  1:2695745/1‑93 (MQ=255)

                                                                  tGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGGGTAAACCATCACCTGACGCTGTTTTATCCGCGCAGcc                      >  1:720973/1‑93 (MQ=255)

                                                                   gCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCa                     >  1:4973015/1‑93 (MQ=255)

                                                                   gCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAACATCACCTGACGCTGTTTCATCAGCGCAGCCa                     >  1:424225/1‑93 (MQ=255)

                                                                    cTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGccct                    >  2:1996770/1‑91 (MQ=255)

                                                                    cTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAt                    >  2:3572106/1‑93 (MQ=255)

                                                                    cTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAt                    <  1:266781/93‑1 (MQ=255)

                                                                     tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   >  1:2833234/1‑93 (MQ=255)

                                                                     tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   >  2:1558319/1‑93 (MQ=255)

                                                                      cAACCGGCTGGTGAGTCGGGACGTCTTCGCCCAGACCGATGGAGTCGTGGGCGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtac                  >  2:2610570/1‑93 (MQ=255)

                                                                       aaCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacg                 <  1:4048376/93‑1 (MQ=255)

                                                                        aCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAACCCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgc                >  2:2951771/1‑92 (MQ=255)

                                                                        aCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgt                >  2:1308904/1‑93 (MQ=255)

                                                                         ccGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgta               <  1:1096694/93‑1 (MQ=255)

                                                                          cGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtac              <  1:2461376/93‑1 (MQ=255)

                                                                          cGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtac              <  1:1655845/93‑1 (MQ=255)

                                                                          cGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtac              <  1:4072054/93‑1 (MQ=255)

                                                                             cTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGc           >  1:5425109/1‑93 (MQ=255)

                                                                               ggTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGt         >  1:1798814/1‑93 (MQ=255)

                                                                                gTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGtt        >  2:4609624/1‑93 (MQ=255)

                                                                                gTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGtt        <  1:5650427/93‑1 (MQ=255)

                                                                                   aGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACg     >  2:1159499/1‑93 (MQ=255)

                                                                                    gTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGt    <  1:1918363/93‑1 (MQ=255)

                                                                                    gTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGt    <  2:5868866/93‑1 (MQ=255)

                                                                                     tCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTg   <  2:5610147/93‑1 (MQ=255)

                                                                                     tCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTg   <  1:1369589/93‑1 (MQ=255)

                                                                                      cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  2:3867131/93‑1 (MQ=255)

                                                                                      cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  2:6047565/93‑1 (MQ=255)

                                                                                      cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  2:2165992/93‑1 (MQ=255)

GCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078156‑3078334

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