BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,845,852	T→C	T61A (ACC→GCC)	selD ←	selenophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,845,852	0	T	C	97.7%	395.9 / ‑1.0	129	T61A (ACC→GCC)	selD	selenophosphate synthase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (74/52); total (76/53)
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.48e-01

Only 100 of 133 total aligned reads displayed.

CGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAT  >  NC_000913/1845762‑1845941

cGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGgcgc                                                                                         >  2:4521593/1‑93 (MQ=255)

 gCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGcc                                                                                        >  1:4259820/1‑93 (MQ=255)

  cTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATGAAGAAGTCGGTGGTACTGATAACGCTGGCGCCa                                                                                       <  1:4066404/93‑1 (MQ=255)

   tGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAt                                                                                      >  2:1625241/1‑93 (MQ=255)

   tGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAt                                                                                      >  2:1004525/1‑93 (MQ=255)

        gCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGccc                                                                                 <  1:3376527/93‑1 (MQ=255)

          gTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                               >  1:2904218/1‑93 (MQ=255)

           ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              >  2:681491/1‑93 (MQ=255)

           ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              <  1:3868407/93‑1 (MQ=255)

           ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              <  1:378849/93‑1 (MQ=255)

           ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              >  2:2842611/1‑93 (MQ=255)

            tAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGAt                                                                             <  2:1083623/93‑1 (MQ=255)

             aGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATc                                                                            >  1:4030993/1‑93 (MQ=255)

                cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGTTACTGATAACGCTGGCGCCATTGCCCAGATCGTa                                                                         <  2:2974216/93‑1 (MQ=255)

                cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTa                                                                         <  1:1119776/93‑1 (MQ=255)

                cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGACGCCATTGCCCAGATCGTa                                                                         <  2:1966589/93‑1 (MQ=255)

                 gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                        >  1:1402387/1‑93 (MQ=255)

                 gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                        >  2:4581272/1‑93 (MQ=255)

                   cgcAATGCGGCCAAAATCGAAAGGACTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTacac                                                                      <  1:3445359/93‑1 (MQ=255)

                    gcAATGGGGCCAAAATCGAAAGGAATATGAACGATTGGCATAAAGAGGTGGGTGGTAGTAATAACGCGGGCGCCATAGTTCAGATCGTacacg                                                                     >  2:3977645/1‑92 (MQ=255)

                    gcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAcc                                                                     >  1:3745414/1‑93 (MQ=255)

                        tGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgccg                                                                 >  2:286422/1‑93 (MQ=255)

                        tGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgccg                                                                 >  2:1933933/1‑93 (MQ=255)

                          cGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCa                                                               >  2:3543066/1‑93 (MQ=255)

                           ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                              <  1:701066/93‑1 (MQ=255)

                            gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                             <  1:3398428/93‑1 (MQ=255)

                            gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                             <  1:678916/93‑1 (MQ=255)

                            gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                             <  2:1267390/93‑1 (MQ=255)

                              cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                           >  2:659884/1‑93 (MQ=255)

                              cAAAATCGAAAGGATCATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                           <  2:4396982/93‑1 (MQ=255)

                               aaaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtc                                                          >  1:897642/1‑93 (MQ=255)

                               aaaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtc                                                          >  1:4426147/1‑93 (MQ=255)

                                 aaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTcgc                                                        <  2:3098603/93‑1 (MQ=255)

                                   acGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCgg                                                      <  1:3521087/92‑1 (MQ=255)

                                   acGAAAGGATTATCAACGATCGGCAGAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCgg                                                      <  1:2745592/92‑1 (MQ=255)

                                    cGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCggg                                                     <  1:2494017/93‑1 (MQ=255)

                                     gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  1:816353/1‑93 (MQ=255)

                                     gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  1:816354/1‑93 (MQ=255)

                                     gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  1:2182472/1‑93 (MQ=255)

                                     gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  2:1340356/1‑93 (MQ=255)

                                     gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    <  2:689367/93‑1 (MQ=255)

                                     gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    <  1:2820687/93‑1 (MQ=255)

                                      aaaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGtt                                                   >  2:1063892/1‑93 (MQ=255)

                                      aaaGGATTATCAACGATCGGCATAAAGAAGTCGGGGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGtt                                                   >  1:1801133/1‑93 (MQ=255)

                                        aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 >  2:1900582/1‑93 (MQ=255)

                                        aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 >  2:4399152/1‑93 (MQ=255)

                                        aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 <  2:3363301/93‑1 (MQ=255)

                                        aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 >  1:916093/1‑93 (MQ=255)

                                        aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 <  1:3102752/93‑1 (MQ=255)

                                        aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 >  2:2292556/1‑93 (MQ=255)

                                         ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                                <  2:1943759/93‑1 (MQ=255)

                                         ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                                <  2:1727116/93‑1 (MQ=255)

                                         ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                                >  1:2180912/1‑93 (MQ=255)

                                          gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                               >  2:3543123/1‑93 (MQ=255)

                                          gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                               >  1:4391102/1‑93 (MQ=255)

                                          gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                               >  1:1337062/1‑93 (MQ=255)

                                           aTTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAtt                                              <  1:3655915/93‑1 (MQ=255)

                                           aTTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAtt                                              >  2:3086698/1‑93 (MQ=255)

                                           aTTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAtt                                              >  2:559391/1‑93 (MQ=255)

                                           aTTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGATTCAtt                                              >  1:4535359/1‑93 (MQ=255)

                                             tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                            <  2:1994726/93‑1 (MQ=255)

                                              aTCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAcc                                           >  1:3337381/1‑93 (MQ=255)

                                               tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                          >  2:3754383/1‑93 (MQ=255)

                                                cAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCca                                         >  2:4690198/1‑93 (MQ=255)

                                                 aaCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcac                                        >  2:3624548/1‑93 (MQ=255)

                                                   cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                      <  1:1002409/93‑1 (MQ=255)

                                                   cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                      <  1:1952915/93‑1 (MQ=255)

                                                    gATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAg                                     >  2:2067332/1‑93 (MQ=255)

                                                         gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                >  2:1690737/1‑93 (MQ=255)

                                                           aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                              >  1:487006/1‑93 (MQ=255)

                                                            tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                             >  2:4540403/1‑93 (MQ=255)

                                                            tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                             >  2:4540402/1‑93 (MQ=255)

                                                            tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                             >  1:1011375/1‑93 (MQ=255)

                                                             aaaGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCgg                            <  2:3909710/93‑1 (MQ=255)

                                                              aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                           >  1:4254010/1‑93 (MQ=255)

                                                               agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                          <  1:1164629/93‑1 (MQ=255)

                                                               agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                          >  2:1342812/1‑93 (MQ=255)

                                                               agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                          <  1:2487565/93‑1 (MQ=255)

                                                                gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                         <  1:3186831/93‑1 (MQ=255)

                                                                gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                         >  2:4276031/1‑93 (MQ=255)

                                                                  agTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaa                       >  1:796728/1‑93 (MQ=255)

                                                                    tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                     >  2:766462/1‑93 (MQ=255)

                                                                      ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACaaa                   <  2:2270033/93‑1 (MQ=255)

                                                                        tggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACt                 >  2:3650096/1‑93 (MQ=255)

                                                                         ggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt                >  2:2445143/1‑93 (MQ=255)

                                                                          gtACTGATAACGCTGGCGCCCTTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc               <  2:3400131/93‑1 (MQ=255)

                                                                            aCTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGc             >  1:3536123/1‑93 (MQ=255)

                                                                             cTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGcc            <  1:1412041/93‑1 (MQ=255)

                                                                             cTGATAACCCCGGCGCCCATGCCCAGATCGTACCCCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGcc            <  1:982441/93‑1 (MQ=255)

                                                                              tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt           <  2:2408272/93‑1 (MQ=255)

                                                                              tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt           <  2:4450766/93‑1 (MQ=255)

                                                                              tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt           <  2:2238586/93‑1 (MQ=255)

                                                                               gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTg          >  1:246184/1‑93 (MQ=255)

                                                                                 tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt        >  1:3629034/1‑93 (MQ=255)

                                                                                 tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt        >  2:4483096/1‑93 (MQ=255)

                                                                                 tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt        >  1:2015975/1‑93 (MQ=255)

                                                                                  aaCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTc       >  1:4409487/1‑93 (MQ=255)

                                                                                     gCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACt    <  2:2225349/93‑1 (MQ=255)

                                                                                      cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa   >  1:2186854/1‑93 (MQ=255)

                                                                                      cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa   >  2:439680/1‑93 (MQ=255)

CGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAT  >  NC_000913/1845762‑1845941

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