BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	853,690	A→G	N274S (AAC→AGC)	ybiR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	853,690	0	A	G	86.1%	656.9 / 88.4	274	N274S (AAC→AGC)	ybiR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (26/12); new base G (111/125); total (137/137)
Fisher's exact test for biased strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 274 total aligned reads displayed.

TGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATATTGTTGC  >  NC_000913/853600‑853779

tggactggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACg                                                                                         >  2:193464/1‑93 (MQ=255)

 ggactggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGt                                                                                        <  1:1813046/93‑1 (MQ=255)

  gactggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTg                                                                                       >  1:7008961/1‑93 (MQ=255)

   actggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGa                                                                                      <  1:4032441/93‑1 (MQ=255)

   actggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGa                                                                                      <  1:6197621/93‑1 (MQ=255)

     tggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGt                                                                                    >  2:2763219/1‑93 (MQ=255)

      ggacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTc                                                                                   <  1:1023169/93‑1 (MQ=255)

       gacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCa                                                                                  <  2:5461023/93‑1 (MQ=255)

       gacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCa                                                                                  <  2:5406084/93‑1 (MQ=255)

       gacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCa                                                                                  >  1:3782634/1‑93 (MQ=255)

       gacGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCa                                                                                  >  1:972498/1‑93 (MQ=255)

        acGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCat                                                                                 <  1:1542750/93‑1 (MQ=255)

         cGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCatc                                                                                <  1:4737023/93‑1 (MQ=255)

           ctgctgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTTGCGTGAGTCatcta                                                                              <  2:5190111/93‑1 (MQ=255)

             gctgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCatctatc                                                                            <  1:5579422/93‑1 (MQ=255)

             gctgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCatctatc                                                                            >  1:7278901/1‑93 (MQ=255)

             gctgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCatctatc                                                                            <  1:7477940/93‑1 (MQ=255)

             gctgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCatctatc                                                                            <  2:2538485/93‑1 (MQ=255)

              ctgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGGAGCGTGGGGCatctatct                                                                           >  1:3628547/1‑93 (MQ=255)

               tgctgGTGTTTTTGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTg                                                                          >  2:1536540/1‑93 (MQ=255)

               tgctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTg                                                                          >  2:2348228/1‑93 (MQ=255)

                 ctgGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGaa                                                                        <  2:6777344/93‑1 (MQ=255)

                   ggTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAAcc                                                                      <  2:6363268/93‑1 (MQ=255)

                   ggTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAAcc                                                                      <  1:195480/93‑1 (MQ=255)

                     tgtTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCg                                                                    <  2:4603258/93‑1 (MQ=255)

                     tgtTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCg                                                                    >  1:4452215/1‑93 (MQ=255)

                     tgtTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCg                                                                    <  2:5036022/93‑1 (MQ=255)

                       tttATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCggg                                                                  <  2:1859511/93‑1 (MQ=255)

                        ttATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGt                                                                 <  2:7287051/93‑1 (MQ=255)

                         tATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGtt                                                                >  1:4135695/1‑93 (MQ=255)

                           tGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTAt                                                              >  2:613280/1‑93 (MQ=255)

                           tGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTAt                                                              <  1:34015/93‑1 (MQ=255)

                            ggCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATg                                                             >  1:7331933/1‑93 (MQ=255)

                            ggCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATg                                                             <  2:471266/93‑1 (MQ=255)

                             gCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATgg                                                            <  1:1772711/93‑1 (MQ=255)

                              cGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGt                                                           <  1:2879989/93‑1 (MQ=255)

                              cGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGt                                                           >  2:3903944/1‑93 (MQ=255)

                               gATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGtt                                                          <  2:4820614/93‑1 (MQ=255)

                               gATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGtt                                                          >  2:1380614/1‑93 (MQ=255)

                               gATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGtt                                                          >  2:2051630/1‑93 (MQ=255)

                                 tGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTaa                                                        <  1:2504701/93‑1 (MQ=255)

                                 tGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTaa                                                        <  2:5309469/93‑1 (MQ=255)

                                  gTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAAc                                                       <  1:5753431/93‑1 (MQ=255)

                                        cGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAAt                                                 >  2:592948/1‑93 (MQ=255)

                                        cGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAAt                                                 <  2:2199197/93‑1 (MQ=255)

                                         gACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATc                                                <  2:697345/93‑1 (MQ=255)

                                         gACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATc                                                <  2:3365713/93‑1 (MQ=255)

                                          aCGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCg                                               >  1:6486585/1‑93 (MQ=255)

                                           cGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCgg                                              <  1:5790583/93‑1 (MQ=255)

                                           cGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCgg                                              >  2:2904031/1‑93 (MQ=255)

                                            gTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGt                                             >  1:50492/1‑93 (MQ=255)

                                             tCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGtt                                            <  1:1129329/93‑1 (MQ=255)

                                             tCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGtt                                            <  2:4760312/93‑1 (MQ=255)

                                             tCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGtt                                            >  1:3098070/1‑93 (MQ=255)

                                                aTTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTAt                                         >  2:6298432/1‑93 (MQ=255)

                                                aTTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGCTATGGTTAACGGCAATCGGCTTAt                                         >  1:1704271/1‑93 (MQ=255)

                                                 tttACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATc                                        >  2:4794969/1‑93 (MQ=255)

                                                 tttACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATc                                        <  1:3740305/93‑1 (MQ=255)

                                                  ttACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCg                                       <  1:931893/93‑1 (MQ=255)

                                                      tGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAgg                                   <  2:3622823/93‑1 (MQ=255)

                                                       gACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGt                                  >  2:4248095/1‑93 (MQ=255)

                                                       gACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGt                                  <  1:428070/93‑1 (MQ=255)

                                                        aCCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAgggg                                 >  1:6236790/1‑91 (MQ=255)

                                                        aCCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTg                                 <  1:2079655/93‑1 (MQ=255)

                                                          ccAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGAt                               >  2:3781313/1‑93 (MQ=255)

                                                          ccAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGAt                               >  2:4427285/1‑93 (MQ=255)

                                                           cAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATc                              <  2:238017/93‑1 (MQ=255)

                                                           cAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATc                              >  1:6333835/1‑93 (MQ=255)

                                                            aGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCa                             <  1:2155153/93‑1 (MQ=255)

                                                            aGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCa                             >  1:6040383/1‑93 (MQ=255)

                                                            aGCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCa                             <  1:7134360/93‑1 (MQ=255)

                                                            aGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCa                             >  1:73990/1‑93 (MQ=255)

                                                             gCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTTAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATcnn                            >  1:520289/1‑91 (MQ=255)

                                                             gCTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAg                            <  2:5107904/93‑1 (MQ=255)

                                                              cTTCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGt                           <  2:6480572/93‑1 (MQ=255)

                                                               ttCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTa                          >  1:3196417/1‑93 (MQ=255)

                                                               ttCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTa                          <  1:2846794/93‑1 (MQ=255)

                                                                tCCAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTaa                         <  1:2443132/93‑1 (MQ=255)

                                                                 ccAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAAc                        <  2:6273517/93‑1 (MQ=255)

                                                                  cAGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACg                       <  1:5278502/93‑1 (MQ=255)

                                                                   aGCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGt                      <  2:770184/93‑1 (MQ=255)

                                                                    gCGTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTg                     >  2:1056422/1‑93 (MQ=255)

                                                                      gTTGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGcc                   <  1:3378988/93‑1 (MQ=255)

                                                                      gTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGcc                   <  1:276752/93‑1 (MQ=255)

                                                                       ttGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCg                  <  2:4460586/93‑1 (MQ=255)

                                                                       ttGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCg                  >  2:5744434/1‑93 (MQ=255)

                                                                        tGCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGa                 <  2:2954401/93‑1 (MQ=255)

                                                                        tGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGa                 <  2:4137834/93‑1 (MQ=255)

                                                                         gCAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAg                >  1:6669391/1‑93 (MQ=255)

                                                                          cAAGGCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGt               >  1:550139/1‑93 (MQ=255)

                                                                           aaGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGaagt              >  1:3719422/1‑90 (MQ=255)

                                                                              gCGTGTTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCa           <  1:2351419/93‑1 (MQ=255)

                                                                                 tgtTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCatat        >  1:5332570/1‑93 (MQ=255)

                                                                                 tgtTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCatat        <  2:428670/93‑1 (MQ=255)

                                                                                 tgtTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCatat        >  1:1893917/1‑93 (MQ=255)

                                                                                  gtTGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATAtt       <  2:5514438/93‑1 (MQ=255)

                                                                                    tGGGTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATAttgt     <  1:4956609/93‑1 (MQ=255)

                                                                                     gggTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATAttgtt    <  1:3501428/93‑1 (MQ=255)

                                                                                     gggTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATAttgtt    <  2:4041373/93‑1 (MQ=255)

                                                                                      ggTAGCGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATAttgttg   >  1:2237221/1‑93 (MQ=255)

TGGACTGGACGCTGCTGCTGGTGTTTATGGCGATGTTTATCGACGTCCATTTACTGACCCAGCTTCCAGCGTTGCAAGGCGTGTTGGGTAACGTGAGTCATCTATCTGAACCCGGGTTATGGTTAACGGCAATCGGTTTATCGCAGGTGATCAGTAACGTGCCGAGTACCATATTGTTGC  >  NC_000913/853600‑853779

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