BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_012417	2,287	C→G	11.1%	E130Q (GAA→CAA)	USA300HOU_RS14900 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012417	2,287	0	C	G	11.1%	41068.2 / inf	12769	E130Q (GAA→CAA)	USA300HOU_RS14900	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (5759/5586); new base G (688/730); total (6449/6320)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 12929 total aligned reads displayed.

GCTACGCCCTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGTCGTT  >  NC_012417/2196‑2377

 cTACGCCCTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTGCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTc                                                                                            >  2:922014/1‑91 (MQ=255)

   aCGCCCTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAg                                                                                          <  2:2517155/91‑1 (MQ=255)

   aCGCCCTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAg                                                                                          >  1:2444944/1‑91 (MQ=255)

     gCCCTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGaa                                                                                        >  2:2748411/1‑91 (MQ=255)

        cTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAg                                                                                     <  2:1615690/91‑1 (MQ=255)

         tttAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGa                                                                                    >  2:1746424/1‑91 (MQ=255)

         tttAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATAAGAGTTTAAAACTTGAGAATAGa                                                                                    >  1:1901727/1‑91 (MQ=255)

           tAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGaaa                                                                                  <  1:1416516/91‑1 (MQ=255)

           tAAGTTTTTGGCTTTAACATCTTCAAAAACTACTGTGCTTATGTTTCCAGTTAACATTTCTCCACATGAGTTTAAAACTTCAGAATAGaaa                                                                                  >  2:3513015/1‑91 (MQ=255)

             AGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACa                                                                                <  2:1182262‑M1/91‑1 (MQ=255)

                TTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTc                                                                             <  2:2588188‑M1/91‑1 (MQ=255)

                 TTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCa                                                                            >  1:1222615‑M1/1‑91 (MQ=255)

                   tgaagatgttaaagccaaaaactGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCAcc                                                                          >  2:1036846‑M2/24‑91 (MQ=255)

                    gaagatgttaaagccaaaaactGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCg                                                                         >  1:408356‑M2/23‑91 (MQ=255)

                       tttAACATCTTCAAAAACTGCTGTGATTATGTTTCCCGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCtc                                                                      >  2:853464/1‑90 (MQ=255)

                         tAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTGAGAATAGAAACAGTCACCGCTAtt                                                                    <  1:2303633/91‑1 (MQ=255)

                           aCATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTAtttt                                                                  <  2:1685831/91‑1 (MQ=255)

                             aTCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTg                                                                >  2:237496/1‑91 (MQ=255)

                               cTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTa                                                              >  1:2144937/1‑91 (MQ=255)

                               cTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTa                                                              >  2:1824317/1‑91 (MQ=255)

                                ttCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAg                                                             >  2:1258363/1‑91 (MQ=255)

                                ttCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAg                                                             >  2:650166/1‑91 (MQ=255)

                                ttCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAg                                                             >  1:3567037/1‑91 (MQ=255)

                                  caaaaactGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGtt                                                           <  1:3369567‑M2/83‑1 (MQ=255)

                                    aaaactGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTcc                                                         >  1:1957946‑M2/7‑91 (MQ=255)

                                    aaaactGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTcc                                                         >  2:3312244‑M2/7‑91 (MQ=255)

                                      aactGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCcca                                                       >  2:3880994‑M2/5‑91 (MQ=255)

                                       actGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCccac                                                      <  2:2032443‑M2/88‑1 (MQ=255)

                                        ctGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCccacc                                                     <  2:2161020‑M2/89‑1 (MQ=255)

                                        ctGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCccacc                                                     >  1:1566590‑M2/3‑91 (MQ=255)

                                         tGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCccacca                                                    >  2:444027‑M2/2‑91 (MQ=255)

                                          gctgTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCaa                                                   >  1:1352017/1‑91 (MQ=255)

                                           ctgTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAAc                                                  <  1:971889/91‑1 (MQ=255)

                                           ctgTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAAc                                                  >  2:3241595/1‑91 (MQ=255)

                                           ctgTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAAc                                                  >  2:786710/1‑91 (MQ=255)

                                            tgtgATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACt                                                 >  2:3181756/1‑91 (MQ=255)

                                              tgATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAATAGTCACCGCTATTTTTGTAGTTCCCACCAACTGc                                               <  1:2649789/91‑1 (MQ=255)

                                              tgATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTGCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGc                                               <  1:3077401/91‑1 (MQ=255)

                                               gATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGct                                              >  2:1966396/1‑91 (MQ=255)

                                                aTTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGctc                                             <  1:3837201/91‑1 (MQ=255)

                                                  tATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCATCGCTATTTTTGTAGTTCCCACCAACTGCTCTa                                           >  1:3300620/1‑91 (MQ=255)

                                                  tATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTa                                           <  2:2890357/91‑1 (MQ=255)

                                                  tATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTCCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTa                                           >  2:3594603/1‑91 (MQ=255)

                                                   aTGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAg                                          <  2:2665226/91‑1 (MQ=255)

                                                   aTGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAAAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAg                                          <  1:1127381/91‑1 (MQ=255)

                                                     gTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTGAGAATAGAAACAGTCACCGGTATTTTTGTAGTTCCCACCAACTGCTCTAGGt                                        <  2:3811050/91‑1 (MQ=255)

                                                     gTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGt                                        <  1:2860641/91‑1 (MQ=255)

                                                     gTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGt                                        >  2:1401350/1‑91 (MQ=255)

                                                      tttCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTg                                       >  1:2794007/1‑91 (MQ=255)

                                                        tCCAGTTAAAAGTTCTCCATATGTGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATCTTTGTAGTTCCCACCAACTACTCTAGGTGat                                     <  2:2968671/91‑1 (MQ=255)

                                                        tCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGat                                     >  1:1707141/1‑91 (MQ=255)

                                                           aGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATAtt                                  <  2:3117116/91‑1 (MQ=255)

                                                            gTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTg                                 <  1:3773807/91‑1 (MQ=255)

                                                             ttAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGa                                >  2:3661077/1‑91 (MQ=255)

                                                              tAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGaa                               >  1:828079/1‑91 (MQ=255)

                                                               aaaaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAAc                              <  1:1000741/91‑1 (MQ=255)

                                                               aaaaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAAc                              >  1:3159719/1‑91 (MQ=255)

                                                               aaaaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAAc                              <  1:67967/91‑1 (MQ=255)

                                                               aaaaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAAc                              <  1:2721391/91‑1 (MQ=255)

                                                                aaaGTTCTCCATATGAGTTTAAAACTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACt                             >  1:2931128/1‑91 (MQ=255)

                                                                aaaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCACCTGCTCTAGGTGAAATTGAAcc                             >  2:2812963/1‑90 (MQ=255)

                                                                aaaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACt                             <  1:54062/91‑1 (MQ=255)

                                                                 aaGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTa                            <  1:3434306/91‑1 (MQ=255)

                                                                  aGTTCTCCATATGAGTTTAAAACTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGACCTAg                           >  1:3231203/1‑91 (MQ=255)

                                                                   gTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGt                          <  1:1057100/91‑1 (MQ=255)

                                                                   gTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGt                          >  1:1287765/1‑91 (MQ=255)

                                                                    ttCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGtt                         <  2:71527/91‑1 (MQ=255)

                                                                     tctcCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGttt                        >  1:1490441/1‑91 (MQ=255)

                                                                      ctcGATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACCAGtttt                       <  1:3011325/91‑1 (MQ=255)

                                                                      ctcCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGtttt                       <  2:3440663/91‑1 (MQ=255)

                                                                      ctcCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGtttt                       >  1:3636692/1‑91 (MQ=255)

                                                                      ctcCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGtttt                       <  2:245277/91‑1 (MQ=255)

                                                                       tcCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTc                      <  2:591335/91‑1 (MQ=255)

                                                                        ccATATGAGTTTAAAACTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCt                     >  1:1631681/1‑91 (MQ=255)

                                                                        ccATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCt                     >  1:2518091/1‑91 (MQ=255)

                                                                        ccATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCt                     <  1:615126/91‑1 (MQ=255)

                                                                         cATATGAGTTTAAAACTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCtt                    <  2:3042832/91‑1 (MQ=255)

                                                                         cATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCtt                    <  2:1894215/91‑1 (MQ=255)

                                                                          atatGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCttt                   <  2:1467743/91‑1 (MQ=255)

                                                                             tGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTa                <  2:2232723/91‑1 (MQ=255)

                                                                                 tttAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATtct            >  1:2372504/1‑91 (MQ=255)

                                                                                 tttAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATtct            <  1:3713923/91‑1 (MQ=255)

                                                                                  ttAAAACTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATtctc           <  1:2075284/91‑1 (MQ=255)

                                                                                  ttAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATtctc           >  2:2303803/1‑91 (MQ=255)

                                                                                  ttAAAACTTCAGAATAGAAACAGTCACCGATATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATtctc           >  2:3382640/1‑91 (MQ=255)

                                                                                   tAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTctct          >  1:2586716/1‑91 (MQ=255)

                                                                                   tAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTctct          <  1:1839513/91‑1 (MQ=255)

                                                                                     aaaCTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTAt        <  2:1328744/91‑1 (MQ=255)

                                                                                     aaaCTTGAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTAt        <  1:3505114/91‑1 (MQ=255)

                                                                                     aaaCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTAt        <  1:1778689/91‑1 (MQ=255)

                                                                                     aaaCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTAt        <  2:1534502/91‑1 (MQ=255)

                                                                                      aaCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATg       <  2:3738367/91‑1 (MQ=255)

                                                                                      aaCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATg       >  1:754152/1‑91 (MQ=255)

                                                                                       aCTTCAGTATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGt      <  1:2014233/91‑1 (MQ=255)

                                                                                       aCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCCACTGCTCTCGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGt      >  1:47583/1‑91 (MQ=255)

                                                                                       aCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGt      >  1:222245/1‑91 (MQ=255)

                                                                                       aCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGt      <  2:519559/91‑1 (MQ=255)

                                                                                       aCTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGt      <  1:888414/91‑1 (MQ=255)

                                                                                          tCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGTCGt   <  2:2101329/91‑1 (MQ=255)

                                                                                          tCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGTCGt   >  2:1111083/1‑91 (MQ=255)

GCTACGCCCTTTAAGTTTTTGGCTTTAACATCTTCAAAAACTGCTGTGATTATGTTTCCAGTTAAAAGTTCTCCATATGAGTTTAAAACTTCAGAATAGAAACAGTCACCGCTATTTTTGTAGTTCCCACCAACTGCTCTAGGTGATATTGAACTAGTTTTCTTTTTATTCTCTATGTCGTT  >  NC_012417/2196‑2377

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