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.2%	410.8 / ‑6.1	121	noncoding (56/76 nt)	pheV	tRNA‑Phe
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (52/68); minor base T (0/1); total (52/69)
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 123 total aligned reads displayed.

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAA  >  NC_000913/3108348‑3108526

  agacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAtt                                                                                      <  1:3562179/93‑1 (MQ=255)

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

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                     >  2:3553884/1‑93 (MQ=255)

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

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

       agGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAg                                                                                 >  2:861209/1‑93 (MQ=255)

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

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

          cGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTcc                                                                              <  2:707585/93‑1 (MQ=255)

           gAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                             <  2:2025162/93‑1 (MQ=255)

            aaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGTGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                            >  1:2658852/1‑93 (MQ=255)

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

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

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

                aGGTTTAATGCACCCCAATGCCCGAATCGCTCCGTCAATCGCGCAGGGATTAGAGATTCCCCGTGTCCTTGGTTTGGTTCCGTGTACGGGcac                                                                        <  1:3618708/93‑1 (MQ=2)

                   tttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                                     >  1:1751801/1‑93 (MQ=14)

                   tttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                                     >  1:912789/1‑93 (MQ=14)

                       aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAt                                                                 <  1:2310630/93‑1 (MQ=25)

                       aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAt                                                                 >  2:2310125/1‑93 (MQ=25)

                        tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGACGGGGATTGAAAATCCCCGTGTCCTTGGCTAGAATCCGAGTCCGGGCACCACTAAtt                                                                >  2:3244066/1‑93 (MQ=11)

                          cgcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCt                                                              <  2:62207/93‑1 (MQ=255)

                           gcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGCTTGAAAATCCCCGCGTCCTTGGTTAGATTCCGAGTCCGGGCAGCACTAATTCtt                                                             >  1:1994219/1‑93 (MQ=14)

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

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

                             cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                           >  1:2921072/1‑93 (MQ=255)

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

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                         >  2:3235978/1‑93 (MQ=255)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                         <  1:3535473/93‑1 (MQ=255)

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

                                gTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                                        >  2:3500115/1‑93 (MQ=255)

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

                                 tttcccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAc                                                       <  2:3170913/90‑1 (MQ=255)

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

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

                                      cGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                                  <  2:238539/93‑1 (MQ=255)

                                       ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                                                 <  2:3308888/93‑1 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                >  2:4592676/1‑93 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                >  2:1847023/1‑93 (MQ=255)

                                          tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                                              >  1:477504/1‑93 (MQ=255)

                                          tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                                              <  1:4466574/93‑1 (MQ=255)

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

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

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTACTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                             >  2:3389105/1‑93 (MQ=255)

                                              tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                                          >  1:736864/1‑93 (MQ=255)

                                              tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                                          >  2:4088807/1‑93 (MQ=255)

                                              tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                                          >  1:4115047/1‑93 (MQ=255)

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

                                               cAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCg                                         <  1:335354/93‑1 (MQ=255)

                                                aGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCgg                                        <  1:1486171/93‑1 (MQ=255)

                                                 gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                                       >  1:2369785/1‑93 (MQ=255)

                                                  tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGACCCCGCCCACAAGGCGGGt                                      >  1:3467300/1‑93 (MQ=255)

                                                    ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttt                                    >  1:2703876/1‑93 (MQ=255)

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

                                                     gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCAGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtttt                                   <  2:1706503/93‑1 (MQ=255)

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

                                                      tAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGGCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                                  <  1:3861221/93‑1 (MQ=255)

                                                        gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                >  2:4467623/1‑93 (MQ=255)

                                                        gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                >  2:2806039/1‑93 (MQ=255)

                                                         agCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCt                               >  1:3911979/1‑93 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              >  1:1850738/1‑93 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              >  2:3323530/1‑93 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              >  2:3323528/1‑93 (MQ=255)

                                                           cAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt                             >  1:1949575/1‑93 (MQ=255)

                                                           cAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt                             >  2:3213875/1‑93 (MQ=255)

                                                            aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                            >  1:839459/1‑93 (MQ=255)

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

                                                             ggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                           <  1:1003063/93‑1 (MQ=255)

                                                             ggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                           <  1:2925169/93‑1 (MQ=255)

                                                             ggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                           <  1:477995/93‑1 (MQ=255)

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

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

                                                                 aTTGAAAATCCCCGTGTGCTTGGTTAGATTCCGAGTCCGTGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                       <  2:4257400/93‑1 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                       <  2:2877441/93‑1 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                       >  2:2141898/1‑93 (MQ=255)

                                                                   tGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                     <  1:769773/93‑1 (MQ=255)

                                                                   tGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                     <  2:4159482/93‑1 (MQ=255)

                                                                   tGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                     <  1:2408500/93‑1 (MQ=255)

                                                                    gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGa                    <  2:650745/93‑1 (MQ=255)

                                                                    gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGa                    <  2:2130584/93‑1 (MQ=255)

                                                                    gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGa                    <  2:3500976/93‑1 (MQ=255)

                                                                     aaaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAc                   >  2:1594083/1‑93 (MQ=255)

                                                                     aaaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGC‑TTTGGATCTGACa                  <  1:2757924/93‑1 (MQ=255)

                                                                      aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACa                  >  2:4248160/1‑93 (MQ=255)

                                                                       aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAACCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGActa                 >  1:406848/1‑91 (MQ=255)

                                                                        aTCCCCGTGTCCTTGGTTAGATTCCGAGTTCGGGCACCACTAATTCTTAAGTACCCGCCTACAAGGCGGTGTTGTGCCTGTGGATCTGACAAt                >  1:4208172/1‑93 (MQ=21)

                                                                        aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAAt                >  2:3210440/1‑93 (MQ=255)

                                                                         tCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATa               <  1:1185715/93‑1 (MQ=255)

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

                                                                          ccccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATaa              <  1:2701152/93‑1 (MQ=255)

                                                                            ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAcc            >  1:986596/1‑93 (MQ=255)

                                                                            ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAcc            <  2:411196/93‑1 (MQ=255)

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

                                                                              gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCtt          <  1:1026120/93‑1 (MQ=255)

                                                                              gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCtt          >  1:1150591/1‑93 (MQ=255)

                                                                              gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCtt          >  2:3912080/1‑93 (MQ=255)

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

                                                                                    ttGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaa    >  1:268692/1‑93 (MQ=255)

                                                                                    ttGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaa    >  1:4138923/1‑93 (MQ=255)

                                                                                    ttGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAAAAACCTTCACGaa    >  2:1005068/1‑93 (MQ=255)

                                                                                      ggTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaa  >  1:2820430/1‑93 (MQ=255)

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAA  >  NC_000913/3108348‑3108526

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