BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,456,098	A→G	Q41Q (CAA→CAG)	yheG →	YheG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,456,098	0	A	G	98.8%	545.7 / ‑2.8	170	Q41Q (CAA→CAG)	yheG	YheG
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (79/89); total (79/91)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 178 total aligned reads displayed.

GGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGG  >  NC_000913/3456008‑3456187

ggCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCg                                                                                         <  1:4667046/93‑1 (MQ=255)

  cgcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTc                                                                                       <  2:4145843/93‑1 (MQ=255)

   gcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCa                                                                                      <  1:3850570/93‑1 (MQ=255)

   gcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCa                                                                                      <  1:2736349/93‑1 (MQ=255)

     tGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAgg                                                                                    >  1:244705/1‑93 (MQ=255)

     tGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAgg                                                                                    >  2:1801055/1‑93 (MQ=255)

      ggCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGt                                                                                   >  1:2171626/1‑93 (MQ=255)

       gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGtt                                                                                  <  2:4519246/93‑1 (MQ=255)

         acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTgg                                                                                <  1:4331128/93‑1 (MQ=255)

         acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTgg                                                                                >  1:4485288/1‑93 (MQ=255)

         acaACGAATTGGTTATCTCTATTCAGACGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTgc                                                                                >  1:4096315/1‑92 (MQ=255)

          caACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGc                                                                               >  1:2753217/1‑93 (MQ=255)

          caACGAATTGGTTATCTCTATTCAGAGGGCGAGATTAATTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGGTGGc                                                                               >  1:792617/1‑93 (MQ=255)

           aaCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGcc                                                                              <  2:4591310/93‑1 (MQ=255)

           aaCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGcc                                                                              >  1:4458427/1‑93 (MQ=255)

            aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCa                                                                             <  1:4274421/93‑1 (MQ=255)

             agAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAg                                                                            <  1:3927280/92‑1 (MQ=255)

              gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGt                                                                           <  2:838021/93‑1 (MQ=255)

              gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGt                                                                           <  2:1123937/93‑1 (MQ=255)

              gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGt                                                                           >  2:1702135/1‑93 (MQ=255)

                 ttGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGcc                                                                        >  1:1161171/1‑93 (MQ=255)

                 ttGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGcc                                                                        <  2:2486913/93‑1 (MQ=255)

                  tGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCa                                                                       <  2:1422744/93‑1 (MQ=255)

                   ggTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAg                                                                      >  1:1075066/1‑93 (MQ=255)

                   ggTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAg                                                                      <  2:1560039/93‑1 (MQ=255)

                        tctcTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCAt                                                                 >  2:4128717/1‑93 (MQ=255)

                         ctctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATg                                                                <  2:1657468/93‑1 (MQ=255)

                           ctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAc                                                              >  1:1064124/1‑93 (MQ=255)

                            tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             >  2:470397/1‑93 (MQ=255)

                            tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             >  2:3973036/1‑93 (MQ=255)

                            tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             <  2:2953914/93‑1 (MQ=255)

                             aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            <  1:1475188/93‑1 (MQ=255)

                             aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            >  1:4584773/1‑93 (MQ=255)

                              ttCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAg                                                           <  2:4089183/93‑1 (MQ=255)

                              tgcAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAg                                                           <  1:1857299/91‑1 (MQ=255)

                                cAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCt                                                         <  2:173640/93‑1 (MQ=255)

                                 agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        <  2:3604171/93‑1 (MQ=255)

                                 agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        >  2:1300071/1‑93 (MQ=255)

                                 agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCAATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        >  1:1657904/1‑93 (MQ=255)

                                  gagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCAAGTGCCAGACCATGACCAGCTTg                                                       >  1:993804/1‑93 (MQ=255)

                                    gggCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTc                                                     <  1:891190/93‑1 (MQ=255)

                                      gCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcac                                                   >  2:422696/1‑93 (MQ=255)

                                       ggagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGGCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcaca                                                  <  1:4002358/92‑1 (MQ=255)

                                         agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAgg                                                <  2:4389811/93‑1 (MQ=255)

                                         agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAgg                                                <  2:4244734/93‑1 (MQ=255)

                                         agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAgg                                                <  2:3367873/93‑1 (MQ=255)

                                          gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                               <  1:3249893/93‑1 (MQ=255)

                                           attattTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCt                                              <  2:4253786/93‑1 (MQ=255)

                                           attattTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCt                                              <  2:3199123/93‑1 (MQ=255)

                                            ttatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTg                                             >  2:1469890/1‑93 (MQ=255)

                                             tatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTgt                                            >  1:3548122/1‑92 (MQ=255)

                                               tttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtt                                          <  2:3185655/93‑1 (MQ=255)

                                                  atCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTg                                       <  2:2282870/93‑1 (MQ=255)

                                                   tctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                      <  1:2935313/93‑1 (MQ=255)

                                                    ctcGCCGATGCGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                     <  1:3292360/93‑1 (MQ=255)

                                                     tcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAg                                    >  1:2743675/1‑93 (MQ=255)

                                                     tcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAg                                    >  2:1533193/1‑93 (MQ=255)

                                                      cGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGCTTTTGAAGc                                   >  2:1216094/1‑93 (MQ=255)

                                                      cGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGc                                   <  2:2860723/93‑1 (MQ=255)

                                                        ccGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTc                                 >  1:4193683/1‑93 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                               <  2:4287298/93‑1 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                               >  1:3243517/1‑93 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                               >  2:2884416/1‑93 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgc                               >  1:3667999/1‑92 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGCCACAGGCTGATTTTGAAGCCCgg                               >  2:3441074/1‑93 (MQ=255)

                                                            tACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCGCAGGCTGATTTTGAAGCTCGGCt                             >  2:2616997/1‑93 (MQ=255)

                                                             aCGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTg                            <  1:635687/93‑1 (MQ=255)

                                                              cGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTgg                           >  1:4712033/1‑93 (MQ=255)

                                                                ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                         <  1:967562/93‑1 (MQ=255)

                                                                ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                         >  1:4159106/1‑93 (MQ=255)

                                                                ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                         <  2:17437/93‑1 (MQ=255)

                                                                ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                         <  1:2444821/93‑1 (MQ=255)

                                                                ccATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                         >  1:1951136/1‑93 (MQ=255)

                                                                 cATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAg                        >  2:3670048/1‑93 (MQ=255)

                                                                 cATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAg                        >  2:750881/1‑93 (MQ=255)

                                                                   ttCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                      <  1:4550294/93‑1 (MQ=255)

                                                                   ttCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                      >  1:3497883/1‑93 (MQ=255)

                                                                     cGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGt                    <  1:2636633/93‑1 (MQ=255)

                                                                     cGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGt                    <  1:1571072/93‑1 (MQ=255)

                                                                      gAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                   <  1:2143894/93‑1 (MQ=255)

                                                                       ggCGGTTACTCGATATTCAGCGTGAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTGGGCTGGAAGCGGTAt                  <  1:3786939/92‑1 (MQ=255)

                                                                       aGCGGTTACTCGATATTCAGCGTCAGGTTTGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGGAt                  >  2:2549098/1‑93 (MQ=255)

                                                                       aGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAt                  >  1:701284/1‑93 (MQ=255)

                                                                          ggTTGCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc               <  2:4035061/93‑1 (MQ=255)

                                                                          ggTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc               <  2:3386358/93‑1 (MQ=255)

                                                                          ggTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc               <  2:2946517/93‑1 (MQ=255)

                                                                           gTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCa              >  1:2014417/1‑93 (MQ=255)

                                                                           gTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCa              >  2:2244692/1‑93 (MQ=255)

                                                                              aCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa           <  2:3641519/93‑1 (MQ=255)

                                                                              aCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa           <  2:1932280/93‑1 (MQ=255)

                                                                               cTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAg          >  1:4744809/1‑93 (MQ=255)

                                                                                  gATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAAt       <  1:2535737/93‑1 (MQ=255)

                                                                                  gATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAAt       >  2:3748312/1‑93 (MQ=255)

                                                                                  cagatTCACCTTCACGTTGGCCAGTTCCAGACCATGGCCAGCTTGTCACAGGCTGATTTTGACGCTCGGCTGGAAGCGGTATTCCATAAGAAt       <  2:3681256/90‑1 (MQ=255)

                                                                                   atatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGTTCGGCTGGAAGCGGTATTCCATCAGAATa      <  2:2757813/93‑1 (MQ=255)

                                                                                   atatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATa      >  1:141910/1‑93 (MQ=255)

                                                                                   atatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATa      <  2:2757812/93‑1 (MQ=255)

                                                                                     atTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAcc    <  2:1172783/93‑1 (MQ=255)

                                                                                     atTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAcc    >  2:4621124/1‑93 (MQ=255)

                                                                                       tCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCgg  <  1:1489727/93‑1 (MQ=255)

GGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGG  >  NC_000913/3456008‑3456187

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