BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,676,843	T→C	R131R (CGT→CGC)	b1604 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,676,843	0	T	C	99.2%	378.9 / ‑4.3	119	R131R (CGT→CGC)	b1604	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (56/61); total (58/61)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 125 total aligned reads displayed.

CGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCA  >  NC_000913/1676754‑1676931

cgcgTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagc                                                                                       >  1:1099652/1‑93 (MQ=255)

 gcgTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagcc                                                                                      >  1:4174454/1‑93 (MQ=255)

 gcgTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTAAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagcc                                                                                      >  2:3926383/1‑93 (MQ=255)

   gTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccag                                                                                    >  1:302451/1‑93 (MQ=255)

   gTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccag                                                                                    >  1:593214/1‑93 (MQ=255)

     aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGTGCTTCTATCGCagccagcc                                                                                  >  2:2657788/1‑93 (MQ=255)

     aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccagcc                                                                                  <  2:1021784/93‑1 (MQ=255)

     aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACAGTCCAGGGATTCTATCGCagccagcc                                                                                  >  1:3304116/1‑93 (MQ=255)

      aTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAgccagcca                                                                                 <  1:584392/93‑1 (MQ=255)

      aTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAgccagcca                                                                                 >  1:4078494/1‑93 (MQ=255)

      aTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAgccagcca                                                                                 <  2:2408844/93‑1 (MQ=255)

       ttAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGccagccag                                                                                >  1:2053020/1‑93 (MQ=255)

         aaCGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGaa                                                                              <  1:3551099/93‑1 (MQ=255)

           cGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGt                                                                            >  1:1675223/1‑93 (MQ=255)

           cGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGt                                                                            >  2:449855/1‑93 (MQ=255)

               gttgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTGTATCGCAGCCAGCCAGAAGTAAat                                                                        >  2:2535683/1‑93 (MQ=255)

               gttgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAat                                                                        <  1:1042699/93‑1 (MQ=255)

               gttgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAat                                                                        >  2:2535682/1‑93 (MQ=255)

                 tgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatga                                                                      <  1:1661225/93‑1 (MQ=255)

                   tCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatgatg                                                                    <  1:4334459/93‑1 (MQ=255)

                   tCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatgatg                                                                    >  1:2478950/1‑93 (MQ=255)

                    cGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGc                                                                   <  2:361622/93‑1 (MQ=255)

                     gAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGcc                                                                  <  2:4241390/93‑1 (MQ=255)

                     gAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGcc                                                                  >  2:2252943/1‑93 (MQ=255)

                      aaCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCa                                                                 >  1:4628543/1‑93 (MQ=255)

                      aaCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCa                                                                 >  1:3676718/1‑93 (MQ=255)

                       aCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCCGAGCAAGCCAGAAGTCAATGATGCCAt                                                                >  1:2458775/1‑93 (MQ=255)

                        cTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATc                                                               >  1:261187/1‑93 (MQ=255)

                            cGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCa                                                           <  1:3848456/93‑1 (MQ=255)

                             gAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCaa                                                          >  1:2059044/1‑93 (MQ=255)

                               aaGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAg                                                        >  2:134908/1‑93 (MQ=255)

                                 gATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCg                                                      >  1:2776404/1‑93 (MQ=255)

                                 gATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCg                                                      >  2:1612012/1‑93 (MQ=255)

                                 gATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCg                                                      >  2:3323464/1‑93 (MQ=255)

                                 aaTCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCg                                                      <  1:1947484/92‑1 (MQ=255)

                                  aTCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCgg                                                     <  1:2145479/93‑1 (MQ=255)

                                   tCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGc                                                    >  2:731080/1‑93 (MQ=255)

                                   tCAGGCTGTTCTGATTGAACCGTCTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGc                                                    <  2:2584430/93‑1 (MQ=255)

                                     aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                  <  1:2272005/93‑1 (MQ=255)

                                     aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                  <  1:2272007/93‑1 (MQ=255)

                                     aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                  <  2:3076673/93‑1 (MQ=255)

                                     aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                  <  2:3311714/93‑1 (MQ=255)

                                       gCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaaa                                                <  1:4370702/93‑1 (MQ=255)

                                       gCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaaa                                                <  2:1068467/93‑1 (MQ=255)

                                       gCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaaa                                                <  2:4261409/93‑1 (MQ=255)

                                         tGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGc                                              <  2:179999/93‑1 (MQ=255)

                                          gTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCg                                             >  2:555225/1‑93 (MQ=255)

                                          gTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCg                                             <  1:4502622/93‑1 (MQ=255)

                                           ttCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGa                                            <  1:1618659/93‑1 (MQ=255)

                                            tCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAGGCGGCAAAAGCGaa                                           <  2:1849220/93‑1 (MQ=255)

                                              tGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAg                                         <  2:1369092/93‑1 (MQ=255)

                                              tGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAg                                         <  2:1369093/93‑1 (MQ=255)

                                               gATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAgg                                        <  2:219247/93‑1 (MQ=255)

                                               gATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAgg                                        <  1:3748879/93‑1 (MQ=255)

                                                aTTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGt                                       <  2:4266636/93‑1 (MQ=255)

                                                   taaCCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGcc                                    <  1:2517367/92‑1 (MQ=255)

                                                   gAACCGTTTGATGCGGGCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGcc                                    <  2:1959693/93‑1 (MQ=255)

                                                      ccGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTAc                                 <  2:3321503/93‑1 (MQ=255)

                                                       cGTTTGACACGCTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAATGTGCCTAct                                <  2:4076773/93‑1 (MQ=255)

                                                         tttGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActct                              >  1:2346515/1‑93 (MQ=255)

                                                         tttGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActct                              >  1:234684/1‑93 (MQ=255)

                                                          ttGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCtt                             <  1:2183854/93‑1 (MQ=255)

                                                           tGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCttt                            >  2:4039076/1‑93 (MQ=255)

                                                           tGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCttt                            <  1:1659906/93‑1 (MQ=255)

                                                            gACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTc                           <  2:2680973/93‑1 (MQ=255)

                                                               acGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTAc                        >  2:4117132/1‑93 (MQ=255)

                                                                cGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACa                       >  2:4249183/1‑93 (MQ=255)

                                                                cGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACa                       >  2:1029660/1‑93 (MQ=255)

                                                                 ggTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAt                      <  2:4544583/93‑1 (MQ=255)

                                                                  gTCACCGTCCAGGGGTTCTATCGCAGCCAGCCAGAAGTTAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTGCTCTTTCTACAtg                     >  2:1384425/1‑92 (MQ=255)

                                                                  gTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtc                     >  2:604974/1‑93 (MQ=255)

                                                                    cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                   <  2:2240262/93‑1 (MQ=255)

                                                                    cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                   >  2:180640/1‑93 (MQ=255)

                                                                    cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                   >  2:3203284/1‑93 (MQ=255)

                                                                    cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                   >  2:3203285/1‑93 (MQ=255)

                                                                    cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                   <  1:810353/93‑1 (MQ=255)

                                                                     aCCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgtt                  >  2:2241828/1‑93 (MQ=255)

                                                                     aCCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgtt                  >  1:3500717/1‑93 (MQ=255)

                                                                     aCCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgtt                  >  2:3415337/1‑93 (MQ=255)

                                                                      ccGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgttc                 <  1:2767042/93‑1 (MQ=255)

                                                                      ccGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgttc                 <  1:1876527/93‑1 (MQ=255)

                                                                      ccGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgttc                 <  1:4116289/93‑1 (MQ=255)

                                                                            aGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAATGTGCCTACTCTTTCTACATCGTTCGTCaaa           <  1:1862446/93‑1 (MQ=255)

                                                                            aGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCaaa           <  1:1462989/93‑1 (MQ=255)

                                                                            aGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCaaa           <  1:1751474/93‑1 (MQ=255)

                                                                             gggCTTCTATCGACGCCAGCCAGAAGTCAATAATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGACTACTCTTTCTAGATCGTTCGTCAAAt          >  2:1639851/1‑93 (MQ=255)

                                                                              ggCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATc         >  2:10030/1‑93 (MQ=255)

                                                                                cTTCTATCGCAGCCAGCCAGAAGTGAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa       <  2:2587483/93‑1 (MQ=255)

                                                                                cTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa       <  2:681704/93‑1 (MQ=255)

                                                                                cTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa       <  2:754903/93‑1 (MQ=255)

                                                                                cTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa       <  2:1865933/93‑1 (MQ=255)

                                                                                 ttCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGAt      >  1:4320925/1‑93 (MQ=255)

                                                                                  tCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATg     <  1:3464852/93‑1 (MQ=255)

                                                                                  tCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATg     <  2:2366501/93‑1 (MQ=255)

                                                                                   cTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGc    >  2:4496179/1‑93 (MQ=255)

                                                                                   cTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGc    >  1:1410571/1‑93 (MQ=255)

                                                                                    tATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGcc   <  1:3060413/93‑1 (MQ=255)

                                                                                    tATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGcc   <  1:4073501/93‑1 (MQ=255)

                                                                                    caTCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGcc   <  1:3808863/92‑1 (MQ=255)

                                                                                     aTCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCa  <  2:1836010/93‑1 (MQ=255)

CGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCA  >  NC_000913/1676754‑1676931

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