BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,729,644	T→C	M845T (ATG→ACG)	lhr →	Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,729,644	0	T	C	99.6%	775.8 / ‑4.9	239	M845T (ATG→ACG)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (137/101); total (137/102)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

Only 100 of 257 total aligned reads displayed.

CTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATG  >  NC_000913/1729554‑1729733

cTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGc                                                                                         >  2:7095943/1‑93 (MQ=255)

 tCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGAcgcg                                                                                        >  2:2945216/1‑93 (MQ=255)

  cAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGt                                                                                       <  2:6848577/93‑1 (MQ=255)

   aGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTc                                                                                      <  2:2823506/93‑1 (MQ=255)

   aGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTc                                                                                      >  2:6741898/1‑93 (MQ=255)

    gggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGCTCGTTATGATCTTCCCGCACTGGAACGTTAGACGAGTcc                                                                                     >  1:683673/1‑92 (MQ=255)

     ggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGc                                                                                    <  2:9276912/93‑1 (MQ=255)

     ggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGc                                                                                    >  2:5460173/1‑93 (MQ=255)

      gATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGCTTGACGCGTCGcc                                                                                   <  1:764337/93‑1 (MQ=255)

       atatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCCCTGGAACGTTTGACGCGTCGCCt                                                                                  >  2:6938742/1‑93 (MQ=255)

        tatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTg                                                                                 >  1:4871697/1‑93 (MQ=255)

        tatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTg                                                                                 >  1:3094045/1‑93 (MQ=255)

        tatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTg                                                                                 >  1:9466785/1‑93 (MQ=255)

         atCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGa                                                                                >  2:7301464/1‑93 (MQ=255)

         atCCGGATTTTCCGGTCAGTCTCGAAACCCTACGCGAATGTATGC‑AGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGaa                                                                               <  1:3815323/93‑1 (MQ=255)

          tCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGaa                                                                               >  1:1547786/1‑93 (MQ=255)

           ccGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAAc                                                                              >  2:8183518/1‑93 (MQ=255)

             ggATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACgg                                                                            <  1:4809636/93‑1 (MQ=255)

              gATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGt                                                                           <  1:5395967/93‑1 (MQ=255)

                ttttCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTgg                                                                         >  1:1468489/1‑93 (MQ=255)

                ttttCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTGCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTgg                                                                         <  2:9846188/93‑1 (MQ=255)

                ttttCCGGTCATTCTCGAAACCCTACGAGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTgg                                                                         <  1:1872580/93‑1 (MQ=255)

                    ccGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAAAGGAGGCgca                                                                     >  2:8219242/1‑91 (MQ=255)

                        tCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTc                                                                 >  1:9255891/1‑93 (MQ=255)

                        tCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTc                                                                 >  1:856811/1‑93 (MQ=255)

                            tctcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAat                                                             <  2:3240226/93‑1 (MQ=255)

                             ctcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAata                                                            <  1:3018342/93‑1 (MQ=255)

                             ctcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAata                                                            <  1:2232695/93‑1 (MQ=255)

                                 aaaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCg                                                        >  1:40135/1‑93 (MQ=255)

                                 aaaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCg                                                        >  2:84979/1‑93 (MQ=255)

                                   aCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGAt                                                      <  2:5537696/93‑1 (MQ=255)

                                    cccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATg                                                     <  2:7572838/93‑1 (MQ=255)

                                     ccTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGt                                                    <  1:6081616/93‑1 (MQ=255)

                                      cTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTa                                                   >  1:5814049/1‑93 (MQ=255)

                                       tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTaa                                                  <  1:8333719/93‑1 (MQ=255)

                                        aCGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAAc                                                 >  2:5148233/1‑93 (MQ=255)

                                        aCGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAAc                                                 <  2:2340462/93‑1 (MQ=255)

                                         cgcgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACg                                                <  2:7061735/93‑1 (MQ=255)

                                         cgcgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACg                                                >  1:6301326/1‑93 (MQ=255)

                                          gcgAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGa                                               >  2:7668318/1‑93 (MQ=255)

                                            gAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAATGAcc                                             <  2:359283/93‑1 (MQ=255)

                                             aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCa                                            >  2:7909792/1‑93 (MQ=255)

                                             aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCa                                            <  2:4111913/93‑1 (MQ=255)

                                             aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCa                                            <  1:4367739/93‑1 (MQ=255)

                                              aTGTGTGCAAGATGTTTATGATGTTCGCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGAccgc                                           <  1:3519642/93‑3 (MQ=255)

                                               tGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAACTTCAAATATTCGATGTAACGACCACt                                          >  2:9931390/1‑93 (MQ=255)

                                               tGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACt                                          >  1:1182683/1‑93 (MQ=255)

                                                gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         <  1:324368/93‑1 (MQ=255)

                                                gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         >  1:6388483/1‑93 (MQ=255)

                                                gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                         <  1:881119/93‑1 (MQ=255)

                                                 tCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTAc                                        <  1:8121790/93‑1 (MQ=255)

                                                  cTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACg                                       >  1:9191044/1‑93 (MQ=255)

                                                   tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATc                                                          >  2:3178428/1‑73 (MQ=255)

                                                   tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                      >  1:1070173/1‑93 (MQ=255)

                                                   tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                      >  1:2129762/1‑93 (MQ=255)

                                                   tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                      <  1:1407343/93‑1 (MQ=255)

                                                   tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                      >  2:2465755/1‑93 (MQ=255)

                                                   tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTAAAATATCCGATGTAACGACCACAACGc                                      >  1:5661479/1‑93 (MQ=255)

                                                    gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGcc                                     >  2:1743774/1‑93 (MQ=255)

                                                    gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGcc                                     <  1:4683986/93‑1 (MQ=255)

                                                      aaGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACCCCCt                                   >  2:2393233/1‑93 (MQ=255)

                                                       aGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                  >  1:1389757/1‑93 (MQ=255)

                                                       aGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                  >  2:5849664/1‑93 (MQ=255)

                                                       aGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                  >  2:6441938/1‑93 (MQ=255)

                                                       aGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                  <  2:4264882/93‑1 (MQ=255)

                                                       aGATGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                  >  1:7118029/1‑93 (MQ=255)

                                                         aTGTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGc                                <  1:2565757/93‑1 (MQ=255)

                                                           gTTTATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCt                              >  2:8915147/1‑93 (MQ=255)

                                                             ttATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCttt                            <  2:4609901/93‑1 (MQ=255)

                                                             ttATGATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCttt                            <  2:10038795/93‑1 (MQ=255)

                                                                 gATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGcc                        >  1:6492414/1‑93 (MQ=255)

                                                                 gATCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGcc                        >  2:1439993/1‑93 (MQ=255)

                                                                  aTCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCca                       <  1:3542632/93‑1 (MQ=255)

                                                                   tCTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCcac                      >  1:5646904/1‑93 (MQ=255)

                                                                    ccgcccGCACAGGAACGGATGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGTCCACTACGCCCTCGCCTTTCGCcaca                     <  1:7882724/90‑1 (MQ=255)

                                                                    cTTCCCGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGCTGTAAGGACCACTACTCCCTCGCCTTTCGCcaaa                     >  1:1046105/1‑91 (MQ=255)

                                                                       cccGAACTGGAACGTTTTACGCGTCGCCTGAACGGTGGCGAAATTCAAATATTCGATGTAAGGTCCACTACGCCCGCGACTTTCACCACAAGt                  >  2:7060846/1‑93 (MQ=255)

                                                                        ccGCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGtt                 <  2:2129499/93‑1 (MQ=255)

                                                                          gCACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTa               <  2:8012869/93‑1 (MQ=255)

                                                                          cgccTGGGACGGTTGACCGGTTGCCTGAACGGGGGGCAAAATCAAAAAACCGGAGTAACCGCCCCTTCCCCCTCGCCTTTCGCCACAAGTTTa               <  2:475318/90‑1 (MQ=255)

                                                                           cACTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAt              >  2:3818001/1‑93 (MQ=255)

                                                                            aCTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAtt             >  1:8709619/1‑93 (MQ=255)

                                                                            aCTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAtt             >  2:9759601/1‑93 (MQ=255)

                                                                             cTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTg            >  1:1294181/1‑93 (MQ=255)

                                                                             cTGGAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTg            <  1:6013142/93‑1 (MQ=255)

                                                                               ggAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGtt          >  2:5383152/1‑93 (MQ=255)

                                                                               ggAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGtt          >  2:1595357/1‑93 (MQ=255)

                                                                                gAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGTCACAAGTTTATTGTTc         >  1:8986621/1‑93 (MQ=255)

                                                                                gAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc         <  2:1827542/93‑1 (MQ=255)

                                                                                gAACGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc         <  2:9422496/93‑1 (MQ=255)

                                                                                 aaCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg        <  1:2721394/93‑1 (MQ=255)

                                                                                 aaCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg        >  1:708506/1‑93 (MQ=255)

                                                                                 aaCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg        <  1:7723801/93‑1 (MQ=255)

                                                                                  aCGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCgg       <  2:6518777/93‑1 (MQ=255)

                                                                                   cGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGc      <  2:8052222/93‑1 (MQ=255)

                                                                                   cGTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGc      <  1:9568487/93‑1 (MQ=255)

                                                                                    gTTTGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCt     >  2:8818479/1‑93 (MQ=255)

                                                                                     tttGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTa    >  1:459459/1‑93 (MQ=255)

                                                                                     tttGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTa    <  1:4831120/93‑1 (MQ=255)

                                                                                       tGACGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATg  <  1:7259195/93‑1 (MQ=255)

CTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATG  >  NC_000913/1729554‑1729733

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