BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,751,393	G→A	T154T (ACC→ACT)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,751,393	0	G	A	98.6%	507.8 / ‑3.8	144	T154T (ACC→ACT)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (75/67); total (77/67)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 146 total aligned reads displayed.

GCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTAAT  >  NC_000913/1751303‑1751482

gccgcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGt                                                                                         <  1:2966639/93‑1 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                      >  1:891518/1‑93 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                      >  2:1958349/1‑93 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                      >  2:2122683/1‑93 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                      >  2:4442770/1‑93 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                      <  1:2018498/93‑1 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                      >  2:4442769/1‑93 (MQ=255)

    cTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAAt                                                                                     <  2:1664491/93‑1 (MQ=255)

       tgGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGa                                                                                  >  1:3916965/1‑93 (MQ=255)

        ggTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGAc                                                                                 >  1:2020205/1‑93 (MQ=255)

         gTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACa                                                                                >  2:4602950/1‑93 (MQ=255)

          ttACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACaa                                                                               >  1:1640585/1‑93 (MQ=255)

             cGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAtt                                                                            <  1:1538683/93‑1 (MQ=255)

             cGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAtt                                                                            <  1:2843431/93‑1 (MQ=255)

              ggCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAttt                                                                           <  1:1988658/93‑1 (MQ=255)

              ggCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAttt                                                                           <  2:3903359/93‑1 (MQ=255)

              ggCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAttt                                                                           <  2:1138291/93‑1 (MQ=255)

                cTATTCAACGTGCGAGTGAGATGAACAAGGTTTGCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttct                                                                         <  2:2276297/93‑1 (MQ=255)

                 tATTCAACGTGCCAGCGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGACGCCACACAGGTTTCCGGACTGGTCAATCGACAAATttctc                                                                        >  2:2862983/1‑92 (MQ=255)

                 tATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttctt                                                                        >  2:3652139/1‑93 (MQ=255)

                      aaCGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtc                                                                   >  2:2311528/1‑93 (MQ=255)

                      aaCGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtc                                                                   >  1:4395403/1‑93 (MQ=255)

                        cgtGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtcgt                                                                 <  1:1354311/91‑1 (MQ=255)

                         gtGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGtt                                                                >  2:2313079/2‑93 (MQ=255)

                           gCCAGAGAGATGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgc                                                              <  2:466153/93‑1 (MQ=255)

                              agagagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcgcg                                                           <  1:2241063/93‑1 (MQ=255)

                              agagagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcgcg                                                           >  1:4008634/1‑93 (MQ=255)

                               gagagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcgcgc                                                          >  2:2850233/1‑93 (MQ=255)

                               gagagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcgcgc                                                          >  2:4717732/1‑93 (MQ=255)

                                 gagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGcgcgcgc                                                        >  2:1960028/1‑92 (MQ=255)

                                   gagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTcc                                                      >  1:646673/1‑93 (MQ=255)

                                   gagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTcc                                                      <  2:3531638/93‑1 (MQ=255)

                                   gagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTcc                                                      >  1:752397/1‑93 (MQ=255)

                                     ggAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCt                                                    <  1:1254191/93‑1 (MQ=255)

                                       aaCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCttt                                                  >  1:2543624/1‑93 (MQ=255)

                                       aaCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCttt                                                  <  1:3185624/93‑1 (MQ=255)

                                        aCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCtttt                                                 >  1:2478577/1‑93 (MQ=255)

                                        aCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCtttt                                                 >  2:430846/1‑93 (MQ=255)

                                        aCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCtttt                                                 >  1:299655/1‑93 (MQ=255)

                                         cAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTc                                                >  1:1498007/1‑93 (MQ=255)

                                          aaGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTcc                                               >  2:3994679/1‑93 (MQ=255)

                                          aaGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGGTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTcc                                               >  1:1033702/1‑93 (MQ=255)

                                           aGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTccc                                              >  2:898302/1‑93 (MQ=255)

                                            ggTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTcccc                                             <  2:4687147/93‑1 (MQ=255)

                                            ggTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTcccc                                             >  2:2543873/1‑93 (MQ=255)

                                             gTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCa                                            <  1:1647256/93‑1 (MQ=255)

                                             gTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCa                                            <  2:4560062/93‑1 (MQ=255)

                                             gTTTCCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCa                                            <  2:2002008/93‑1 (MQ=255)

                                              ttttCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCAt                                           >  2:2201089/1‑93 (MQ=255)

                                              ttttCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCAt                                           >  2:3937415/1‑93 (MQ=255)

                                                tcccccAGCCTGACCAAAAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATaa                                         <  1:4094051/91‑1 (MQ=255)

                                                 tCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAg                                        <  1:3670667/93‑1 (MQ=255)

                                                  ccccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGa                                       >  2:1765212/1‑93 (MQ=255)

                                                   cccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaa                                      >  1:3848795/1‑93 (MQ=255)

                                                   cccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaa                                      >  1:1733835/1‑93 (MQ=255)

                                                    ccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaaa                                     <  1:4134103/93‑1 (MQ=255)

                                                     cAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAAt                                    <  2:3905417/93‑1 (MQ=255)

                                                      aGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATc                                   <  1:2334912/93‑1 (MQ=255)

                                                      aGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATc                                   >  1:2284322/1‑93 (MQ=255)

                                                      aGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATc                                   <  1:4035957/93‑1 (MQ=255)

                                                      aGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATc                                   <  1:4130190/93‑1 (MQ=255)

                                                       gCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCg                                  >  2:3112159/1‑93 (MQ=255)

                                                       gCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCg                                  <  1:2739436/93‑1 (MQ=255)

                                                         cTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTGCCTTTTCCCCATAAGAAATCGGc                                >  2:3232806/1‑93 (MQ=255)

                                                          tGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCt                               >  1:4143163/1‑93 (MQ=255)

                                                           ggccAAGAGCCGCCACACAGGTTTCCGGACTAGTCAAGCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtt                              <  1:1162787/91‑1 (MQ=255)

                                                           gACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtt                              <  1:3144507/93‑1 (MQ=255)

                                                           gACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtt                              >  2:3625255/1‑93 (MQ=255)

                                                            aCCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCttt                             <  1:3871910/93‑1 (MQ=255)

                                                             ccAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtttt                            <  1:2825060/93‑1 (MQ=255)

                                                             ccAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtttt                            >  2:4204745/1‑93 (MQ=255)

                                                             ccAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtttt                            <  2:680080/93‑1 (MQ=255)

                                                              cAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCttttt                           >  1:3105597/1‑93 (MQ=255)

                                                               aaTAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTc                          <  1:4388777/93‑1 (MQ=255)

                                                                 tAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCa                        <  2:266433/93‑1 (MQ=255)

                                                                 tAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCa                        >  1:520464/1‑93 (MQ=255)

                                                                  aGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAg                       >  1:1741187/1‑93 (MQ=255)

                                                                  aGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAg                       <  1:4312492/93‑1 (MQ=255)

                                                                    ccgccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGc                     <  2:3046771/93‑1 (MQ=255)

                                                                    ccgccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGc                     <  2:478496/93‑1 (MQ=255)

                                                                    ccgccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGc                     >  2:1376056/1‑93 (MQ=255)

                                                                       ccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTaa                  <  2:4351490/93‑1 (MQ=255)

                                                                        cacacaGGTTTCCGGACTAGTCCATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAAc                 <  1:3239454/93‑1 (MQ=255)

                                                                        cacacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAAc                 >  2:2744887/1‑93 (MQ=255)

                                                                        cacacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAAc                 >  2:2195427/1‑93 (MQ=255)

                                                                        cacacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAAc                 >  2:4179224/1‑93 (MQ=255)

                                                                        cacacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAAc                 <  1:8365/93‑1 (MQ=255)

                                                                         acacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACt                >  1:3562059/1‑93 (MQ=255)

                                                                         acacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACt                >  1:4694042/1‑93 (MQ=255)

                                                                          cacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACtt               <  1:1145051/93‑1 (MQ=255)

                                                                           acaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACttt              >  2:3088865/1‑93 (MQ=255)

                                                                            caGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTg             <  1:2407261/93‑1 (MQ=255)

                                                                             aGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGt            >  1:4419039/1‑93 (MQ=255)

                                                                              ggTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTc           <  2:3595137/93‑1 (MQ=255)

                                                                              ggTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTc           <  2:3277039/93‑1 (MQ=255)

                                                                              ggTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTc           <  2:259714/93‑1 (MQ=255)

                                                                               gTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCa          >  2:4316972/1‑93 (MQ=255)

                                                                                      gACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTaa   <  2:2087507/93‑1 (MQ=255)

                                                                                      gACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTaa   <  2:893882/93‑1 (MQ=255)

                                                                                       aCTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTAAt  <  2:646673/93‑1 (MQ=255)

GCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTAAT  >  NC_000913/1751303‑1751482

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