BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,579,794	T→C	Y101H (TAT→CAT)	yhhZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,579,794	0	T	C	99.2%	377.5 / ‑4.9	120	Y101H (TAT→CAT)	yhhZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (92/27); total (92/28)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 132 total aligned reads displayed.

CCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATT  >  NC_000913/3579704‑3579883

ccATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcat                                                                                         >  1:2584966/1‑93 (MQ=255)

 cATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcatt                                                                                        <  1:2968862/93‑1 (MQ=255)

     atttattAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcattatta                                                                                    >  1:4219968/1‑93 (MQ=255)

     atttattAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcattatta                                                                                    >  2:2227230/1‑93 (MQ=255)

      tttattAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcattattat                                                                                   >  1:299810/1‑93 (MQ=255)

      tttattAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcattattat                                                                                   <  2:813132/93‑1 (MQ=255)

      tttattAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGcattattat                                                                                   >  1:3581930/1‑93 (MQ=255)

       ttattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtata                                                                                  >  2:2253329/1‑93 (MQ=255)

       ttattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtata                                                                                  <  2:3212755/93‑1 (MQ=255)

       ttattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtata                                                                                  >  2:2069843/1‑93 (MQ=255)

       ttattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtata                                                                                  >  1:4391702/1‑93 (MQ=255)

        tattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtatat                                                                                 <  2:3912370/93‑1 (MQ=255)

        tattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtatat                                                                                 >  2:2746394/1‑93 (MQ=255)

         attaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATtatata                                                                                >  2:1860586/1‑93 (MQ=255)

          ttaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATAc                                                                               >  2:1632464/1‑93 (MQ=255)

             aTTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAAc                                                                            >  1:3973751/1‑93 (MQ=255)

              ttCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACt                                                                           >  1:4133410/1‑93 (MQ=255)

              ttCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACt                                                                           >  1:3919416/1‑93 (MQ=255)

                ccATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTaa                                                                         >  2:3874657/1‑93 (MQ=255)

                  aTTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAaga                                                                       >  2:1636453/1‑93 (MQ=255)

                   ttAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAagag                                                                      <  1:1651004/93‑1 (MQ=255)

                      aCAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGgcg                                                                   >  2:307087/1‑93 (MQ=255)

                      aCAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGgcg                                                                   >  2:1481486/1‑93 (MQ=255)

                       cAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAag                                                                        >  2:416618/1‑87 (MQ=255)

                       cAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAag                                                                        <  1:416618/87‑1 (MQ=255)

                       cAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGgcgc                                                                  >  2:3004366/1‑93 (MQ=255)

                         ataatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCtt                                                                >  2:2544849/1‑93 (MQ=255)

                         ataatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCtt                                                                >  1:2121428/1‑93 (MQ=255)

                         ataatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCtt                                                                >  1:3725338/1‑93 (MQ=255)

                          taatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCttt                                                               >  1:174134/1‑93 (MQ=255)

                            atGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCttttt                                                             <  2:279492/93‑1 (MQ=255)

                             tGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCtttttt                                                            >  2:1051087/1‑93 (MQ=255)

                             tGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCtttttt                                                            >  1:3956729/1‑93 (MQ=255)

                               aaCAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTAt                                                          >  2:2387907/1‑93 (MQ=255)

                                aCAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATc                                                         >  2:349304/1‑93 (MQ=255)

                                aCAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATc                                                         >  1:1770504/1‑93 (MQ=255)

                                 cAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCg                                                        >  2:597747/1‑93 (MQ=255)

                                 cAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCg                                                        >  2:910141/1‑93 (MQ=255)

                                  aaTTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCgg                                                       <  1:411065/93‑1 (MQ=255)

                                  aaTTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCgg                                                       >  2:2960048/1‑93 (MQ=255)

                                  aaTTATTTATGGGATTTGACTTCTATCGAATAAATAGATATGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCgg                                                       >  1:116473/1‑93 (MQ=255)

                                   attattTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGc                                                      >  1:3546492/1‑93 (MQ=255)

                                   attattTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGc                                                      >  1:1454512/1‑93 (MQ=255)

                                    ttatttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCt                                                     >  2:136718/1‑93 (MQ=255)

                                     tatttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTa                                                    >  1:1377277/1‑93 (MQ=255)

                                      atttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTAt                                                   >  2:4055738/1‑93 (MQ=255)

                                       tttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTAtt                                                  >  2:2891902/1‑93 (MQ=255)

                                       tttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTAtt                                                  >  1:2924030/1‑93 (MQ=255)

                                        ttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtc                                                 >  2:462808/1‑93 (MQ=255)

                                           tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtcatc                                              >  1:2501533/1‑93 (MQ=255)

                                           tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtcatc                                              >  2:4265283/1‑93 (MQ=255)

                                           tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtcatc                                              <  1:49570/93‑1 (MQ=255)

                                           tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATAAAACTAAAAGGCGCTTTTTTATCGGCTATtcatc                                              >  2:3731296/1‑93 (MQ=255)

                                               aTTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCa                                          >  1:4495476/1‑93 (MQ=255)

                                                tttGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAg                                         >  1:2581946/1‑93 (MQ=255)

                                                  tGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGAt                                       <  1:2641993/93‑1 (MQ=255)

                                                     cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAt                                    >  2:1302678/1‑93 (MQ=255)

                                                     cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAt                                    >  2:211492/1‑93 (MQ=255)

                                                      ttCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAtt                                   <  1:4517694/93‑1 (MQ=255)

                                                       tCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTg                                  >  1:565081/1‑93 (MQ=255)

                                                       tCTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTg                                  >  1:94436/1‑93 (MQ=255)

                                                        cTATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGa                                 >  1:3331807/1‑93 (MQ=255)

                                                         tATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGaa                                >  2:1921792/1‑93 (MQ=255)

                                                         tATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGaa                                >  1:3879802/1‑93 (MQ=255)

                                                         tATCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGaa                                >  1:3617141/1‑93 (MQ=255)

                                                           tCGAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGaaaa                              <  2:1554126/93‑1 (MQ=255)

                                                             gAATAAATAGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAAcc                            >  1:24974/1‑93 (MQ=255)

                                                              aataaataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTAGCGGCTATTCATCACCAGATCATTGAAAACCa                           <  2:4049862/93‑1 (MQ=255)

                                                              aataaataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCa                           <  2:4402387/93‑1 (MQ=255)

                                                                taaataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAAc                         >  1:4209626/1‑93 (MQ=255)

                                                                taaataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAAc                         >  2:4489051/1‑93 (MQ=255)

                                                                 aaataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACt                        >  1:2289921/1‑93 (MQ=255)

                                                                 aaataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACt                        <  1:626887/93‑1 (MQ=255)

                                                                  aataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTg                       >  2:777026/1‑93 (MQ=255)

                                                                  aataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTg                       >  2:4096548/1‑93 (MQ=255)

                                                                  aataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTg                       >  1:3048938/1‑93 (MQ=255)

                                                                  aataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTg                       >  2:1816112/1‑93 (MQ=255)

                                                                  aataGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTg                       >  1:4653546/1‑93 (MQ=255)

                                                                    taGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGa                     <  2:794035/93‑1 (MQ=255)

                                                                     aGATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGAt                    <  1:2474228/93‑1 (MQ=255)

                                                                      gATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATa                   >  2:3477079/1‑93 (MQ=255)

                                                                      gATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATa                   >  2:2174888/1‑93 (MQ=255)

                                                                      gATTTGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATa                   >  1:2251240/1‑93 (MQ=255)

                                                                        tttGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACa                 >  1:50510/1‑93 (MQ=255)

                                                                        tttGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACa                 >  1:1760702/1‑93 (MQ=255)

                                                                          tGGTAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGa               <  2:4273406/93‑1 (MQ=255)

                                                                             tAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAAc            >  1:1197193/1‑93 (MQ=255)

                                                                             tAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAAc            >  1:4068250/1‑93 (MQ=255)

                                                                             tAGATTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAAc            <  1:1212033/93‑1 (MQ=255)

                                                                                aTTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAAt         >  2:2435824/1‑93 (MQ=255)

                                                                                aTTGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAAt         >  2:587762/1‑93 (MQ=255)

                                                                                 ttGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATa        >  1:3833406/1‑93 (MQ=255)

                                                                                 ttGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATa        >  1:2994238/1‑93 (MQ=255)

                                                                                 ttGGAAAAGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATa        >  1:1611540/1‑93 (MQ=255)

                                                                                     aaaaGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTa    >  2:1442057/1‑93 (MQ=255)

                                                                                     aaaaGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTa    >  1:3984846/1‑93 (MQ=255)

                                                                                      aaaGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTAt   >  2:2794202/1‑93 (MQ=255)

                                                                                      aaaGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTAt   <  1:561133/93‑1 (MQ=255)

                                                                                       aaGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTAtt  <  2:3762530/93‑1 (MQ=255)

                                                                                       aaGCATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTAtt  >  1:807974/1‑93 (MQ=255)

CCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATT  >  NC_000913/3579704‑3579883

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