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.1%	379.3 / ‑5.1	108	D29N (GAT→AAT)	hemA	HemA
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (66/41); total (67/41)
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.96e-01

Only 100 of 117 total aligned reads displayed.

CCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262929‑1263110

ccGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTGTCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCg                                                                                           <  2:3987726/93‑1 (MQ=255)

  gCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAAt                                                                                         >  1:2732484/1‑93 (MQ=255)

    aGACATGACCTTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  1:2767683/93‑1 (MQ=255)

    aGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCTTTTTCGCCGGATAAGCTCAATCa                                                                                       >  2:2603323/1‑93 (MQ=255)

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

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

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

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

     gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  1:4074395/1‑93 (MQ=255)

     gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  1:2369381/1‑93 (MQ=255)

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

        aTGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAggag                                                                                   >  1:1311586/1‑91 (MQ=255)

        aTGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgcg                                                                                   >  2:2995483/1‑93 (MQ=255)

          gACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCt                                                                                 <  2:3734211/93‑1 (MQ=255)

          gACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCt                                                                                 >  1:2295729/1‑93 (MQ=255)

          gACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCt                                                                                 <  1:3215652/93‑1 (MQ=255)

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

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

            cccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTg                                                                               >  1:4470498/1‑93 (MQ=255)

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

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

               ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGTGAACGTGTATCGTTTTCGCCGGATAATCTCAATCAGGCGCTTGata                                                                            >  2:2379524/1‑91 (MQ=255)

               ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGGACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            >  2:3202371/1‑93 (MQ=255)

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

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

                tttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAg                                                                           <  1:4185784/93‑1 (MQ=255)

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

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

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

                      ccTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     <  1:3488547/92‑1 (MQ=255)

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

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

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

                       cTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCtt                                                                    <  1:2716643/93‑1 (MQ=255)

                        tCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACCTCCTGCTTg                                                                   >  1:3027074/1‑93 (MQ=255)

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

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

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

                        tCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACACCCTGCTTg                                                                   >  1:1983924/1‑93 (MQ=255)

                         cGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgc                                                                  >  2:1766022/1‑93 (MQ=255)

                           gTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTgcgc                                                                >  2:3109435/1‑93 (MQ=255)

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

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

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

                              tCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGc                                                             <  1:649663/93‑1 (MQ=255)

                              tCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGc                                                             <  1:3068782/93‑1 (MQ=255)

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

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

                                 aCCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGa                                                          <  2:3003836/93‑1 (MQ=255)

                                 aCCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGa                                                          <  1:1308680/93‑1 (MQ=255)

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

                                     tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      <  2:3712094/93‑1 (MQ=255)

                                     aaaaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      <  2:1195343/92‑1 (MQ=255)

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

                                      aaaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTg                                                     <  2:4273629/93‑1 (MQ=255)

                                      aaaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTCCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTg                                                     <  1:828772/93‑1 (MQ=255)

                                        accGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   <  2:3082721/91‑1 (MQ=255)

                                        aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   >  2:4624993/1‑93 (MQ=255)

                                        aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCCTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   >  2:3567584/1‑93 (MQ=255)

                                         aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  <  1:1958132/93‑1 (MQ=255)

                                          cGGCACCTGTATCGTTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 <  2:1245882/93‑1 (MQ=255)

                                            gccccTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggc                                               <  1:689364/90‑1 (MQ=255)

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

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

                                              aCCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATCGTGCAGggcgg                                             >  2:2284186/1‑93 (MQ=255)

                                               ccTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcggc                                            >  2:1390179/1‑93 (MQ=255)

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

                                                 tGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgt                                          <  1:4508819/93‑1 (MQ=255)

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

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

                                                   tATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtgg                                        >  1:3365061/1‑93 (MQ=255)

                                                    accGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGACGATGGTGCAGGGCGGCgtggt                                       <  1:3807156/91‑1 (MQ=255)

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

                                                        cTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTg                                   <  2:2354612/93‑1 (MQ=255)

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

                                                          gCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGTCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTc                                 >  1:3081536/1‑93 (MQ=255)

                                                             agaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGAc                              <  1:289456/93‑1 (MQ=255)

                                                             agaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGAc                              >  1:131590/1‑93 (MQ=255)

                                                             agaACGGGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGTCGTCGag                              >  2:1011394/1‑92 (MQ=255)

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

                                                               aaCGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGt                            >  1:3406471/1‑93 (MQ=255)

                                                                    gtATCTTTGTCGCCGGATAAGCTCAATCCGGCGCTTGACAGCGTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       <  1:2498728/93‑1 (MQ=255)

                                                                    gtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       <  2:893085/93‑1 (MQ=255)

                                                                    gtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       <  1:2498607/93‑1 (MQ=255)

                                                                         gTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCAc                  >  1:3641831/1‑93 (MQ=255)

                                                                           tttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACgg                >  2:2116363/1‑93 (MQ=255)

                                                                             tCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGaa              >  2:4298908/1‑93 (MQ=255)

                                                                              cGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAAc             >  1:1127529/1‑93 (MQ=255)

                                                                               cccGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACt            <  2:4328776/92‑1 (MQ=255)

                                                                                 cGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTTTCGACGTGCAACCGCACGGAACttt          >  1:3285054/1‑93 (MQ=255)

                                                                                 cGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACttt          >  2:4487796/1‑93 (MQ=255)

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

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

                                                                                  ggATAAGCTCAATCAGGAGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         >  2:2956640/1‑93 (MQ=255)

                                                                                    aTAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATc       >  1:4551893/1‑93 (MQ=255)

                                                                                      aaGCTCAATGAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGGAACCGCACGGAACTTTATCtt     <  2:1508488/93‑1 (MQ=255)

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

                                                                                        gCTCAATCAGGCGCTTAACGGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAg   <  2:3612225/93‑1 (MQ=255)

                                                                                         cTCAATCAGGCGCTTGATAGCCTGCTTGCACAGCCGATGGTGCAGGACGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGc  <  2:2454806/93‑1 (MQ=255)

                                                                                         cTCAATCAGGCGCTTGACAGCGTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGc  <  1:2497836/93‑1 (MQ=255)

CCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262929‑1263110

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