BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,572,272	G→A	G190E (GGG→GAG)	yjiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,572,272	0	G	A	99.5%	1541.5 / ‑4.1	433	G190E (GGG→GAG)	yjiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (205/226); total (206/227)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 447 total aligned reads displayed.

TTACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCAAA  >  NC_000913/4572182‑4572361

ttACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAgaga                                                                                         <  1:1330322/93‑1 (MQ=255)

ttACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAgaga                                                                                         >  2:2362266/1‑93 (MQ=255)

  aCCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACa                                                                                       <  2:5233895/93‑1 (MQ=255)

    cAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGc                                                                                     >  1:4967732/1‑93 (MQ=255)

    cAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGc                                                                                     <  2:3163467/93‑1 (MQ=255)

        gatgaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGa                                                                                 <  2:1508611/93‑1 (MQ=255)

         atgaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGaa                                                                                <  1:7340341/93‑1 (MQ=255)

          tgaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAAt                                                                               >  2:1294119/1‑93 (MQ=255)

          tgaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAAt                                                                               <  2:6528026/93‑1 (MQ=255)

          tgaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAAt                                                                               <  1:2788620/93‑1 (MQ=255)

          agaAGGCAAAGAGGGCGTCGGCCCAGATTCGTTGCCAGGACGAGGTCGTGGACTCATCACCTTTACGGACAGTCGACAAGAGACAGCCCGAAt                                                                               <  1:5667310/92‑1 (MQ=255)

           gaAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATg                                                                              <  2:8328390/93‑1 (MQ=255)

            aaGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATgg                                                                             >  1:3989189/1‑93 (MQ=255)

            aaGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATgg                                                                             >  2:5819979/1‑93 (MQ=255)

              ggCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCg                                                                           <  2:1623028/93‑1 (MQ=255)

              ggCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCg                                                                           <  2:2956147/93‑1 (MQ=255)

                cAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGt                                                                         <  1:5452110/93‑1 (MQ=255)

                 aaaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTg                                                                        >  1:3416249/1‑93 (MQ=255)

                  aaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTgg                                                                       >  2:3438347/1‑92 (MQ=255)

                  aaGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGc                                                                       >  2:2077102/1‑93 (MQ=255)

                   agagGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCg                                                                      <  2:7158865/93‑1 (MQ=255)

                   agagGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCg                                                                      <  1:7020899/93‑1 (MQ=255)

                     agGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCggt                                                                    >  1:6193584/1‑91 (MQ=255)

                      gggCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTAt                                                                   <  1:186035/93‑1 (MQ=255)

                      gggCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTAt                                                                   >  1:1976031/1‑93 (MQ=255)

                       ggCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGCGCGTATg                                                                  >  1:2526105/1‑93 (MQ=255)

                         cGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgca                                                                <  2:8058723/93‑1 (MQ=255)

                          gTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgcag                                                               >  2:8774806/1‑93 (MQ=255)

                           tCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGCCAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgcaga                                                              >  1:6286421/1‑92 (MQ=255)

                           tCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgcagc                                                              >  1:2977148/1‑93 (MQ=255)

                            cGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATgcagca                                                             >  2:6629303/1‑93 (MQ=255)

                             ggACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaa                                                            >  2:7930911/1‑93 (MQ=255)

                             ggACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaa                                                            >  1:6316787/1‑93 (MQ=255)

                              gACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaag                                                           >  1:768427/1‑93 (MQ=255)

                              gACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaag                                                           <  2:4541864/93‑1 (MQ=255)

                                ccAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaagaa                                                         <  2:3108079/93‑1 (MQ=255)

                                ccAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCaagaa                                                         <  1:4288960/93‑1 (MQ=255)

                                  aGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTGCAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGc                                                       <  2:4064319/93‑1 (MQ=255)

                                  aGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGc                                                       <  2:6755016/93‑1 (MQ=255)

                                  aGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGc                                                       >  2:4710221/1‑93 (MQ=255)

                                  aGATTCGTTGCCAGGACGACGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGc                                                       <  1:8570976/93‑1 (MQ=255)

                                   gATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACGAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCa                                                      >  1:6677897/1‑93 (MQ=255)

                                   gATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCa                                                      <  2:444362/93‑1 (MQ=255)

                                   gATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCa                                                      >  1:7998625/1‑93 (MQ=255)

                                     ttCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGa                                                    <  2:3719711/93‑1 (MQ=255)

                                         ttGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGc                                                >  2:5427965/1‑93 (MQ=255)

                                          tGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCa                                               <  1:1740414/93‑1 (MQ=255)

                                           gCCAGGACGATGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCTGTGCGTATGCAGCAAGAAGCAGAACGCAg                                              <  2:5636364/93‑1 (MQ=255)

                                            ccAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGt                                             <  2:7451891/93‑1 (MQ=255)

                                            ccAGGACGAGGTCGTCAACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGt                                             <  1:2335769/93‑1 (MQ=255)

                                             cAGGACGAGGTTGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            <  1:567692/93‑1 (MQ=255)

                                             cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            >  2:3202709/1‑93 (MQ=255)

                                             cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            <  1:1018751/93‑1 (MQ=255)

                                             cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            <  2:4995077/93‑1 (MQ=255)

                                             cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            <  1:6819274/93‑1 (MQ=255)

                                             cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            <  1:4682275/93‑1 (MQ=255)

                                             cAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTc                                            >  1:1565494/1‑93 (MQ=255)

                                               ggACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGc                                          <  1:5802060/93‑1 (MQ=255)

                                                gACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCt                                         <  2:4130961/93‑1 (MQ=255)

                                                gACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCt                                         <  2:5353208/93‑1 (MQ=255)

                                                 aCGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCTCtt                                        <  2:1376690/93‑1 (MQ=255)

                                                  agAGGTCGTCGAATCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTa                                       <  2:6334543/92‑1 (MQ=255)

                                                    aGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcg                                     >  2:4221482/1‑93 (MQ=255)

                                                     ggTCGTCGACTCATCACCTTTACAGACATTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcgc                                    <  2:4856099/93‑1 (MQ=255)

                                                     ggTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcgc                                    <  1:1183402/93‑1 (MQ=255)

                                                     ggTCGTCGACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTAcgc                                    >  1:5990094/1‑93 (MQ=255)

                                                        cgtcgACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGa                                 <  1:8627686/93‑1 (MQ=255)

                                                         gtcgACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGaa                                >  2:6164503/1‑93 (MQ=255)

                                                          tcgACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAg                               >  2:2684500/1‑93 (MQ=255)

                                                            gACTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCg                             >  1:3263270/1‑93 (MQ=255)

                                                             aCTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGt                            <  1:3548107/93‑1 (MQ=255)

                                                              cTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTa                           <  1:4560567/93‑1 (MQ=255)

                                                              cTCATCACCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTa                           <  1:7509872/93‑1 (MQ=255)

                                                               tcatcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAg                          >  2:1050300/1‑93 (MQ=255)

                                                                catcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                         <  1:4067359/93‑1 (MQ=255)

                                                                catcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                         <  1:8202163/93‑1 (MQ=255)

                                                                catcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGt                         >  1:5254528/1‑93 (MQ=255)

                                                                 atcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGtg                        >  1:5531815/1‑92 (MQ=255)

                                                                  tcaCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCg                       >  2:2444335/1‑93 (MQ=255)

                                                                    aCCTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGaa                     <  1:491474/93‑1 (MQ=255)

                                                                     ccTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGaaa                    >  1:2636015/1‑93 (MQ=255)

                                                                      cTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAAt                   <  1:3188862/93‑1 (MQ=255)

                                                                      cTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAAt                   >  1:7179550/1‑93 (MQ=255)

                                                                      cTTTACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAAt                   >  2:3296011/1‑93 (MQ=255)

                                                                        ttACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGTAGCAGAAAGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtc                 >  1:8408499/1‑93 (MQ=255)

                                                                         tACAGACAGTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATtct                <  2:7688104/93‑1 (MQ=255)

                                                                            agacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAg             >  2:3069380/1‑93 (MQ=255)

                                                                            agacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAg             <  1:8705656/93‑1 (MQ=255)

                                                                             gacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGc            >  2:5090301/1‑93 (MQ=255)

                                                                             gacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGc            >  1:5153505/1‑93 (MQ=255)

                                                                             gacagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGc            >  2:6936821/1‑93 (MQ=255)

                                                                               cagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTg          >  1:4826253/1‑93 (MQ=255)

                                                                               cagTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTg          <  1:3771045/93‑1 (MQ=255)

                                                                                agTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTgg         >  1:641888/1‑93 (MQ=255)

                                                                                agTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTgg         <  1:7910241/93‑1 (MQ=255)

                                                                                 gTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGc        >  2:7779492/1‑93 (MQ=255)

                                                                                 gTCGACAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGc        <  2:2004355/93‑1 (MQ=255)

                                                                                      cAAGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCaa   >  1:2157412/1‑93 (MQ=255)

                                                                                      cAAGAGACAGCAAGAATGGCGGTGCGTATGCATCAAGAAGCAGAACGAAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCaa   <  2:1176076/93‑1 (MQ=255)

                                                                                       aaGAGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCaaa  <  2:3858006/93‑1 (MQ=255)

TTACCAGTGATGAAGGCAAAGAGGGCGTCGGACCAGATTCGTTGCCAGGACGAGGTCGTCGACTCATCACCTTTACAGACAGTCGACAAGGGACAGCCCGAATGGCGGTGCGTATGCAGCAAGAAGCAGAACGCAGTCGCTTACGCGGAAGCGTAGTCGAAATTCTCAGCTGGCATCAAA  >  NC_000913/4572182‑4572361

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