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	98.5%	468.8 / ‑3.8	135	S396N (AGT→AAT)	ybhF	YbhF
Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (54/79); total (54/81)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 147 total aligned reads displayed.

AAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTACCC  >  NC_000913/826928‑827107

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

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

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

  aaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctc                                                                                       >  2:1900464/1‑93 (MQ=255)

  aaaGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctc                                                                                       <  2:371613/93‑1 (MQ=255)

   aaGCGTAAATTCTGTTCGACCGTCAGGCTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAAttctct                                                                                      >  2:1461204/1‑93 (MQ=255)

    aGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCtt                                                                                     <  1:171432/93‑1 (MQ=255)

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

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

      cGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTc                                                                                   >  2:4010207/1‑93 (MQ=255)

      cGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTc                                                                                   <  2:4080396/93‑1 (MQ=255)

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

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

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

          atttCTGTTTGACTGTCAGCTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGat                                                                               <  2:2822441/91‑1 (MQ=255)

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

           aTTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatc                                                                              <  1:1571057/93‑1 (MQ=255)

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

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

            ttCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcc                                                                             >  1:118454/1‑93 (MQ=255)

             tCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGTCGCGCTTTACCGGAATTCTCTTTCAGatcca                                                                            <  2:529459/93‑1 (MQ=255)

             tCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcca                                                                            <  2:3652191/93‑1 (MQ=255)

                gTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatcgatcc                                                                         >  2:3059873/1‑93 (MQ=255)

                gTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGatccatcc                                                                         >  1:230430/1‑93 (MQ=255)

                 gtCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATccc                                                                        <  2:2478126/92‑1 (MQ=255)

                  tCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATcccc                                                                       <  2:218352/93‑1 (MQ=255)

                  tCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATcccc                                                                       <  2:3584618/93‑1 (MQ=255)

                    gACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAg                                                                     >  2:3353255/1‑93 (MQ=255)

                    gACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAg                                                                     >  2:2745383/1‑93 (MQ=255)

                     aCCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGc                                                                    >  1:311852/1‑93 (MQ=255)

                     aCCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGc                                                                    >  1:2836472/1‑93 (MQ=255)

                      ccGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCa                                                                   >  1:1877403/1‑93 (MQ=255)

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

                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCAcc                                                                 <  2:2346279/93‑1 (MQ=255)

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

                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCAcc                                                                 <  1:1460763/93‑1 (MQ=255)

                          cAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAg                                                               <  2:3589631/93‑1 (MQ=255)

                          cAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAg                                                               >  1:4178389/1‑93 (MQ=255)

                           aGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAgc                                                              <  1:1000133/93‑1 (MQ=255)

                           aGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTCCAGATCCATCCCCAGCACCAgc                                                              >  1:3612444/1‑93 (MQ=255)

                           aGGTGACCGTAGAGCGAAAATTGTTGCGCCATATAGGCGAGATGCTGGCGCGCTTTACCGGAATTCTCATTCAGATCCAGCCCCAGCACCAgc                                                              >  2:3154968/1‑93 (MQ=255)

                              ttACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTAAGATCAATCCCCAGAACCAGCGcc                                                           >  2:1536414/1‑93 (MQ=255)

                                aCCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTg                                                         <  1:1711259/93‑1 (MQ=255)

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

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

                                   gTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGcc                                                      >  2:2313929/1‑93 (MQ=255)

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

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

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

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

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

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

                                           aaaaTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCg                                              >  2:2860266/1‑93 (MQ=255)

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

                                              aTTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCa                                           >  2:561839/1‑93 (MQ=255)

                                               tttttGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCAc                                          >  2:1538805/1‑93 (MQ=255)

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

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

                                                  ttGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAg                                       <  1:831094/93‑1 (MQ=255)

                                                   tGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGc                                      <  1:4017789/93‑1 (MQ=255)

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

                                                      gcCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTAAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCaaa                                   <  1:4331891/93‑1 (MQ=255)

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

                                                        cATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCCCCAGCAAAcc                                 >  1:2217719/1‑93 (MQ=255)

                                                        cATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAAcc                                 >  2:1410145/1‑93 (MQ=255)

                                                        cATATAGCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAAcc                                 <  1:4511728/93‑1 (MQ=255)

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

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

                                                              gCCGAGATGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcat                           <  2:4737956/93‑1 (MQ=255)

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

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

                                                                 gagaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCAcatcat                        <  1:1283313/93‑1 (MQ=255)

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

                                                                  agaTGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACatcatc                       <  2:2703448/93‑1 (MQ=255)

                                                                     tGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTa                    >  2:3019336/1‑93 (MQ=255)

                                                                     tGCTGGCGCGCTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTa                    <  2:3385563/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                            gcgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCCCCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggcgg             >  1:1176976/1‑93 (MQ=255)

                                                                             cgcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtggt            <  2:3428191/93‑1 (MQ=255)

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

                                                                             cgcTTTACCGGAATGCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAggtggt            <  1:153016/93‑1 (MQ=255)

                                                                              gcTTTACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTc           <  2:4416045/93‑1 (MQ=255)

                                                                                tttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGa         >  2:165961/1‑93 (MQ=255)

                                                                                tttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGa         >  2:3471847/1‑93 (MQ=255)

                                                                                 ttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAt        <  1:2232502/93‑1 (MQ=255)

                                                                                 ttACCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAt        <  1:2913297/93‑1 (MQ=255)

                                                                                 ttACCGGAATTCTCTCTCAGATCCACCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAGCCGCTCATCATCTTAAAGGTGGTCGAt        <  1:2232501/93‑1 (MQ=255)

                                                                                   aCCGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGAttt      <  1:3209483/93‑1 (MQ=255)

                                                                                     cGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAc    <  1:3286481/93‑1 (MQ=255)

                                                                                     cGGAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAc    <  1:1023191/93‑1 (MQ=255)

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

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

                                                                                      ggAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAcc   >  1:3600146/1‑93 (MQ=255)

                                                                                       gAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAccc  >  2:4711577/1‑93 (MQ=255)

                                                                                       gAATTCTCTTTCAGATCCATCCCCAGCACCAGCGCCTGGCCGGAAGTCGGCACCAGCAAACCGCACATCATCTTAAAGGTGGTCGATTTAccc  <  2:1372997/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: `‑`