BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,481,375	C→T	I97I (ATC→ATT)	hrpA →	helicase, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,481,375	0	C	T	99.2%	443.6 / ‑5.1	126	I97I (ATC→ATT)	hrpA	helicase, ATP‑dependent
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (63/62); total (63/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 130 total aligned reads displayed.

CGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCA  >  NC_000913/1481286‑1481464

cGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAttgtc                                                                                        >  1:2163247/1‑93 (MQ=255)

  aCCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTTGc                                                                                      >  1:4388852/1‑93 (MQ=255)

  aCCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGc                                                                                      >  2:4001123/1‑93 (MQ=255)

   ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGcc                                                                                     <  2:4061483/93‑1 (MQ=255)

   ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGcc                                                                                     >  2:321959/1‑93 (MQ=255)

    cGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCg                                                                                    >  1:40411/1‑93 (MQ=255)

     ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCgg                                                                                   >  1:2276801/1‑93 (MQ=255)

     ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTAGTGATCACCAGGTGGTGATTGTCGCCgg                                                                                   >  1:3747854/1‑93 (MQ=255)

      gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCggg                                                                                  <  2:1568536/93‑1 (MQ=255)

      gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCggg                                                                                  <  2:4302788/93‑1 (MQ=255)

      gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCggg                                                                                  <  1:1772585/93‑1 (MQ=255)

      gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCggg                                                                                  <  1:569659/93‑1 (MQ=255)

      gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCggg                                                                                  <  1:4635991/93‑1 (MQ=255)

        aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGa                                                                                >  2:2971356/1‑93 (MQ=255)

        aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGa                                                                                <  1:2408200/93‑1 (MQ=255)

           tACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAAc                                                                             >  1:345424/1‑93 (MQ=255)

           tACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAAc                                                                             <  1:1325423/93‑1 (MQ=255)

           tACTTATCCTGACAATTTACCGCTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAAc                                                                             <  1:932486/93‑1 (MQ=255)

             cTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACgg                                                                           >  2:3803219/1‑93 (MQ=255)

             cTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACgg                                                                           <  1:2792583/93‑1 (MQ=255)

             cTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACgg                                                                           <  1:2253031/93‑1 (MQ=255)

             cTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACgg                                                                           <  1:1559755/93‑1 (MQ=255)

              ttATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACggg                                                                          <  1:1593806/93‑1 (MQ=255)

                aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAATCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGtt                                                                        >  2:1233019/1‑93 (MQ=255)

                 tCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGCGGTGATTGTCGCCGGGGACACGGGTTc                                                                       >  2:3025589/1‑93 (MQ=255)

                  ccTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCt                                                                      <  1:2704282/93‑1 (MQ=255)

                   cTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTg                                                                     >  2:2356129/1‑93 (MQ=255)

                    tGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTgg                                                                    >  1:4301074/1‑93 (MQ=255)

                    tGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTgg                                                                    >  2:390450/1‑93 (MQ=255)

                     gACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGt                                                                   >  1:3078629/1‑93 (MQ=255)

                     gACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGt                                                                   >  1:977770/1‑93 (MQ=255)

                       cAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTaa                                                                 >  2:103650/1‑93 (MQ=255)

                        caTTTACCGGTTAGTCAGAAAAAACAGGACATCCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTaaa                                                                <  1:1878973/92‑1 (MQ=255)

                        aaTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTaaa                                                                >  1:1496095/1‑93 (MQ=255)

                         aTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTaaaa                                                               <  2:2622162/93‑1 (MQ=255)

                         aTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTaaaa                                                               <  2:816236/93‑1 (MQ=255)

                               cGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTc                                                         <  1:2942693/93‑1 (MQ=255)

                                ggTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCa                                                        <  2:3345022/93‑1 (MQ=255)

                                ggTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCa                                                        <  2:2947844/93‑1 (MQ=255)

                                 gTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAg                                                       >  2:3177767/1‑93 (MQ=255)

                                   tAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGtt                                                     <  1:3114397/93‑1 (MQ=255)

                                   tAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGtt                                                     <  1:1423043/93‑1 (MQ=255)

                                   tAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGtt                                                     <  1:2303033/93‑1 (MQ=255)

                                     gTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAc                                                   >  2:2557233/1‑93 (MQ=255)

                                       cAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCg                                                 >  2:2177878/1‑93 (MQ=255)

                                        aGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGa                                                >  1:979215/1‑93 (MQ=255)

                                        aGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGa                                                >  1:4533779/1‑93 (MQ=255)

                                          aaaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCgggg                                                                                 <  2:4300098/58‑1 (MQ=255)

                                           aaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                             <  1:1666230/93‑1 (MQ=255)

                                            aaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAACCGACTCAGTTACCGAAAAt                                            >  2:3320273/1‑93 (MQ=255)

                                            aaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAACAt                                            >  1:4192840/1‑93 (MQ=255)

                                            aaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAAt                                            >  1:2203321/1‑93 (MQ=255)

                                             aaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATc                                           <  2:2031591/93‑1 (MQ=255)

                                             aaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATc                                           <  1:1136288/93‑1 (MQ=255)

                                             aaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATc                                           >  1:4206981/1‑93 (MQ=255)

                                              aaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCt                                          >  2:855271/1‑93 (MQ=255)

                                                cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGt                                        >  1:3918403/1‑93 (MQ=255)

                                                cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTGCTGGTAAAACGACTCAGTTACCGAAAGTCTGt                                        >  2:747395/1‑93 (MQ=255)

                                                 aGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                       <  1:2739559/93‑1 (MQ=255)

                                                 aGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                       <  1:268059/93‑1 (MQ=255)

                                                 aGGACATGCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                       <  2:459643/93‑1 (MQ=255)

                                                 aGGACAGTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                       <  2:459485/93‑1 (MQ=255)

                                                  ggACATTCTCGAAGCGATTCGTGATCAGCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                      <  1:1854967/93‑1 (MQ=255)

                                                  ggACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                      <  1:4429043/93‑1 (MQ=255)

                                                  ggACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                      >  1:3571920/1‑93 (MQ=255)

                                                  cgACATTCTCGAAGCGATTCGTGATCACCAGGTGGTCATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                      <  1:1855116/92‑1 (MQ=255)

                                                   gACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAACATCTGTATg                                     >  1:3069743/1‑93 (MQ=255)

                                                   gACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATg                                     <  1:4603413/93‑1 (MQ=255)

                                                    aCATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATgg                                    <  1:1641169/93‑1 (MQ=255)

                                                    aCATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATgg                                    <  2:4014298/93‑1 (MQ=255)

                                                    aCATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATgg                                    <  1:605833/93‑1 (MQ=255)

                                                     cATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGa                                   >  2:375891/1‑93 (MQ=255)

                                                     cATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGa                                   <  1:1227513/93‑1 (MQ=255)

                                                     cATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGGTACCGAAAATCTGTATgcg                                   >  2:2440731/1‑91 (MQ=255)

                                                        tctcGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCt                                <  1:653972/93‑1 (MQ=255)

                                                        tcctgAAGCGATTCGTGAACACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCt                                <  2:799478/89‑1 (MQ=255)

                                                         ctcGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTg                               >  1:2825347/1‑93 (MQ=255)

                                                          tcGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                              >  2:2839621/1‑93 (MQ=255)

                                                          tcGAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                              <  1:2184267/93‑1 (MQ=255)

                                                            gAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                            >  2:1222239/1‑93 (MQ=255)

                                                            gAAGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                            >  1:4043275/1‑93 (MQ=255)

                                                             aaGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgc                           >  1:778665/1‑93 (MQ=255)

                                                              aGCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                          >  1:4662072/1‑93 (MQ=255)

                                                               gCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                         >  1:4332875/1‑93 (MQ=255)

                                                               gCGATTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGGTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                         >  1:26935/1‑93 (MQ=255)

                                                                  aTTCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCggg                      >  1:4702048/1‑93 (MQ=255)

                                                                    tCGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAt                    <  1:717485/93‑1 (MQ=255)

                                                                     cGTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAtt                   <  2:1772363/93‑1 (MQ=255)

                                                                      gTGATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTa                  >  2:1704441/1‑93 (MQ=255)

                                                                        gATCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                <  2:4346483/93‑1 (MQ=255)

                                                                         aTCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg               <  2:23554/93‑1 (MQ=255)

                                                                          tCACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAgg              <  1:1352518/93‑1 (MQ=255)

                                                                           cACCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGa             >  1:3552563/1‑93 (MQ=255)

                                                                            aCCAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc            <  2:4533622/93‑1 (MQ=255)

                                                                              cAGGTGGTGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg          <  1:432775/93‑1 (MQ=255)

                                                                                ggtggtGATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt        <  1:3061153/93‑1 (MQ=255)

                                                                                  tggtgATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg      >  1:274618/1‑93 (MQ=255)

                                                                                  tggtgATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg      >  1:4130003/1‑93 (MQ=255)

                                                                                    gtgATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGAATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGc    >  1:764493/1‑93 (MQ=255)

                                                                                     tgATTGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGcc   >  1:4220069/1‑93 (MQ=255)

CGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCA  >  NC_000913/1481286‑1481464

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: `‑`