BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,795,756	A→G	A394A (GCA→GCG)	rfaL →	O‑antigen ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,795,756	0	A	G	87.6%	669.7 / 83.1	267	A394A (GCA→GCG)	rfaL	O‑antigen ligase
Reads supporting (aligned to +/- strand): ref base A (18/15); new base G (106/128); total (124/143)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 273 total aligned reads displayed.

TTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGT  >  NC_000913/3795664‑3795846

ttttATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCa                                                                                            >  2:3975294/1‑93 (MQ=255)

  ttATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGcgcg                                                                                          <  1:230583/93‑1 (MQ=255)

   tatataGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGgcgcgc                                                                                         >  2:2532829/1‑93 (MQ=255)

   tatataGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACgc                                                                                         >  1:4907990/1‑93 (MQ=255)

    atataGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCgca                                                                                        <  1:6044852/93‑1 (MQ=255)

     tataGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCgcag                                                                                       >  2:2465301/1‑93 (MQ=255)

     tataGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCgcag                                                                                       >  1:2095602/1‑93 (MQ=255)

      ataGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCgcagc                                                                                      <  1:2953408/93‑1 (MQ=255)

       taGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCgcagca                                                                                     >  1:1123583/1‑93 (MQ=255)

       taGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACgcagca                                                                                     >  2:2028278/1‑93 (MQ=255)

        aGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCAt                                                                                    <  2:5129221/93‑1 (MQ=255)

           atatAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTcc                                                                                 <  2:945503/93‑1 (MQ=255)

           atatAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTcc                                                                                 >  2:7033367/1‑93 (MQ=255)

           atatAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTcc                                                                                 >  2:6564215/1‑93 (MQ=255)

            tataAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCa                                                                                >  2:4514860/1‑93 (MQ=255)

            tataAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCa                                                                                >  1:3419545/1‑93 (MQ=255)

            tataAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCa                                                                                >  1:6317474/1‑93 (MQ=255)

            tataAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCa                                                                                >  1:2592150/1‑93 (MQ=255)

              taAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAAt                                                                              >  2:627436/1‑93 (MQ=255)

              taAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAAt                                                                              <  2:5012470/93‑1 (MQ=255)

              taAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAAt                                                                              <  2:1868354/93‑1 (MQ=255)

                aaaaaCGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTa                                                                            >  1:2225296/1‑93 (MQ=255)

                 aaaaCGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTAt                                                                           >  1:6890237/1‑93 (MQ=255)

                  aaaCGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATc                                                                          >  1:1039033/1‑93 (MQ=255)

                  aaaCGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATc                                                                          >  1:1887460/1‑93 (MQ=255)

                   aaCGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCa                                                                         >  2:1297958/1‑93 (MQ=255)

                   aaCGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCa                                                                         >  1:4229359/1‑93 (MQ=255)

                     cGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCAtt                                                                       <  1:5288763/93‑1 (MQ=255)

                      gAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATta                                                                      >  2:3582518/1‑93 (MQ=255)

                        gCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATtata                                                                    <  1:171649/93‑1 (MQ=255)

                        gCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATtata                                                                    <  1:411724/93‑1 (MQ=255)

                        gCTTTGGGTTTGTTCATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATtata                                                                    >  2:6454752/1‑93 (MQ=255)

                         cTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTatat                                                                   <  1:368713/93‑1 (MQ=255)

                         cTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTatat                                                                   <  2:1713129/93‑1 (MQ=255)

                          tttGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATc                                                                  <  2:4141226/93‑1 (MQ=255)

                          tttGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATc                                                                  <  2:5622198/93‑1 (MQ=255)

                             gggTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCTGCATTCCAATTATCATTATATCCGc                                                               <  2:6026882/93‑1 (MQ=255)

                              ggTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCt                                                              >  1:3287176/1‑93 (MQ=255)

                               gTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCta                                                             <  1:4638672/93‑1 (MQ=255)

                                tttGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCtat                                                            <  2:1540281/93‑1 (MQ=255)

                                 ttgttgATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCtata                                                           >  2:4916257/1‑93 (MQ=255)

                                   gttgATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGt                                                         <  1:884066/93‑1 (MQ=255)

                                   gttgATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGt                                                         <  1:7061367/93‑1 (MQ=255)

                                   gttgATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGt                                                         <  1:2421228/93‑1 (MQ=255)

                                      gATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCct                                                      <  1:2805991/93‑1 (MQ=255)

                                       atatTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCctc                                                     <  1:1476782/93‑1 (MQ=255)

                                       atatTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCctc                                                     <  1:563250/93‑1 (MQ=255)

                                        tatTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCctct                                                    >  2:330758/1‑93 (MQ=255)

                                        tatTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCctct                                                    >  1:198669/1‑93 (MQ=255)

                                          ttAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTa                                                  >  1:5223810/1‑93 (MQ=255)

                                          ttAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTa                                                  >  2:2312684/1‑93 (MQ=255)

                                          ttAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTa                                                  <  1:3050509/93‑1 (MQ=255)

                                            aaCGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACt                                                >  2:1447587/1‑93 (MQ=255)

                                             aCGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTc                                               >  2:3064183/1‑93 (MQ=255)

                                             aCGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTc                                               >  2:516714/1‑93 (MQ=255)

                                             aCGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTc                                               <  1:4802460/93‑1 (MQ=255)

                                             aCGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTc                                               >  2:1082133/1‑93 (MQ=255)

                                             aCGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCCATAGTCCTCTTACTc                                               >  1:4884683/1‑93 (MQ=255)

                                              cGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCg                                              >  2:5577857/1‑93 (MQ=255)

                                                    gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATta                                        >  2:3858871/1‑93 (MQ=255)

                                                     cATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATtaa                                       >  1:3132220/1‑93 (MQ=255)

                                                      aTTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATtaat                                      <  2:2974323/93‑1 (MQ=255)

                                                       ttGTGGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATtaata                                     >  1:825108/1‑93 (MQ=255)

                                                         gtgGGGATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTaataat                                   <  1:1433405/93‑1 (MQ=255)

                                                           ggggATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCg                                 <  1:486541/93‑1 (MQ=255)

                                                             ggATTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTa                               <  2:4450189/93‑1 (MQ=255)

                                                              gATTGGACTCAGTGATGGGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTaa                              <  1:533663/93‑1 (MQ=255)

                                                               aTTGGACTCAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAac                             >  1:3529206/1‑93 (MQ=255)

                                                               aTTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTAGAGTCCTCTTACTCGTCATTAATAATCGTAac                             <  1:4489811/93‑1 (MQ=255)

                                                                   gACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAacaata                         <  2:2837244/93‑1 (MQ=255)

                                                                    aCTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAacaatac                        <  1:5409434/93‑1 (MQ=255)

                                                                       cAGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAacaatacaat                     <  2:2593218/93‑1 (MQ=255)

                                                                        aGTGATGTGATCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACaatt                    <  2:6109808/93‑1 (MQ=255)

                                                                         gtgatgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACaatta                   <  1:3885918/93‑1 (MQ=255)

                                                                           gatgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACaattaat                 >  2:4255648/1‑93 (MQ=255)

                                                                            atgtgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCTTCTTACTCGTCATTAATAATCGTAACAATACaattaatt                <  1:7508803/93‑1 (MQ=255)

                                                                            atgtgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACaattaatt                <  2:1991098/93‑1 (MQ=255)

                                                                            atgtgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACaattaatt                >  1:5889738/1‑93 (MQ=255)

                                                                            atgtgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACaattaatt                <  2:5998928/93‑1 (MQ=255)

                                                                              gtgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAAttaattaa              <  2:3172162/93‑1 (MQ=255)

                                                                              gtgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAAttaattaa              <  1:3736447/93‑1 (MQ=255)

                                                                               tgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAg             >  1:6443624/1‑93 (MQ=255)

                                                                               tgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAg             >  2:7137473/1‑93 (MQ=255)

                                                                               tgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAg             <  1:7206835/93‑1 (MQ=255)

                                                                               tgatCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAg             >  1:1190771/1‑93 (MQ=255)

                                                                                 atCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGaa           <  2:4720909/93‑1 (MQ=255)

                                                                                  tCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAAt          >  1:6788090/1‑93 (MQ=255)

                                                                                  tCATATGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAAt          <  1:1508246/93‑1 (MQ=255)

                                                                                  tCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAAt          >  1:6564936/1‑93 (MQ=255)

                                                                                  tCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAAt          >  2:5729313/1‑93 (MQ=255)

                                                                                    atatGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATaa        >  2:5368499/1‑93 (MQ=255)

                                                                                    atatGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATaa        <  1:5975623/93‑1 (MQ=255)

                                                                                    atatGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATaa        >  2:6895406/1‑93 (MQ=255)

                                                                                     tatGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATaaa       <  2:3743384/93‑1 (MQ=255)

                                                                                     tatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATaaa       <  1:6664358/93‑1 (MQ=255)

                                                                                      atGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAAc      <  2:733535/93‑1 (MQ=255)

                                                                                       tGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACa     >  2:435577/1‑93 (MQ=255)

                                                                                       tGGGCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACa     >  1:6182758/1‑93 (MQ=255)

                                                                                       tGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACa     <  1:5127224/93‑1 (MQ=255)

                                                                                        gggCGCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACaa    <  2:7256833/93‑1 (MQ=255)

TTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGT  >  NC_000913/3795664‑3795846

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