BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,263,021	G→A	D29N (GAT→AAT)	hemA →	HemA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,263,021	0	G	A	99.5%	678.7 / ‑6.6	192	D29N (GAT→AAT)	hemA	HemA
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (114/77); total (115/77)
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.99e-01

Only 100 of 198 total aligned reads displayed.

GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110

  aGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  1:2729302/93‑1 (MQ=255)

  aGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  2:1453811/93‑1 (MQ=255)

   gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      <  2:7603202/93‑1 (MQ=255)

    aCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAgg                                                                                     >  1:1751239/1‑93 (MQ=255)

     cATGACCCTTTTAGCACTCGGTATCAACCATAAAGCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    <  1:2964068/93‑1 (MQ=255)

     cATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    >  2:6468133/1‑93 (MQ=255)

     cATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    >  2:3214884/1‑93 (MQ=255)

       tGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgcgc                                                                                  >  1:3764075/1‑93 (MQ=255)

       tGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgcgc                                                                                  >  1:5469646/1‑93 (MQ=255)

         aCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtt                                                                                >  1:3482824/1‑93 (MQ=255)

         aCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtt                                                                                <  2:1649865/93‑1 (MQ=255)

         aCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtt                                                                                <  1:5410869/93‑1 (MQ=255)

          cccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTg                                                                               <  1:3376288/93‑1 (MQ=255)

           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGa                                                                              >  1:1795656/1‑93 (MQ=255)

           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGa                                                                              >  2:4511521/1‑93 (MQ=255)

            cTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGAc                                                                             >  2:7483425/1‑93 (MQ=255)

             ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            >  1:3109544/1‑93 (MQ=255)

             ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            >  2:2799704/1‑93 (MQ=255)

             ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            <  2:3335001/93‑1 (MQ=255)

               ttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGc                                                                          >  1:757475/1‑93 (MQ=255)

               ttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGc                                                                          <  2:4403083/93‑1 (MQ=255)

                tAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGcc                                                                         <  1:1058414/93‑1 (MQ=255)

                tAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGcc                                                                         >  2:6178361/1‑93 (MQ=255)

                  gCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTg                                                                       >  2:1030842/1‑93 (MQ=255)

                   cACTCGGTATTAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCGCAATCAGGCGCTTGACAGCCTGc                                                                      >  1:7013852/1‑93 (MQ=255)

                   cACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGc                                                                      >  2:523836/1‑93 (MQ=255)

                    aCTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     <  2:1990072/93‑1 (MQ=255)

                    aCTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     >  1:7170926/1‑93 (MQ=255)

                    aCTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     <  1:608626/93‑1 (MQ=255)

                     cTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCtt                                                                    >  2:1803640/1‑93 (MQ=255)

                     cTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCtt                                                                    >  1:6733986/1‑93 (MQ=255)

                      tCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTg                                                                   >  1:6472810/1‑93 (MQ=255)

                      tCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTg                                                                   >  1:3268491/1‑93 (MQ=255)

                       cGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgc                                                                  <  1:7369323/93‑1 (MQ=255)

                         gTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgcgc                                                                >  2:7332327/1‑93 (MQ=255)

                          tATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCa                                                               <  1:7319731/93‑1 (MQ=255)

                            tccaCCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGc                                                             <  1:7038436/90‑1 (MQ=255)

                            tCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTTGCCGGATAAGCTCAATCACGCGCTTGACAGCCTGCTTGCGCAGc                                                             >  2:4319607/1‑93 (MQ=255)

                            tCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGc                                                             >  1:4009305/1‑93 (MQ=255)

                            tCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGc                                                             >  1:7905037/1‑93 (MQ=255)

                             cAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGcc                                                            >  2:6283681/1‑93 (MQ=255)

                             cAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGcc                                                            >  2:5973117/1‑93 (MQ=255)

                             cAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGcc                                                            <  2:4164905/93‑1 (MQ=255)

                              aaCCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCg                                                           >  2:4971977/1‑93 (MQ=255)

                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGAt                                                         <  2:3175301/93‑1 (MQ=255)

                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGAt                                                         >  1:1514534/1‑93 (MQ=255)

                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGAt                                                         <  1:5316740/93‑1 (MQ=255)

                                   tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      >  2:4605131/1‑93 (MQ=255)

                                   tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      >  2:1045774/1‑93 (MQ=255)

                                    aaaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTg                                                     >  1:1902889/1‑93 (MQ=255)

                                      aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   <  2:1236260/93‑1 (MQ=255)

                                      aaCGGCACCTGTAGCGCTGAGAGATCGTGTTTCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   <  2:2226640/93‑1 (MQ=255)

                                       aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  >  2:71267/1‑93 (MQ=255)

                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 <  1:3891030/93‑1 (MQ=255)

                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 >  1:14005/1‑93 (MQ=255)

                                          gCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggc                                               >  1:2093733/1‑93 (MQ=255)

                                           cACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcg                                              >  2:7878301/1‑93 (MQ=255)

                                            aCCTGTATCGCTGCGAGAGCGTCTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcgg                                             <  2:4257799/93‑1 (MQ=255)

                                             ccTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcggc                                            <  1:5874804/93‑1 (MQ=255)

                                             ccTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcggc                                            <  1:4472378/93‑1 (MQ=255)

                                              cTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGgcggcg                                           <  2:3080515/93‑1 (MQ=255)

                                              cTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGacggcg                                           <  2:4708829/93‑1 (MQ=255)

                                               tGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgt                                          >  2:2927474/1‑93 (MQ=255)

                                                gTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtg                                         >  1:381988/1‑93 (MQ=255)

                                                gTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtg                                         <  2:798205/93‑1 (MQ=255)

                                                   tCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCACGGCGGCgggggg                                      >  2:5498799/1‑91 (MQ=255)

                                                     gCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCt                                    >  1:3624614/1‑93 (MQ=255)

                                                      cTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTg                                   >  2:5637300/1‑93 (MQ=255)

                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  >  1:2265030/1‑93 (MQ=255)

                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  >  1:2512499/1‑93 (MQ=255)

                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  >  2:5847809/1‑93 (MQ=255)

                                                        gCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTc                                 >  1:2055495/1‑93 (MQ=255)

                                                        gCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTc                                 >  2:6070179/1‑93 (MQ=255)

                                                          gagaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGa                               >  2:2052227/1‑93 (MQ=255)

                                                          gagaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGa                               >  2:6315528/1‑93 (MQ=255)

                                                            gaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACg                             >  1:1000661/1‑93 (MQ=255)

                                                            gaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACg                             >  1:7827510/1‑93 (MQ=255)

                                                            gaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACg                             <  2:7832315/93‑1 (MQ=255)

                                                               cGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGc                          <  1:4828643/93‑1 (MQ=255)

                                                                gtgtTTCTTTTTCCCCGGATAAGCTCAACCAGGCTATTGACACCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCa                         <  2:155746/93‑1 (MQ=255)

                                                                gtgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCa                         >  2:7050282/1‑93 (MQ=255)

                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCaa                        <  1:2632816/93‑1 (MQ=255)

                                                                  gtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       >  2:4862118/1‑93 (MQ=255)

                                                                  gtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       >  1:3419105/1‑93 (MQ=255)

                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGTAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAcc                      <  1:2130105/93‑1 (MQ=255)

                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAcc                      <  2:7463877/93‑1 (MQ=255)

                                                                    aTCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCg                     >  2:994962/1‑93 (MQ=255)

                                                                     tCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCGACCGc                    >  1:3410456/1‑93 (MQ=255)

                                                                     tCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGc                    >  2:5329598/1‑93 (MQ=255)

                                                                          ttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGa               >  2:406665/1‑93 (MQ=255)

                                                                           tCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGaa              <  2:3905963/93‑1 (MQ=255)

                                                                            cGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAAc             <  2:6227453/93‑1 (MQ=255)

                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         >  2:6822215/1‑93 (MQ=255)

                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         <  2:1314174/93‑1 (MQ=255)

                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         <  1:5125072/93‑1 (MQ=255)

                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCt      <  1:4774654/93‑1 (MQ=255)

                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCt      >  1:1205195/1‑93 (MQ=255)

                                                                                    aaGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCtt     >  2:4852579/1‑93 (MQ=255)

                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    <  1:1103638/93‑1 (MQ=255)

                                                                                      gCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAg   <  2:2703774/93‑1 (MQ=255)

GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110

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