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.6%	925.7 / ‑5.4	254	D29N (GAT→AAT)	hemA	HemA
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (140/112); total (141/113)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 266 total aligned reads displayed.

GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110

 cAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc                                                                                        >  1:8490530/1‑93 (MQ=255)

  aggcATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  2:3036057/90‑1 (MQ=255)

   ggcATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  2:3795109/3‑93 (MQ=255)

   ggcATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  1:2095550/3‑93 (MQ=255)

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

    gcATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAgg                                                                                     <  2:1500086/92‑1 (MQ=255)

    gcATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAgg                                                                                     >  1:1569329/2‑93 (MQ=255)

    gcATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAgg                                                                                     <  2:3097346/92‑1 (MQ=255)

     cATGACCCTTTTAGCACTCGGTATCAACCGTAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    >  1:4127951/1‑93 (MQ=255)

     cATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    <  1:6604758/93‑1 (MQ=255)

     cATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    <  1:3186078/93‑1 (MQ=255)

     cATGACACTTTTAGCACTCGGTATCAACCATAAAACGGCAACTGTATCGCTGCGAGAACGTGTGTCGTTTTCGACTGATAAGATTAATCAgac                                                                                    >  2:8236642/1‑91 (MQ=255)

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

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

         aCCCTTTTGGCACGCGGGGTGATCCATAAAACTGCACCTGTATCGCTGTGAGGACGTGTATGGTTTTCGCCGGATAATCTCAATAAGgtgcgt                                                                                >  2:4822139/1‑91 (MQ=255)

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

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

           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTACCGTTTTCGCCTGATAAGCTCAATCAGGCGCTTGa                                                                              >  2:6861823/1‑93 (MQ=255)

            cTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGAc                                                                             <  1:6868245/93‑1 (MQ=255)

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

            cTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGAc                                                                             <  1:1260968/93‑1 (MQ=255)

            cTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGAc                                                                             <  1:9769975/93‑1 (MQ=255)

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

                 tgCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        <  2:7714736/92‑1 (MQ=255)

                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACCGCCt                                                                        >  2:2221141/1‑93 (MQ=255)

                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        >  1:2627431/1‑93 (MQ=255)

                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        <  1:6757815/93‑1 (MQ=255)

                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        >  1:8146971/1‑93 (MQ=255)

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

                  gCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTg                                                                       <  1:4980129/93‑1 (MQ=255)

                   cACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGc                                                                      <  2:8178191/93‑1 (MQ=255)

                    ccTCGGTATCAACCATAAAACGGCACCTGTATGGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     <  1:5978283/92‑1 (MQ=255)

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

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

                     cTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGATCAATCAGGCGCTTGACAGCCTGCtt                                                                    >  1:2381448/1‑93 (MQ=255)

                      tCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTg                                                                   <  2:5989549/93‑1 (MQ=255)

                       cGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCGGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgc                                                                  >  2:2118596/1‑93 (MQ=255)

                        ggTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgcg                                                                 <  1:5617622/93‑1 (MQ=255)

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

                          tATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCa                                                               >  1:9752757/1‑93 (MQ=255)

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

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

                            tCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGCATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGc                                                             >  1:4211543/1‑93 (MQ=255)

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

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

                                 cATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATg                                                        >  1:7568946/1‑93 (MQ=255)

                                  aTAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATgg                                                       >  1:1392870/1‑93 (MQ=255)

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

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

                                     aaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGc                                                    <  2:7505134/93‑1 (MQ=255)

                                      aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCACTCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   >  2:6270151/1‑93 (MQ=255)

                                       aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  <  1:9301941/93‑1 (MQ=255)

                                       aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCATGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  <  2:5770022/93‑1 (MQ=255)

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

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

                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTGGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 <  2:1784222/93‑1 (MQ=255)

                                            aCCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcgg                                             >  2:6199668/1‑93 (MQ=255)

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

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

                                                    cGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGc                                     <  1:10013485/93‑1 (MQ=255)

                                                     gCTGCGAGAACGTGTATCGTTTTCGCCCGATAAGCTCTATCAGACGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCt                                    >  2:1801885/1‑93 (MQ=255)

                                                      cTGCGAGAACGTTTATCGTTTTCGCAGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTg                                   <  2:5196901/93‑1 (MQ=255)

                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  <  1:4601049/93‑1 (MQ=255)

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

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

                                                              aagtgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTg                           <  2:6096552/91‑1 (MQ=255)

                                                              aCGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTg                           >  2:7345930/1‑93 (MQ=255)

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

                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCac                        >  1:3427702/1‑92 (MQ=255)

                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCaa                        >  1:381853/1‑93 (MQ=255)

                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCaa                        >  2:5464294/1‑93 (MQ=255)

                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGAGGTGCTGTCGACGTGCaa                        >  1:9090869/1‑93 (MQ=255)

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

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

                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   <  2:2618144/64‑1 (MQ=255)

                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAcc                      >  1:5412171/1‑93 (MQ=255)

                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCCCTTGCCAGCCTGCTTGCGCAGCCGATGGTGCAGCACGGCCTGGTGCTCTCGACGTGCAccc                      >  1:7176393/1‑93 (MQ=255)

                                                                     tCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCCACCGc                    >  1:5282099/1‑93 (MQ=255)

                                                                      cGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCa                   >  2:6951604/1‑93 (MQ=255)

                                                                       gTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCAc                  <  2:2681089/93‑1 (MQ=255)

                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  1:1943333/1‑93 (MQ=255)

                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCGTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  1:5293625/1‑93 (MQ=255)

                                                                             gCCGGATAGGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACt            <  2:8982978/93‑1 (MQ=255)

                                                                             gCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACt            >  1:3627166/1‑93 (MQ=255)

                                                                              ccGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACtt           <  2:1232260/93‑1 (MQ=255)

                                                                              ccGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTAGACGTGCAACCGCACGGAACtt           >  1:6335764/1‑93 (MQ=255)

                                                                               cGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGGGCttt          >  1:2656074/1‑93 (MQ=255)

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

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

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

                                                                                 gATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTAt        >  2:7626388/1‑93 (MQ=255)

                                                                                 gATAAGCTCAAGCAGGCGCTTGACACCCTTCTTGCCCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTAt        <  1:9563569/93‑1 (MQ=255)

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

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

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

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

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

                                                                                      gCTCAATCAGGCGCTTGACAGGCTGCTTGCGCAGTCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAg   <  1:6885911/93‑1 (MQ=255)

                                                                                       cTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGc  <  2:6323549/93‑1 (MQ=255)

                                                                                       cTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGc  >  1:8671494/1‑93 (MQ=255)

GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110

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