BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	827,018	C→T	S396N (AGT→AAT)	ybhF ←	YbhF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	827,018	0	C	T	99.6%	915.2 / ‑5.5	253	S396N (AGT→AAT)	ybhF	YbhF
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (129/123); total (129/124)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 275 total aligned reads displayed.

AAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTACCC  >  NC_000913/826928‑827107

aaaaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttc                                                                                         >  2:7594610/1‑93 (MQ=255)

aaaaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttc                                                                                         <  1:4631481/93‑1 (MQ=255)

aaaaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttc                                                                                         <  2:2418595/93‑1 (MQ=255)

 aaaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttct                                                                                        >  1:2784959/1‑93 (MQ=255)

  aaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCGTTACCGGAAttctc                                                                                       >  2:1223112/1‑93 (MQ=255)

   aaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctct                                                                                      >  2:6492254/1‑93 (MQ=255)

    aGCGTAAATTCTGTCCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCAATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCtt                                                                                     <  1:243671/93‑1 (MQ=255)

     gCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCttt                                                                                    <  1:6254503/93‑1 (MQ=255)

       gTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCa                                                                                  >  1:66260/1‑93 (MQ=255)

       gTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCa                                                                                  <  1:9147909/93‑1 (MQ=255)

       gTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCa                                                                                  <  1:9642046/93‑1 (MQ=255)

       gTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCa                                                                                  <  1:9906805/93‑1 (MQ=255)

        tAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAg                                                                                 >  1:1284168/1‑93 (MQ=255)

        tAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAg                                                                                 <  1:7535217/93‑1 (MQ=255)

         aaaTTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGa                                                                                <  2:9475544/93‑1 (MQ=255)

          aaTTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGat                                                                               <  2:3208338/93‑1 (MQ=255)

           aTTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatc                                                                              >  1:3034103/1‑93 (MQ=255)

              cTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccat                                                                           <  1:3709500/93‑1 (MQ=255)

              cTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccat                                                                           <  2:3458709/93‑1 (MQ=255)

              cTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccat                                                                           >  1:8160017/1‑93 (MQ=255)

              cTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccat                                                                           >  2:5980013/1‑93 (MQ=255)

               tGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatc                                                                          <  2:1275244/93‑1 (MQ=255)

               tGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGatccatc                                                                          <  1:4894556/93‑1 (MQ=255)

                 ttCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATccc                                                                        >  2:3864482/1‑93 (MQ=255)

                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTTTTTCAGATCCATCCCCa                                                                      <  2:2552283/93‑1 (MQ=255)

                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCa                                                                      <  2:1592899/93‑1 (MQ=255)

                     aCCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGc                                                                    <  1:7189638/93‑1 (MQ=255)

                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCAcc                                                                 >  1:9445166/1‑93 (MQ=255)

                         tCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCa                                                                >  1:9386257/1‑93 (MQ=255)

                         tCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCa                                                                >  1:1040519/1‑93 (MQ=255)

                         tCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCa                                                                >  1:1896678/1‑93 (MQ=255)

                            ggTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgcg                                                             >  2:4325138/1‑93 (MQ=255)

                             gTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTGAGATCCATCCCCAGCACCAgcgc                                                            >  2:7208726/1‑93 (MQ=255)

                             gTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgcgc                                                            >  1:5774306/1‑93 (MQ=255)

                             gTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgcgc                                                            >  1:1690748/1‑93 (MQ=255)

                              ttACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGcc                                                           <  2:4651613/93‑1 (MQ=255)

                               tACCGTAGAACGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCt                                                          >  2:8261872/1‑93 (MQ=255)

                                 ccGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTgg                                                        >  1:4432471/1‑93 (MQ=255)

                                 ccGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTgg                                                        >  2:3932444/1‑93 (MQ=255)

                                  cGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGc                                                       >  2:892043/1‑93 (MQ=255)

                                  cGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGc                                                       <  2:5242159/93‑1 (MQ=255)

                                   gTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGcc                                                      <  2:1457547/93‑1 (MQ=255)

                                    tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     <  2:9262195/93‑1 (MQ=255)

                                    tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     <  1:4115307/93‑1 (MQ=255)

                                    tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     <  1:3692600/93‑1 (MQ=255)

                                    tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     >  2:6942605/1‑93 (MQ=255)

                                    tAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     >  1:6448443/1‑93 (MQ=255)

                                    tAGAGCGAAAATTTTGGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCg                                                     <  1:2196849/93‑1 (MQ=255)

                                     agagCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCgg                                                    >  1:563132/1‑93 (MQ=255)

                                      gagCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGa                                                   >  2:5629452/1‑93 (MQ=255)

                                         cGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGt                                                <  1:1997116/93‑1 (MQ=255)

                                         cGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGt                                                <  1:3836099/93‑1 (MQ=255)

                                          gAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTc                                               <  2:2936821/93‑1 (MQ=255)

                                           aaaaTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCg                                              <  2:1118410/93‑1 (MQ=255)

                                           aaaaTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCg                                              <  1:6836275/93‑1 (MQ=255)

                                            aaaTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCgg                                             >  2:3138196/1‑93 (MQ=255)

                                             aaTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGc                                            <  1:7386019/93‑1 (MQ=255)

                                                ttttGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCAcc                                         <  2:5764557/93‑1 (MQ=255)

                                                 tttGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCa                                        <  2:4821777/93‑1 (MQ=255)

                                                 tttGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCa                                        <  2:1876427/93‑1 (MQ=255)

                                                    gcgcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCa                                     >  2:5128550/1‑93 (MQ=255)

                                                     cgcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCaa                                    <  2:6368259/93‑1 (MQ=255)

                                                      gcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCaaa                                   <  1:3502548/93‑1 (MQ=255)

                                                       ccATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAAc                                  <  1:6333475/93‑1 (MQ=255)

                                                       ccATATAGCCGAGATGCAGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAAc                                  <  1:8980488/93‑1 (MQ=255)

                                                         atatAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCg                                <  1:10035767/93‑1 (MQ=255)

                                                         atatAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCg                                >  1:8652814/1‑93 (MQ=255)

                                                         atatAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCg                                <  1:4489228/93‑1 (MQ=255)

                                                          tataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGGTCCATCCCCAGCACCAGCGGCTGGCCGGAAGTCGGCACCAGCAAACCGc                               >  1:2314304/1‑93 (MQ=255)

                                                          tataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCga                               >  2:8191016/1‑92 (MQ=255)

                                                          tataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGc                               <  1:2881725/93‑1 (MQ=255)

                                                          tataGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGc                               >  1:4913437/1‑93 (MQ=255)

                                                            taGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcac                             >  2:6154809/1‑93 (MQ=255)

                                                             aGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcaca                            <  1:7506751/93‑1 (MQ=255)

                                                             aGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGcaca                            >  1:8286072/1‑93 (MQ=255)

                                                               ccGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatc                          >  1:7955214/1‑93 (MQ=255)

                                                                cGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcc                         >  1:1823868/1‑92 (MQ=255)

                                                                cGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatca                         >  2:8591574/1‑93 (MQ=255)

                                                                 gagtTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        >  2:7687701/1‑93 (MQ=255)

                                                                 gagaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        >  1:8834513/1‑93 (MQ=255)

                                                                 gagaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        >  1:2247896/1‑93 (MQ=255)

                                                                 gagaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        >  1:3405691/1‑93 (MQ=255)

                                                                 gagaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        >  2:7369075/1‑93 (MQ=255)

                                                                  agaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACatcatc                       >  2:405456/1‑93 (MQ=255)

                                                                   gaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCTCATCATCt                      >  2:4611992/1‑93 (MQ=255)

                                                                    aTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCtt                     <  2:8483403/93‑1 (MQ=255)

                                                                      gCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTaa                   <  2:2756677/93‑1 (MQ=255)

                                                                      gCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTaa                   <  2:805456/93‑1 (MQ=255)

                                                                       cTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTaaa                  <  1:5347567/93‑1 (MQ=255)

                                                                        tGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAg                 <  2:5340179/93‑1 (MQ=255)

                                                                         ggCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAgg                <  2:9509215/93‑1 (MQ=255)

                                                                          gcgcgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggt               <  1:8459176/93‑1 (MQ=255)

                                                                          gcgcgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggt               <  2:9307071/93‑1 (MQ=255)

                                                                             cgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtggt            >  1:4951021/1‑93 (MQ=255)

                                                                             cgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtggt            >  2:1739218/1‑93 (MQ=255)

                                                                                  tACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAtt       <  2:776730/93‑1 (MQ=255)

                                                                                   aCCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAttt      >  1:7746249/1‑93 (MQ=255)

                                                                                   aCCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAttt      >  1:9753620/1‑93 (MQ=255)

                                                                                     cGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAc    >  1:6029485/1‑93 (MQ=255)

                                                                                      ggAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAcc   <  2:3106407/93‑1 (MQ=255)

AAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTACCC  >  NC_000913/826928‑827107

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: `‑`