BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,389,803	A→G	P29P (CCT→CCC)	b1327 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,389,803	0	A	G	99.0%	310.2 / ‑4.3	100	P29P (CCT→CCC)	b1327	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (48/51); total (48/52)
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.57e-01

Only 100 of 103 total aligned reads displayed.

GCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGA  >  NC_000913/1389714‑1389892

gCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTgggggg                                                                                        <  2:559732/93‑1 (MQ=255)

 cATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGa                                                                                       <  1:2286503/93‑1 (MQ=255)

  aTAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGaa                                                                                      <  1:2778436/93‑1 (MQ=255)

   tAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGaaa                                                                                     >  1:463519/1‑93 (MQ=255)

   tAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGaaa                                                                                     >  1:379833/1‑93 (MQ=255)

     ggTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATa                                                                                   <  2:3403419/93‑1 (MQ=255)

     ggTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATa                                                                                   <  2:1352477/93‑1 (MQ=255)

       ttCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 <  1:3713132/93‑1 (MQ=255)

       ttCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 <  1:1660075/93‑1 (MQ=255)

       ttCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 <  2:1263590/93‑1 (MQ=255)

       ttCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 <  2:3925245/93‑1 (MQ=255)

       ttCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 <  1:2511476/93‑1 (MQ=255)

       gtCCGGCCGTTTGTTCTGTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGa                                                                                 <  2:4311295/92‑1 (MQ=255)

         ccGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATc                                                                               <  1:4045475/93‑1 (MQ=255)

         ccGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATc                                                                               <  2:1870473/93‑1 (MQ=255)

          cGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCa                                                                              >  1:2869928/1‑93 (MQ=255)

          cGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCa                                                                              >  1:1722486/1‑93 (MQ=255)

            gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACt                                                                            >  1:879603/1‑93 (MQ=255)

              cGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGc                                                                          >  2:1097644/1‑93 (MQ=255)

              cGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGc                                                                          <  1:1378674/93‑1 (MQ=255)

               gTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCa                                                                         <  2:117135/93‑1 (MQ=255)

               gTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCa                                                                         >  1:2215745/1‑93 (MQ=255)

                tttGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAg                                                                        >  1:4285271/1‑93 (MQ=255)

                   gTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACa                                                                     >  2:1191093/1‑93 (MQ=255)

                    ttCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACaa                                                                    >  1:1501078/1‑93 (MQ=255)

                     tCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACaaa                                                                   <  1:3552275/93‑1 (MQ=255)

                     tCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACaaa                                                                   <  2:3022999/93‑1 (MQ=255)

                      cTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAag                                                                  >  2:4123477/1‑93 (MQ=255)

                      cTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAag                                                                  >  1:4076458/1‑93 (MQ=255)

                       tttAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAaga                                                                 <  1:2854769/93‑1 (MQ=255)

                       tttAACGCCCCCACCCGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAaga                                                                 >  2:1220590/1‑93 (MQ=255)

                         tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               <  2:1122449/93‑1 (MQ=255)

                         tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               <  1:4048828/93‑1 (MQ=255)

                         tAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAgaga                                                               >  1:1234541/1‑93 (MQ=255)

                           aCGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATg                                                             <  1:2315184/93‑1 (MQ=255)

                           aCGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATg                                                             <  1:4395110/93‑1 (MQ=255)

                             gCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGt                                                           <  2:1765700/93‑1 (MQ=255)

                             gCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGt                                                           >  1:1279086/1‑93 (MQ=255)

                             gCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGt                                                           <  1:4134989/93‑1 (MQ=255)

                              cccccGCCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGtt                                                          <  1:3866564/93‑1 (MQ=255)

                              cccccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGtt                                                          <  1:191616/93‑1 (MQ=255)

                                cccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTgg                                                        >  1:2597417/1‑93 (MQ=255)

                                cccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGGTgg                                                        >  2:2033003/1‑93 (MQ=255)

                                 ccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGa                                                       >  2:3486909/1‑93 (MQ=255)

                                 ccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGa                                                       >  1:4061756/1‑93 (MQ=255)

                                 ccACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGa                                                       >  1:37652/1‑93 (MQ=255)

                                    ccTGGGTGAGATAACACAATTGCCTGATACTGACGGGGTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGatta                                                    >  1:2718025/1‑93 (MQ=255)

                                     cTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattat                                                   <  2:1923368/93‑1 (MQ=255)

                                     cTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattat                                                   >  2:4469316/1‑93 (MQ=255)

                                      tGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattatt                                                  <  1:1887806/93‑1 (MQ=255)

                                       gggTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTg                                                 <  2:1038830/93‑1 (MQ=255)

                                        ggTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGc                                                <  2:1101175/93‑1 (MQ=255)

                                         gTGAGATAACACAATTGCCTGATACTGATGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGCTGGATTATTGCt                                               >  2:3809477/1‑93 (MQ=255)

                                         gTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCt                                               >  2:898472/1‑93 (MQ=255)

                                          tGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTg                                              <  2:3005757/93‑1 (MQ=255)

                                            agaTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGtt                                            >  1:3221737/1‑93 (MQ=255)

                                            agaTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGtt                                            <  2:4036388/93‑1 (MQ=255)

                                             gaTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTa                                           <  1:2801455/93‑1 (MQ=255)

                                              aTAACACGATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAg                                          <  2:1776224/93‑1 (MQ=255)

                                              aTAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAg                                          <  2:349934/93‑1 (MQ=255)

                                               tAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGt                                         >  2:1486473/1‑93 (MQ=255)

                                               tAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAA‑AAATTACTGCCGACCAAGAGAAGGGTTGGTTATTTCTGGTtgtg                                        >  2:1740767/1‑93 (MQ=255)

                                                 acacAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGa                                       <  2:1913428/93‑1 (MQ=255)

                                                  cacaATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGaa                                      <  2:4049457/93‑1 (MQ=255)

                                                   acaATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAg                                     >  2:426974/1‑93 (MQ=255)

                                                   acaATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAg                                     >  1:4476843/1‑93 (MQ=255)

                                                    caATTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGc                                    <  2:3454459/93‑1 (MQ=255)

                                                      aTTGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTg                                  >  2:2502050/1‑93 (MQ=255)

                                                       ttGCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTgc                                 >  2:3836031/1‑92 (MQ=255)

                                                        tGCCTTATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGTTCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGAc                                <  2:411550/93‑1 (MQ=255)

                                                         gCCTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACt                               >  2:2382615/1‑93 (MQ=255)

                                                           cTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACttt                             >  2:462884/1‑93 (MQ=255)

                                                           cTGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACttt                             <  1:2129692/93‑1 (MQ=255)

                                                            tGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTtta                            >  1:2476791/1‑93 (MQ=255)

                                                            tGATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTtta                            >  2:1163226/1‑93 (MQ=255)

                                                             taTACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttat                           <  2:3399461/92‑1 (MQ=255)

                                                             gATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttat                           <  1:4191261/93‑1 (MQ=255)

                                                             gATACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttat                           >  1:4373246/1‑93 (MQ=255)

                                                               tACTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttatta                         <  2:3122700/93‑1 (MQ=255)

                                                                 cTGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTattatt                       >  2:3102987/1‑93 (MQ=255)

                                                                  tGACGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtc                      <  2:937398/93‑1 (MQ=255)

                                                                     cGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                   >  2:945178/1‑93 (MQ=255)

                                                                     cGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                   >  1:883666/1‑93 (MQ=255)

                                                                     cGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                   <  1:2546243/93‑1 (MQ=255)

                                                                     cGGGTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc                   >  1:3998549/1‑93 (MQ=255)

                                                                      gggTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatca                  <  2:590631/93‑1 (MQ=255)

                                                                      gggTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatca                  <  1:253212/93‑1 (MQ=255)

                                                                       ggTTTCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTcatcat                 <  2:1351454/93‑1 (MQ=255)

                                                                       ggTTTCATCAAATTTTGGGGGAAAATAGATCACTGCAGAAAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTcatcat                 <  2:3130364/93‑1 (MQ=255)

                                                                         tttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc               <  2:4394601/93‑1 (MQ=255)

                                                                         tttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc               >  2:2597782/1‑93 (MQ=255)

                                                                         tttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc               >  2:69991/1‑93 (MQ=255)

                                                                         tttCATCAAATTTTGGGGAGAAATAGGTCGCTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGAGAAGCGGACTTTATTATTCACCATTc               <  2:2664660/93‑1 (MQ=255)

                                                                          ttCATCAAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTcc              >  2:1710871/1‑93 (MQ=255)

                                                                           tcatcaAATTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCg             >  1:2100576/1‑93 (MQ=255)

                                                                                 aaTTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTc       >  1:2012900/1‑93 (MQ=255)

                                                                                 aaTTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTc       <  2:3190258/93‑1 (MQ=255)

                                                                                  aTTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCa      <  2:513513/93‑1 (MQ=255)

                                                                                  aTTTTGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCa      <  2:3919030/93‑1 (MQ=255)

                                                                                      tGGGGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGa  >  1:4564813/1‑93 (MQ=255)

GCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTCCGTTCCTCATGGA  >  NC_000913/1389714‑1389892

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