BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,108,438	C→A	noncoding (56/76 nt)	pheV →	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,108,438	0	C	A	99.5%	666.1 / NA	187	noncoding (56/76 nt)	pheV	tRNA‑Phe
Reads supporting (aligned to +/- strand): ref base C (0/1); new base A (57/129); total (57/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAA  >  NC_000913/3108348‑3108527

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                         <  1:5211392/93‑1 (MQ=255)

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                         <  1:1999056/93‑1 (MQ=255)

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                         <  1:1973596/93‑1 (MQ=255)

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                         <  2:5323026/93‑1 (MQ=255)

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                         <  1:928393/93‑1 (MQ=255)

 gagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAt                                                                                        <  1:4677859/93‑1 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCTCACCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                      <  2:2645782/93‑1 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                      <  1:4062913/93‑1 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                      <  1:2461484/93‑1 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                      <  1:4666010/93‑1 (MQ=255)

      gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGa                                                                                   <  1:1989845/93‑1 (MQ=255)

        ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGt                                                                                 >  1:5336358/1‑93 (MQ=255)

        ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGt                                                                                 <  2:2739646/93‑1 (MQ=255)

        agCGGATCAGGTTTCATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGt                                                                                 <  2:596446/92‑1 (MQ=11)

         gCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTc                                                                                <  1:4147667/93‑1 (MQ=255)

          cGAATCCAGTTTAAAGCCCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTcc                                                                               <  2:1098705/93‑1 (MQ=11)

           gAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                              <  2:452700/93‑1 (MQ=37)

           gAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                              >  1:3070012/1‑93 (MQ=37)

           gAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                              >  2:1183232/1‑93 (MQ=37)

            gaTCAGGTTTAATGCGCACCGTTGCCCGGATAGCTCAGTCGTGAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                             <  2:2344912/92‑1 (MQ=11)

            aaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                             >  2:1803567/1‑93 (MQ=25)

            aaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                             <  1:1340267/93‑1 (MQ=25)

              tCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGc                                                                           >  1:2962861/1‑93 (MQ=18)

                aGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcac                                                                         >  1:5121401/1‑93 (MQ=14)

                aGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcac                                                                         >  2:6067300/1‑93 (MQ=14)

                    ttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACt                                                                     >  1:3368551/1‑93 (MQ=17)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   >  1:4565907/1‑93 (MQ=21)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   <  1:752168/93‑1 (MQ=21)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   <  1:6051698/93‑1 (MQ=21)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   <  2:3331136/93‑1 (MQ=21)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   >  2:3939578/1‑93 (MQ=21)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   <  2:5230907/93‑1 (MQ=21)

                      aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                                   >  2:5570313/1‑93 (MQ=21)

                      aaTGCGCCACGTTGCCGAGATAGCTCAGTCGGTAGAGCAGTGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCCCCACTaa                                                                   <  2:1060811/93‑1 (MQ=11)

                        tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAtt                                                                 <  2:4402649/93‑1 (MQ=255)

                        tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAtt                                                                 <  1:1550514/93‑1 (MQ=37)

                           gcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                                              >  2:4497583/1‑93 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCTAGTCCGGGCACCACTAATTCTTa                                                             >  1:5633969/1‑93 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                             >  2:3367302/1‑93 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                             <  2:4850526/93‑1 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                             <  1:2296604/93‑1 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                             <  1:724759/93‑1 (MQ=255)

                            ccccGTTGCCCGAATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                             <  2:819768/93‑1 (MQ=255)

                             cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCTCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                            <  1:3767071/93‑1 (MQ=255)

                             cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                            <  2:2247454/93‑1 (MQ=255)

                             cccGGTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                            <  1:5339971/93‑1 (MQ=255)

                              ccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                                           >  2:4739913/1‑93 (MQ=255)

                                gTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                                         <  2:2849713/93‑1 (MQ=255)

                                  tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                                       <  2:6138095/93‑1 (MQ=255)

                                  tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                                       <  2:4564973/93‑1 (MQ=255)

                                  tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGCCACCACTAATTCTTAAGAAcc                                                       <  1:4013885/93‑1 (MQ=255)

                                   gCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAccc                                                      <  1:1171155/93‑1 (MQ=255)

                                    cccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCg                                                     >  1:3773066/1‑93 (MQ=255)

                                     ccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGc                                                    <  1:2315862/93‑1 (MQ=255)

                                       ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                                                  >  2:3685208/1‑93 (MQ=255)

                                         aTAGCTCAGTCGGTAGAGCAGGGGATTTAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                                                >  1:1598056/1‑93 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                              <  2:1408295/93‑1 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                              <  2:1860303/93‑1 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                              <  2:1020141/93‑1 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                              <  2:2613319/93‑1 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGACTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                              <  2:5926721/93‑1 (MQ=255)

                                            gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAg                                             <  2:5527202/93‑1 (MQ=255)

                                              tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACCAGGc                                           >  1:4410863/1‑93 (MQ=255)

                                              tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                                           <  2:4415479/93‑1 (MQ=255)

                                                  tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                       <  2:3093118/93‑1 (MQ=255)

                                                  tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                       <  2:5766437/93‑1 (MQ=255)

                                                   cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                                      >  1:4883829/1‑93 (MQ=255)

                                                    ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttt                                     <  2:200296/93‑1 (MQ=255)

                                                     gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtttt                                    <  2:2543645/93‑1 (MQ=255)

                                                      tAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                                   <  2:4020784/93‑1 (MQ=255)

                                                      gagagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                                   <  1:2485124/92‑1 (MQ=255)

                                                           cAGGGGATTGAAAACCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt                              <  1:5588642/93‑1 (MQ=255)

                                                            aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                             <  2:3270434/93‑1 (MQ=255)

                                                            aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                             <  2:348385/93‑1 (MQ=255)

                                                              gggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg                           <  1:1771310/93‑1 (MQ=255)

                                                              gggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg                           <  2:4193416/93‑1 (MQ=255)

                                                                gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGAt                         <  2:4466499/93‑1 (MQ=255)

                                                                gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGAt                         >  1:189984/1‑93 (MQ=255)

                                                                gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGAt                         <  2:2664671/93‑1 (MQ=255)

                                                                gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGAt                         <  1:3605434/93‑1 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAGGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                        >  2:728706/1‑93 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                        <  2:1470409/93‑1 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGGTTTTGCTTTTCGATc                        >  1:5436749/1‑93 (MQ=255)

                                                                  ttGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                       >  2:3833771/1‑93 (MQ=255)

                                                                     aaaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAc                    <  1:4140816/93‑1 (MQ=255)

                                                                       aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACaa                  <  2:5013220/93‑1 (MQ=255)

                                                                        aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAAt                 <  2:2963912/93‑1 (MQ=255)

                                                                         tCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATa                >  1:3221916/1‑93 (MQ=255)

                                                                         tCCCCGTGTCCTTGGTTAGATTCCGAGTCCAGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATa                >  1:5855302/1‑93 (MQ=255)

                                                                           cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc              <  2:160707/93‑1 (MQ=255)

                                                                           cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc              <  2:1564817/93‑1 (MQ=255)

                                                                             cGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCt            <  2:216876/93‑1 (MQ=255)

                                                                             cGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGGTTTTGCTTTTGGATCTGACAATAACCt            >  1:3510825/1‑93 (MQ=255)

                                                                               tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc          >  1:4768024/1‑93 (MQ=255)

                                                                               tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc          <  1:1094190/93‑1 (MQ=255)

                                                                               tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGGTTTTGATTTTGGATATGACAATAACCTTc          >  1:3980548/1‑93 (MQ=38)

                                                                                ttCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGATCCCGCCAACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCa         <  1:1522766/92‑1 (MQ=39)

                                                                                 tCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc        <  2:4094531/93‑1 (MQ=255)

                                                                                    ttGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGTCAATAACCTTCACGaa     >  2:5497359/1‑93 (MQ=255)

                                                                                    ctGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaa     <  2:1180256/92‑1 (MQ=255)

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAA  >  NC_000913/3108348‑3108527

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