BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,740,514	T→C	R51R (CGT→CGC)	yiaK →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,740,514	0	T	C	99.8%	1298.5 / NA	407	R51R (CGT→CGC)	yiaK	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (201/205); total (202/205)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 419 total aligned reads displayed.

TGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGC  >  NC_000913/3740424‑3740603

tGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGctt                                                                                         <  1:1029637/93‑1 (MQ=255)

 gACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttc                                                                                        <  2:23933/93‑1 (MQ=255)

 gACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttc                                                                                        >  2:3639743/1‑93 (MQ=255)

  aCAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttca                                                                                       <  1:5713361/93‑1 (MQ=255)

  aCAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttca                                                                                       <  1:393359/93‑1 (MQ=255)

   cAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcat                                                                                      <  1:5319376/93‑1 (MQ=255)

    aGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcatt                                                                                     <  2:1953855/93‑1 (MQ=255)

     gCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcattc                                                                                    >  1:7551283/1‑93 (MQ=255)

     gCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcattc                                                                                    <  2:4531055/93‑1 (MQ=255)

     gCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttcattc                                                                                    <  1:306063/93‑1 (MQ=255)

       gAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaa                                                                                  >  2:7902706/1‑93 (MQ=255)

         aaCGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaaca                                                                                >  1:4819516/1‑93 (MQ=255)

           cGGCTGACTCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCaacaac                                                                              <  1:2305143/93‑1 (MQ=255)

            ggCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACt                                                                             >  1:4151639/1‑93 (MQ=255)

             gCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTg                                                                            <  1:4317837/93‑1 (MQ=255)

              cTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTTCCTCGCTTCATTCAACAACTgg                                                                           >  2:1066448/1‑93 (MQ=255)

               tGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGa                                                                          >  1:1518119/1‑93 (MQ=255)

               agACGCCTGTGGAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGa                                                                          <  2:3287751/92‑1 (MQ=255)

                gACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTggca                                                                         >  1:6108513/1‑91 (MQ=255)

                gACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGaa                                                                         >  2:6717109/1‑93 (MQ=255)

                  cGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGaaaa                                                                       <  1:7913610/93‑1 (MQ=255)

                  cGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGaaaa                                                                       <  1:5836271/93‑1 (MQ=255)

                   gCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAAc                                                                      <  2:4401290/93‑1 (MQ=255)

                   gCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAAc                                                                      >  2:2013260/1‑93 (MQ=255)

                   gCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAAc                                                                      <  1:7326115/93‑1 (MQ=255)

                    ccTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACg                                                                     <  1:5063034/93‑1 (MQ=255)

                    ccTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACg                                                                     >  1:4115704/1‑93 (MQ=255)

                      tgtgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGc                                                                   >  1:7090945/1‑93 (MQ=255)

                      tgtgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGc                                                                   <  1:1793787/93‑1 (MQ=255)

                        tgCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGa                                                                 <  1:3226241/93‑1 (MQ=255)

                         gCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGat                                                                >  2:1197994/1‑93 (MQ=255)

                         gCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGat                                                                <  2:5592507/93‑1 (MQ=255)

                          cAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGata                                                               <  1:6930381/93‑1 (MQ=255)

                           agagATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGatat                                                              <  2:3713057/93‑1 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             >  2:2074521/1‑93 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             <  2:6616300/93‑1 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             <  1:4631403/93‑1 (MQ=255)

                             agaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCa                                                            <  2:4569713/93‑1 (MQ=255)

                             agaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCa                                                            <  2:726294/93‑1 (MQ=255)

                              gaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAt                                                           <  1:2857763/93‑1 (MQ=255)

                               aTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCCACAACTGGAAAACGGCGATATCAtt                                                          <  1:6032783/93‑1 (MQ=255)

                               aTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAtt                                                          <  1:4422677/93‑1 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         >  1:4507187/1‑93 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         >  2:5298373/1‑93 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         >  1:6738942/1‑93 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         >  2:6187692/1‑93 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         <  1:3288949/93‑1 (MQ=255)

                                 gTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTcc                                                        <  1:2519186/93‑1 (MQ=255)

                                  ttCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCt                                                       >  1:3610581/1‑93 (MQ=255)

                                    cGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGa                                                     >  2:562403/1‑93 (MQ=255)

                                    cGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGa                                                     <  1:5174713/93‑1 (MQ=255)

                                    cGCCCGCACCACCGAATCCGGCGATTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGa                                                     <  2:5661636/93‑1 (MQ=255)

                                     gCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGAt                                                    <  1:3767804/93‑1 (MQ=255)

                                         gcgccaccGAATCCGGCGTTTATTCTCACGGCATTAATCGTTTCCCTCGCTTCATTCAACAAATGGAAAACGGCGATATCATTCCTGATGccc                                                <  1:3439827/90‑1 (MQ=255)

                                         gCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccc                                                >  1:126073/1‑93 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            >  2:375797/1‑93 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            <  1:5231671/93‑1 (MQ=255)

                                              accGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaaccc                                           <  2:3980647/93‑1 (MQ=255)

                                                cGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacccaa                                         >  1:5896642/1‑93 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        <  2:1812000/93‑1 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        >  1:2485620/1‑93 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        >  2:6518231/1‑93 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        <  1:5465201/93‑1 (MQ=255)

                                                   aTCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACg                                      >  2:745789/1‑93 (MQ=255)

                                                    tCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGt                                     <  1:3489542/93‑1 (MQ=255)

                                                     ccGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGta                                    <  1:949280/93‑1 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   <  2:4994408/93‑1 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   <  2:4703397/93‑1 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   >  2:2570251/1‑93 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   >  2:5856950/1‑93 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  >  1:5502110/1‑93 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  >  1:1581197/1‑93 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  <  2:2785474/93‑1 (MQ=255)

                                                        gCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATaa                                 <  2:5086065/93‑1 (MQ=255)

                                                           tttATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCa                              >  1:6322593/1‑93 (MQ=255)

                                                             tATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGc                            <  1:4969776/93‑1 (MQ=255)

                                                                tctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTc                         <  2:2124671/93‑1 (MQ=255)

                                                                 ctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATCCCTGATGCCCAACCCAAACATATAACCAGCCTCg                        <  2:7169158/93‑1 (MQ=255)

                                                                    aCGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgcg                     >  1:6241197/1‑93 (MQ=255)

                                                                    aCGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgcg                     >  1:2077618/1‑93 (MQ=255)

                                                                       gCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCaa                  <  2:1462870/93‑1 (MQ=255)

                                                                        cGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAt                 >  2:4186126/1‑93 (MQ=255)

                                                                         gTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt                <  2:2946076/93‑1 (MQ=255)

                                                                         gTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt                <  1:6635910/93‑1 (MQ=255)

                                                                         gTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt                >  1:2392314/1‑93 (MQ=255)

                                                                          ttAATCGTCTCCCTCGCTTCATTGAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTg               <  2:5394287/93‑1 (MQ=255)

                                                                           tAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGa              <  2:338899/93‑1 (MQ=255)

                                                                            aaTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGaa             <  2:6868587/93‑1 (MQ=255)

                                                                            aaTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGaa             <  2:7578398/93‑1 (MQ=255)

                                                                             aTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAAc            >  1:4756771/1‑93 (MQ=255)

                                                                              tCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACa           >  1:774691/1‑93 (MQ=255)

                                                                              tCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACa           >  2:3498161/1‑93 (MQ=255)

                                                                              acGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACa           <  1:7645678/92‑1 (MQ=255)

                                                                                gTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGt         >  2:4780286/1‑93 (MQ=255)

                                                                                 tttCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTg        <  2:6004731/93‑1 (MQ=255)

                                                                                     ccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGag    >  1:2658362/1‑92 (MQ=255)

                                                                                     ccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGAc    <  2:2673642/93‑1 (MQ=255)

                                                                                      cTCGCTTCGTTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACg   >  1:486104/1‑93 (MQ=255)

                                                                                      cTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACg   >  2:3558716/1‑93 (MQ=255)

                                                                                       tCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGc  >  2:166627/1‑93 (MQ=255)

TGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGC  >  NC_000913/3740424‑3740603

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