BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,596,699	T→C	Q498R (CAG→CGG)	ydeV ←	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,596,699	0	T	C	99.4%	525.2 / ‑4.7	162	Q498R (CAG→CGG)	ydeV	putative kinase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (76/85); total (76/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 169 total aligned reads displayed.

CCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTCCCAGCGA  >  NC_000913/1596609‑1596788

ccACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGa                                                                                         <  1:18115/93‑1 (MQ=255)

 cACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGAt                                                                                        <  1:1656605/93‑1 (MQ=255)

  aCGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATc                                                                                       >  1:2200094/1‑93 (MQ=255)

   cGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATcc                                                                                      >  1:2030811/1‑93 (MQ=255)

    ggTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCt                                                                                     <  1:3750515/93‑1 (MQ=255)

      tcatcaTGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGa                                                                                   >  2:1309036/1‑93 (MQ=255)

      tcatcaTGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGa                                                                                   >  2:2993617/1‑93 (MQ=255)

        atcatGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATa                                                                                 <  1:1924165/93‑1 (MQ=255)

          catGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATaaa                                                                               <  1:733726/93‑1 (MQ=255)

            tgatgTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAAct                                                                             <  1:1479751/93‑1 (MQ=255)

             gatgTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAActg                                                                            <  1:3246960/93‑1 (MQ=255)

             gatgTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAActg                                                                            >  1:2392950/1‑93 (MQ=255)

               tgtgGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAActgcc                                                                          >  2:1932242/1‑93 (MQ=255)

               tggggCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAActgcc                                                                          <  2:1065088/90‑1 (MQ=255)

                  ggCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAActgcctgc                                                                       >  1:4577393/1‑93 (MQ=255)

                  ggCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAActgcctgc                                                                       >  1:3163203/1‑93 (MQ=255)

                   gCTGTCAATGAAACTATACCCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACtgcctgcc                                                                      <  1:1875957/93‑1 (MQ=255)

                     tGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCAc                                                                    >  2:3512523/1‑93 (MQ=255)

                     tGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCAc                                                                    >  2:2678051/1‑93 (MQ=255)

                     tGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCAc                                                                    <  2:4159169/93‑1 (MQ=255)

                      gTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCCAGCCGATCCTGATAAACTGCCTGCCAc                                                                    <  2:1433513/93‑1 (MQ=255)

                       tCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACtt                                                                  >  2:4381305/1‑93 (MQ=255)

                        cAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTa                                                                 <  1:897726/93‑1 (MQ=255)

                        cAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTa                                                                 <  2:3407391/93‑1 (MQ=255)

                         aaTGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTAt                                                                >  1:2454985/1‑93 (MQ=255)

                         aaTGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTAt                                                                >  2:3885063/1‑93 (MQ=255)

                          aTGAAAAAAAAACCCAGGCACTTTCCATAACGACGTCGTCAGTCATTGATCATCCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATc                                                               <  2:3432406/93‑1 (MQ=255)

                            gAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATcgc                                                             <  1:4309913/93‑1 (MQ=255)

                             aaaCTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATcgcg                                                            >  1:2704342/1‑93 (MQ=255)

                               aCTATAACCCAGGCGTTTTCCATAACGACGTCGTCAGCCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                                          <  2:1669463/93‑1 (MQ=255)

                               aCTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                                          <  2:3478337/93‑1 (MQ=255)

                                cTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGa                                                         >  1:3440203/1‑93 (MQ=255)

                                 tataACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                                        >  2:2115773/1‑93 (MQ=255)

                                 tataACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                                        <  1:4745332/93‑1 (MQ=255)

                                 tataACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTCCCTGCCCCTTCTCGCGTGaa                                                        >  1:1047463/1‑93 (MQ=255)

                                  ataACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGaac                                                       >  1:848034/1‑92 (MQ=255)

                                  ataACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAAt                                                       <  1:3360101/93‑1 (MQ=255)

                                  ataACCCAGGCGATTTCGATAACCACGTCGTCAGTCCACGATCAACCAG‑CTCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAAt                                                       >  2:2492220/1‑93 (MQ=255)

                                   taACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATc                                                      >  2:22908/1‑93 (MQ=255)

                                   taACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATc                                                      >  1:4063460/1‑93 (MQ=255)

                                    aaCCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATcc                                                     >  1:3929203/1‑93 (MQ=255)

                                      cccAGGCGCTTTCCATAACGACGCCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTg                                                   >  2:1133403/1‑93 (MQ=255)

                                       ccAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGa                                                  <  1:1155826/93‑1 (MQ=255)

                                       ccAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGa                                                  <  2:3656451/93‑1 (MQ=255)

                                       ccAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGa                                                  <  2:3058163/93‑1 (MQ=255)

                                         aGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCTCGTCAATCCTGATa                                                >  2:2441127/1‑93 (MQ=255)

                                          ggCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaa                                               <  2:583401/93‑1 (MQ=255)

                                          ggCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaa                                               >  1:4684741/1‑93 (MQ=255)

                                          ggCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaa                                               >  2:3837065/1‑93 (MQ=255)

                                           gcgcTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaaa                                              <  1:2066176/93‑1 (MQ=255)

                                             gcTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGt                                            >  1:2303761/1‑93 (MQ=255)

                                              cTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGtt                                           <  1:1415497/93‑1 (MQ=255)

                                              cTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGtt                                           >  1:3770494/1‑93 (MQ=255)

                                              cTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGtt                                           >  2:1262219/1‑93 (MQ=255)

                                               tttCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTc                                          <  2:1434206/93‑1 (MQ=255)

                                                ttCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCa                                         >  2:3111174/1‑93 (MQ=255)

                                                ttCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCa                                         >  1:2994025/1‑93 (MQ=255)

                                                 tcaaTAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                                        <  1:2026795/90‑1 (MQ=255)

                                                 tCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                                        <  2:4702219/93‑1 (MQ=255)

                                                 tCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                                        >  1:3548140/1‑93 (MQ=255)

                                                 tCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                                        <  1:3147580/93‑1 (MQ=255)

                                                  ccATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATg                                       >  2:1035769/1‑93 (MQ=255)

                                                    aTAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCt                                     <  2:4218289/93‑1 (MQ=255)

                                                     tAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCtt                                    >  2:1962358/1‑93 (MQ=255)

                                                      aaCGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCttt                                   <  2:4138747/93‑1 (MQ=255)

                                                      aaCGACGTCGTCAGCCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCttt                                   <  1:2798564/93‑1 (MQ=255)

                                                       acgacgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCtttt                                  <  2:3468747/93‑1 (MQ=255)

                                                        ctacgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                                 <  1:585964/91‑1 (MQ=255)

                                                        cgccgtcgtCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                                 <  1:1235082/90‑1 (MQ=255)

                                                        cgacgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                                 >  2:3501679/1‑93 (MQ=255)

                                                        cgacgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                                 <  2:1107531/93‑1 (MQ=255)

                                                         gacgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTcc                                >  1:980254/1‑93 (MQ=255)

                                                         gacgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTcc                                >  2:4278267/1‑93 (MQ=255)

                                                          acgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCg                               <  2:331430/93‑1 (MQ=255)

                                                          acgTCGTCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCg                               <  1:3497168/93‑1 (MQ=255)

                                                           cgtcgtCAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCgg                              <  1:3057687/93‑1 (MQ=255)

                                                            gtcgtcAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCggg                             <  2:757040/93‑1 (MQ=255)

                                                               gtcAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCt                          >  1:1045339/1‑93 (MQ=255)

                                                               gtcAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCt                          <  1:4452884/93‑1 (MQ=255)

                                                                tgaGTCCAAGAACAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg                         <  1:2929678/91‑1 (MQ=255)

                                                                tcAGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg                         >  1:2855208/1‑93 (MQ=255)

                                                                  aGTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgt                       <  2:1369886/93‑1 (MQ=255)

                                                                   gTCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtg                      >  1:1876258/1‑93 (MQ=255)

                                                                    tCCATGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgt                     >  1:36856/1‑93 (MQ=255)

                                                                     tcATGATCAACCAG‑CCCCATCCGATCCTGATAAACTGCCTGCCACTTATCGAGTGAATCCTGATAAAGGGCATGCTGTTCCGGGTCTGgtgtg                    >  1:4196376/2‑93 (MQ=255)

                                                                        tGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGc                 <  1:3114576/93‑1 (MQ=255)

                                                                        tGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGc                 <  2:3372763/93‑1 (MQ=255)

                                                                        tGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGc                 >  1:2929492/1‑93 (MQ=255)

                                                                        tGATCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGc                 >  1:2519284/1‑93 (MQ=255)

                                                                          aTCAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGt               <  1:1182559/93‑1 (MQ=255)

                                                                            cAACCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTcc             <  2:4721483/93‑1 (MQ=255)

                                                                             aaCCAG‑CCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCg            <  2:2803348/93‑1 (MQ=255)

                                                                             aaCCAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCg            >  2:441327/1‑93 (MQ=255)

                                                                                cAG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATACTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTc         <  1:2933828/93‑1 (MQ=255)

                                                                                 aG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTcc        <  1:4249933/93‑1 (MQ=255)

                                                                                 aG‑CCCCAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTcc        <  1:2697863/93‑1 (MQ=255)

                                                                                     cccAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTCCCAg     <  2:3350879/93‑1 (MQ=255)

                                                                                      ccAGCCGATTCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTTATGCTCTTCCGGGTCTGGTGTGTGCGTCCGCGCCCAGc    >  2:3329873/1‑93 (MQ=255)

                                                                                      ccAGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTCCCAGc    <  2:2719729/93‑1 (MQ=255)

                                                                                        aGCCGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTCCCAGCGa  >  2:3151935/1‑93 (MQ=255)

CCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAG‑CCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGTCCGTTCCCAGCGA  >  NC_000913/1596609‑1596788

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