BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,822,652	T→C	Y171H (TAC→CAC)	spoU →	putative RNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,822,652	0	T	C	98.9%	611.8 / ‑3.3	188	Y171H (TAC→CAC)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (87/99); total (88/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 192 total aligned reads displayed.

ATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822562‑3822741

  gatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCt                                                                                       <  1:321571/93‑1 (MQ=255)

   atCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTg                                                                                      <  2:2022309/93‑1 (MQ=255)

    tCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGc                                                                                     >  2:1852440/1‑93 (MQ=255)

     cGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCg                                                                                    <  1:1510494/93‑1 (MQ=255)

      ggCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGt                                                                                   <  2:1943706/93‑1 (MQ=255)

      ggCATGGTGAGGTCGCTGATTGTTTCCGTTACCTCAGCCCTCATTCTTTACAAAACCCAGCGTCAGCGGCAAAATGAAGGCCTGCCCCTGCGt                                                                                   <  2:4348766/93‑1 (MQ=255)

        cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGa                                                                                 >  2:4037233/1‑93 (MQ=255)

        cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGa                                                                                 <  2:11573/93‑1 (MQ=255)

           ggTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaaaa                                                                              <  1:3857400/93‑1 (MQ=255)

             tGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACa                                                                            <  1:3016347/93‑1 (MQ=255)

                aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCATTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         <  2:4750148/93‑1 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        <  2:907968/93‑1 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        <  1:2693373/93‑1 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  1:3173175/1‑93 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  2:739194/1‑93 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  2:1154862/1‑93 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  2:761269/1‑93 (MQ=255)

                   cGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGt                                                                      <  2:3101928/93‑1 (MQ=255)

                   cGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGt                                                                      <  1:4087578/93‑1 (MQ=255)

                   cGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGt                                                                      <  2:2808166/93‑1 (MQ=255)

                    gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     <  1:1635428/93‑1 (MQ=255)

                    gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     >  1:3378912/1‑93 (MQ=255)

                      tGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGc                                                                   <  2:451914/93‑1 (MQ=255)

                       gAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGcc                                                                  >  2:3629237/1‑93 (MQ=255)

                       gAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGcc                                                                  <  1:135250/93‑1 (MQ=255)

                       gAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGcc                                                                  >  1:4557832/1‑93 (MQ=255)

                        aaTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCg                                                                 >  1:3931643/1‑93 (MQ=255)

                        aaTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCg                                                                 >  1:598941/1‑93 (MQ=255)

                          tGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGa                                                               <  2:4434100/93‑1 (MQ=255)

                           gTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCACAATGCAGGCATGCACCTGACTGAAAACAGCATGTTGCCGGaa                                                              >  2:3154883/1‑93 (MQ=255)

                           gTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGaa                                                              >  2:3584894/1‑93 (MQ=255)

                           gTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGaa                                                              <  1:4744979/93‑1 (MQ=255)

                            tttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAGATGCAGGCATGCACCTGCGTGCAAACAGCATGTTGCCGGAAg                                                             >  2:3262723/1‑93 (MQ=255)

                            tttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAg                                                             <  2:2299207/93‑1 (MQ=255)

                            tttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAg                                                             <  1:602085/93‑1 (MQ=255)

                             ttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGc                                                            >  1:445562/1‑93 (MQ=255)

                             ttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGc                                                            <  1:1101554/93‑1 (MQ=255)

                              tCCGTTGCCTCAGCCCTCATCAGTTACGCACCACGGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCa                                                           <  2:942604/93‑1 (MQ=255)

                               ccGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAACCAGCATGTTGCCGGAAGCag                                                          <  1:3617349/93‑1 (MQ=255)

                               ccGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCag                                                          >  2:3241391/1‑93 (MQ=255)

                                cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         <  2:1481916/93‑1 (MQ=255)

                                cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         <  1:4522156/93‑1 (MQ=255)

                                 gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        >  1:1576257/1‑93 (MQ=255)

                                 gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        >  1:2356636/1‑93 (MQ=255)

                                  ttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGc                                                       <  1:3135536/93‑1 (MQ=255)

                                    gCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaa                                                     <  1:402426/93‑1 (MQ=255)

                                      cTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaaca                                                   >  2:64339/1‑93 (MQ=255)

                                       tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  >  1:3250782/1‑93 (MQ=255)

                                       tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  >  1:3357467/1‑93 (MQ=255)

                                       tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  >  2:3012376/1‑93 (MQ=255)

                                         aGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACg                                                <  1:4053489/93‑1 (MQ=255)

                                           cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              >  2:1687528/1‑93 (MQ=255)

                                           cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              >  1:1050740/1‑93 (MQ=255)

                                            ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACACCGCCt                                             >  2:1376014/1‑93 (MQ=255)

                                            ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             >  1:2793679/1‑93 (MQ=255)

                                              tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAGGCAGAGACACAACGCCtgt                                           >  1:1458922/1‑93 (MQ=255)

                                              tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                           >  2:4083340/1‑93 (MQ=255)

                                               cATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgtt                                          >  2:4300667/1‑93 (MQ=255)

                                                aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                         <  1:2497682/93‑1 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                        >  1:447978/1‑93 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                        <  1:289486/93‑1 (MQ=255)

                                                 ctcgttACGAAGCCCAGCGGCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAATCAGAGCAACAACCCCtgttgt                                        <  2:1931769/89‑1 (MQ=255)

                                                   cTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGttt                                      <  2:178383/93‑1 (MQ=255)

                                                    tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGttta                                     >  1:3529532/1‑92 (MQ=255)

                                                    tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                     >  1:2725135/1‑93 (MQ=255)

                                                     ttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGa                                    <  1:1039858/93‑1 (MQ=255)

                                                     ttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGa                                    >  1:4136039/1‑93 (MQ=255)

                                                     ttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGa                                    <  1:2977156/93‑1 (MQ=255)

                                                      tACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGaa                                   <  2:4465380/93‑1 (MQ=255)

                                                       aCGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAg                                  <  2:4495641/93‑1 (MQ=255)

                                                        cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                 <  1:4327475/93‑1 (MQ=255)

                                                        cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                 >  1:2312122/1‑93 (MQ=255)

                                                        cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                 <  1:983187/93‑1 (MQ=255)

                                                           aGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcgg                              >  2:1291328/1‑93 (MQ=255)

                                                            gCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcggc                             >  1:2091805/1‑93 (MQ=255)

                                                            gCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcggc                             <  1:4239936/93‑1 (MQ=255)

                                                             cccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                            >  2:2742000/1‑93 (MQ=255)

                                                             cccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                            <  2:240630/93‑1 (MQ=255)

                                                             cccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                            <  1:3652034/93‑1 (MQ=255)

                                                              ccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                           >  1:1740944/1‑93 (MQ=255)

                                                              ccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                           >  1:2241136/1‑93 (MQ=255)

                                                                   gTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCgg                      >  2:2589038/1‑93 (MQ=255)

                                                                    tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  2:1895480/1‑93 (MQ=255)

                                                                    tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  1:3682855/1‑93 (MQ=255)

                                                                     cAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTg                    >  1:2718066/1‑93 (MQ=255)

                                                                      aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   <  1:71080/93‑1 (MQ=255)

                                                                      aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   <  1:755172/93‑1 (MQ=255)

                                                                      aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   >  2:2043432/1‑93 (MQ=255)

                                                                       gCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGACTATCCGGTGCt                  >  2:22159/1‑93 (MQ=255)

                                                                           cAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCg              <  1:3580818/93‑1 (MQ=255)

                                                                               aTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAg          >  2:3917956/1‑93 (MQ=255)

                                                                                tGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt         <  1:1993292/93‑1 (MQ=255)

                                                                                 gCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTc        <  2:2073814/93‑1 (MQ=255)

                                                                                   aGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGc      <  1:3438138/93‑1 (MQ=255)

                                                                                    ggCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCa     >  1:403467/1‑93 (MQ=255)

                                                                                    ggCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCa     >  2:3808739/1‑93 (MQ=255)

                                                                                     gCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaa    <  1:1264524/93‑1 (MQ=255)

                                                                                     gCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaa    <  2:816972/93‑1 (MQ=255)

                                                                                       aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  1:4406162/1‑93 (MQ=255)

                                                                                       aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  2:2642729/1‑93 (MQ=255)

ATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822562‑3822741

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