BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,149	G→A	G75S (GGC→AGC)	yihI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,149	0	G	A	99.7%	1210.3 / NA	320	G75S (GGC→AGC)	yihI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (162/157); minor base T (0/1); total (162/158)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

cgcgCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGc                                                                                         >  2:1167376/1‑93 (MQ=255)

  cgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGt                                                                                       <  2:354189/93‑1 (MQ=255)

   gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                                      <  2:1380318/93‑1 (MQ=255)

   gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                                      >  1:2352116/1‑93 (MQ=255)

     aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGAc                                                                                    <  2:3289546/93‑1 (MQ=255)

      gcgggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACt                                                                                   >  2:2643562/1‑93 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 >  1:3829008/1‑93 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 <  2:1646106/93‑1 (MQ=255)

           cggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaa                                                                              >  1:6325280/1‑93 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            >  2:761117/1‑93 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            <  2:1552035/93‑1 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            >  1:2310268/1‑93 (MQ=255)

              aaaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGt                                                                           <  1:467052/93‑1 (MQ=255)

              aaaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGt                                                                           <  2:6626998/93‑1 (MQ=255)

              aaaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGt                                                                           >  2:6329358/1‑93 (MQ=255)

              aaaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGt                                                                           >  2:6292079/1‑93 (MQ=255)

               aaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTc                                                                          <  1:2755575/93‑1 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         <  2:2136681/93‑1 (MQ=255)

                 caccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCAc                                                                        >  2:109535/1‑93 (MQ=255)

                     acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                                    >  1:6410130/1‑93 (MQ=255)

                      cGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAAc                                                                   <  1:3144041/93‑1 (MQ=255)

                      cGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAAc                                                                   <  1:5611544/93‑1 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  <  2:5168132/93‑1 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  >  2:2045225/1‑93 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  <  1:5985165/93‑1 (MQ=255)

                        tCAGGCAGCAAAGGCCAGAGCGCACCAAAAGATCCACGCATTTGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                                 >  1:6522184/1‑93 (MQ=255)

                        tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                                 >  1:4004560/1‑93 (MQ=255)

                          aGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGca                                                               >  2:4779551/1‑93 (MQ=255)

                              agcaAAGGCCAGAACGCACCAAAAGATCCACGGATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaaa                                                           >  2:3931414/1‑93 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          >  2:2303954/1‑93 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          >  1:1170094/1‑93 (MQ=255)

                                caAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAcc                                                         >  1:2646572/1‑93 (MQ=255)

                                caAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAcc                                                         <  2:6070811/93‑1 (MQ=255)

                                caAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAcc                                                         <  1:3876372/93‑1 (MQ=255)

                                 aaaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCg                                                        <  2:4238504/93‑1 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      >  1:2396077/1‑93 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      >  1:825898/1‑93 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      >  1:5306567/1‑93 (MQ=255)

                                    ggCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAag                                                     >  1:6725764/1‑93 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    >  2:657931/1‑93 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    <  1:4968501/93‑1 (MQ=255)

                                      ccAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAagag                                                   <  2:4343139/93‑1 (MQ=255)

                                        aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                                 >  1:762055/1‑93 (MQ=255)

                                        aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCACTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                                 <  2:6452251/93‑1 (MQ=255)

                                         gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                                >  1:4334530/1‑93 (MQ=255)

                                         gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                                <  1:3590998/93‑1 (MQ=255)

                                           aCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgaga                                              <  2:3082524/93‑1 (MQ=255)

                                            cGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaa                                             <  1:6152410/93‑1 (MQ=255)

                                            cGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaa                                             <  2:2428129/93‑1 (MQ=255)

                                            cGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaa                                             <  1:1727285/93‑1 (MQ=255)

                                             gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                            >  2:4148594/1‑93 (MQ=255)

                                             gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                            >  1:4816535/1‑93 (MQ=255)

                                              cACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAc                                           <  1:2873512/93‑1 (MQ=255)

                                               aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                          >  1:1386186/1‑93 (MQ=255)

                                                ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                         <  2:2554416/93‑1 (MQ=255)

                                                ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                         >  1:6105626/1‑93 (MQ=255)

                                                ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGCGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcca                                         >  2:6233495/1‑92 (MQ=255)

                                                 cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                                        <  1:4272016/93‑1 (MQ=255)

                                                 cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                                        <  2:1769225/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAGCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  2:1240743/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  2:532829/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  2:4479995/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  1:449279/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCCGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       >  1:5710983/1‑93 (MQ=255)

                                                   aaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                                      >  1:1818814/1‑93 (MQ=255)

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                                     <  1:414807/93‑1 (MQ=255)

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                                     <  2:3445477/93‑1 (MQ=255)

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                                     >  1:1892232/1‑93 (MQ=255)

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCACTGAGCGTGACTGAAAAAGTGACCAAACAGAACAAACCGAAGAGAGAGAAACCTATGc                                     >  1:3206369/1‑93 (MQ=255)

                                                     aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACACACCGAAGAGTGAGAAACCTATGCt                                    >  1:329598/1‑93 (MQ=255)

                                                     aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCt                                    <  1:2085407/93‑1 (MQ=255)

                                                       aTCCACGTATTGGCAGTAAACCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                                  <  1:1235106/93‑1 (MQ=255)

                                                       aTCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                                  <  2:5368339/93‑1 (MQ=255)

                                                        tCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTc                                 >  2:1593085/1‑93 (MQ=255)

                                                         ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                                <  2:5915230/93‑1 (MQ=255)

                                                            cGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCg                             <  2:4841346/93‑1 (MQ=255)

                                                              tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                           <  1:353520/93‑1 (MQ=255)

                                                               aTTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACGAAACAGCACAACCCGAAGAGTGAGAAACCTATGCTTTCACCGCAg                          <  1:5101280/93‑1 (MQ=255)

                                                                ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                         <  2:3806721/93‑1 (MQ=255)

                                                                  ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg                       <  2:5561678/93‑1 (MQ=255)

                                                                     aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg                    >  2:2723532/1‑93 (MQ=255)

                                                                     aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg                    <  2:5252110/93‑1 (MQ=255)

                                                                     aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg                    <  1:5806882/93‑1 (MQ=255)

                                                                      gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCATACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt                   >  2:3558570/1‑93 (MQ=255)

                                                                      gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt                   <  1:5446373/93‑1 (MQ=255)

                                                                      gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt                   <  1:3402847/93‑1 (MQ=255)

                                                                      gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt                   <  1:3024185/93‑1 (MQ=255)

                                                                      atAAACCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt                   <  2:508574/92‑1 (MQ=255)

                                                                        aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg                 >  1:3061369/1‑93 (MQ=255)

                                                                         aaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTgg                >  1:2461192/1‑93 (MQ=255)

                                                                          aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa               >  2:153810/1‑93 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  2:1815345/93‑1 (MQ=255)

                                                                            ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt             <  1:4043972/93‑1 (MQ=255)

                                                                            ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt             >  1:4876728/1‑93 (MQ=255)

                                                                               ctgttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTAc          <  1:1007072/90‑1 (MQ=255)

                                                                                 gttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTg        >  1:5281617/2‑93 (MQ=255)

                                                                                     cATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaaa    <  2:3665447/93‑1 (MQ=255)

                                                                                      aTTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAAc   <  1:3729204/93‑1 (MQ=255)

                                                                                      aTTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAAc   >  2:2858603/1‑93 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACg  <  1:1529281/93‑1 (MQ=255)

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

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