BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,041,340	A→G	A698A (GCT→GCC)	yccC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,041,340	0	A	G	99.3%	436.3 / ‑6.2	138	A698A (GCT→GCC)	yccC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (73/64); total (73/65)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAG  >  NC_000913/1041250‑1041429

gTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGgc                                                                                         <  1:293046/93‑1 (MQ=255)

gTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGgc                                                                                         >  2:3457640/1‑93 (MQ=255)

 ttttACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGgcg                                                                                        >  2:839603/1‑93 (MQ=255)

 ttttACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGgcg                                                                                        >  1:3503974/1‑93 (MQ=255)

 ttttACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGgcg                                                                                        >  2:2437574/1‑93 (MQ=255)

    tACTCTTTCTCGGAGTAACTATAATCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCt                                                                                     <  1:265856/93‑1 (MQ=255)

      ctctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCttt                                                                                   >  2:4296900/1‑93 (MQ=255)

       tctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTa                                                                                  >  2:2034442/1‑93 (MQ=255)

       tctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTCGCGCGTTTGATCACACCATTGAGGATGGCGCCTTTa                                                                                  >  2:1765952/1‑93 (MQ=255)

        ctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTaa                                                                                 >  1:488734/1‑93 (MQ=255)

        ctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTaa                                                                                 <  1:4052548/93‑1 (MQ=255)

           tctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAatat                                                                              <  2:1776203/93‑1 (MQ=255)

           tctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAatat                                                                              >  1:715421/1‑93 (MQ=255)

            ctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATAtt                                                                             <  1:1460219/93‑1 (MQ=255)

            ctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATAtt                                                                             <  1:4411528/93‑1 (MQ=255)

            ctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATAtt                                                                             <  1:1083970/93‑1 (MQ=255)

             tcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGGCTTTAATATTg                                                                            >  1:431632/1‑93 (MQ=255)

                gAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACg                                                                         <  2:4601305/93‑1 (MQ=255)

                 aGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgc                                                                        <  2:2347658/93‑1 (MQ=255)

                   tAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgcct                                                                      >  2:369605/1‑93 (MQ=255)

                    aaCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgcctg                                                                     >  1:128018/1‑93 (MQ=255)

                     aCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgcctgc                                                                    >  1:4681817/1‑93 (MQ=255)

                     aCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgcctgc                                                                    <  1:3517086/93‑1 (MQ=255)

                      cTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgcctgcc                                                                   <  2:534301/93‑1 (MQ=255)

                      cTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACgcctgcc                                                                   <  2:2814923/93‑1 (MQ=255)

                        ataACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCGTTTGATATTGACGcctgcctg                                                                 >  1:1944644/1‑93 (MQ=255)

                        ataACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGcctgcctg                                                                 >  1:995807/1‑93 (MQ=255)

                        ataACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGcctgcctg                                                                 >  2:1469587/1‑93 (MQ=255)

                         taACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGt                                                                >  1:278955/1‑93 (MQ=255)

                         taACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGt                                                                >  1:3671553/1‑93 (MQ=255)

                           aCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTCATATTGACGCCTGCCTGTTc                                                              >  1:4349862/1‑93 (MQ=255)

                           aCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTc                                                              >  1:2825611/1‑93 (MQ=255)

                            ccGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTcc                                                             >  2:3547451/1‑93 (MQ=255)

                            ccGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTcc                                                             <  2:1252374/93‑1 (MQ=255)

                             cGTAATAGTTATAGCCGTAACTGTAAGCGGTTCTGGCGCGTTTAATCACACCATTGAGGAGGGCGCCTTTAATATTAACGCCTGCCTGTTCCa                                                            >  2:2326272/1‑93 (MQ=255)

                             cGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCa                                                            <  1:831206/93‑1 (MQ=255)

                             cGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCa                                                            <  1:4129053/93‑1 (MQ=255)

                              gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTcctg                                                           >  1:1318867/1‑91 (MQ=255)

                              gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAg                                                           >  1:3467990/1‑93 (MQ=255)

                              gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAg                                                           <  1:1989083/93‑1 (MQ=255)

                               tAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                          >  2:1247798/1‑93 (MQ=255)

                               tAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                          <  1:61385/93‑1 (MQ=255)

                                aaTAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGAc                                                         >  1:1361699/1‑93 (MQ=255)

                                  tAGTTATAGCCGTAGCTGTAAGCGGTGCTGGCGCGTTTAATCACGCCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGc                                                       <  2:1612374/93‑1 (MQ=255)

                                  tAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGc                                                       >  1:3268071/1‑93 (MQ=255)

                                   aGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCt                                                      >  1:1363492/1‑93 (MQ=255)

                                   aGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCt                                                      <  2:2965065/93‑1 (MQ=255)

                                   aGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCt                                                      <  1:272465/93‑1 (MQ=255)

                                     ttATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGc                                                    <  1:2483137/93‑1 (MQ=255)

                                     ttATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGc                                                    >  1:878321/1‑93 (MQ=255)

                                     ttATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGc                                                    <  2:3200903/93‑1 (MQ=255)

                                      tataGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCa                                                   >  1:1503020/1‑93 (MQ=255)

                                       ataGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAt                                                  <  2:439275/93‑1 (MQ=255)

                                        taGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                 >  2:969607/1‑93 (MQ=255)

                                        taGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                 <  2:4499953/93‑1 (MQ=255)

                                          gCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGa                                               >  1:1172260/1‑93 (MQ=255)

                                            cGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGAAGGCGCCTTTTATATTGACGCCTGCCCGTTCCAGTTGCTTCATTGACa                                             >  2:1699163/1‑93 (MQ=255)

                                              tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACaaa                                           >  1:791544/1‑93 (MQ=255)

                                              tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACaaa                                           >  2:3967637/1‑93 (MQ=255)

                                              tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGAAGCTGCATTGAaaaa                                           >  1:992775/1‑93 (MQ=255)

                                               aaCTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAAc                                          >  1:2026900/1‑93 (MQ=255)

                                                aCTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACt                                         >  1:2584155/1‑93 (MQ=255)

                                                aCTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACt                                         >  1:3911922/1‑93 (MQ=255)

                                                  tGTAAGCGGTGCTGGCGCGTTTAATCACAGCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCa                                       >  2:2468630/1‑93 (MQ=255)

                                                   gTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCAc                                      >  1:4164844/1‑93 (MQ=255)

                                                    tAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCAcc                                     <  2:2663386/93‑1 (MQ=255)

                                                    tAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCAcc                                     <  1:1231308/93‑1 (MQ=255)

                                                      aGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACctc                                   <  2:704844/93‑1 (MQ=255)

                                                      aGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACctc                                   <  1:3527981/93‑1 (MQ=255)

                                                       gCGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACctct                                  >  1:4411279/1‑93 (MQ=255)

                                                        cGGTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCtt                                 <  1:228195/93‑1 (MQ=255)

                                                         ggTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCttt                                <  1:3211281/93‑1 (MQ=255)

                                                          gTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTg                               <  1:3603213/93‑1 (MQ=255)

                                                          gTGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTg                               >  2:1861098/1‑93 (MQ=255)

                                                           tGCTGGCGCGTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTgg                              <  1:493661/93‑1 (MQ=255)

                                                                 cgcgTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtgt                        <  1:1289024/93‑1 (MQ=255)

                                                                  gcgTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGtt                       <  1:2244863/93‑1 (MQ=255)

                                                                  gcgTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGtt                       >  1:4160408/1‑93 (MQ=255)

                                                                    gTTTAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCa                     <  1:4117344/93‑1 (MQ=255)

                                                                     tttAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCaa                    >  1:1019318/1‑93 (MQ=255)

                                                                     tttAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCaa                    >  1:3014130/1‑93 (MQ=255)

                                                                     tttAATCACACCATTGAGGATGGCGCCCTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGGGTTCaa                    >  1:1845028/1‑93 (MQ=255)

                                                                     ttaaaTCACACCATTGAGGATGGCGCCTTTAATATGGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCaa                    <  2:972316/90‑1 (MQ=255)

                                                                       tAATCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGc                  >  2:4127564/1‑93 (MQ=255)

                                                                        aaTCACACCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGcc                 >  2:3144030/1‑93 (MQ=255)

                                                                           cacaCCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaa              >  1:2527135/1‑93 (MQ=255)

                                                                           cacaCCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaa              >  2:4328747/1‑93 (MQ=255)

                                                                           cacaCCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaa              <  2:2563322/93‑1 (MQ=255)

                                                                            acacCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaaa             >  2:1352432/1‑93 (MQ=255)

                                                                            acacCATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaaa             >  1:787255/1‑93 (MQ=255)

                                                                               ccATTGAGGATGGCGCTTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc          <  2:549539/93‑1 (MQ=255)

                                                                               ccATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc          >  1:2796543/1‑93 (MQ=255)

                                                                               ccATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc          >  1:229552/1‑93 (MQ=255)

                                                                               ccATTGAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCGCTTTGGCGGTGTTGAAGCCAAAAcgc          >  1:38418/1‑93 (MQ=255)

                                                                                   tGAGGAGGGCCCCTTTAATTTTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCaa      <  1:2568935/93‑1 (MQ=255)

                                                                                    gAGGATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAAc     <  2:3742181/93‑1 (MQ=255)

                                                                                      ggATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACca   >  1:4073834/1‑93 (MQ=255)

                                                                                      ggATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACca   >  1:1792997/1‑93 (MQ=255)

                                                                                      ggATGGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACca   >  2:2564686/1‑93 (MQ=255)

                                                                                       gATGGCGCCTTTAATATTGCCGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcag  <  2:4293537/93‑1 (MQ=255)

GTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAG  >  NC_000913/1041250‑1041429

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